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Abstract. The final period of decay of magneto-turbulence in an external, homo- 
geneous magnetic field is considered and it is shown that it develops pronounced axi- 
symmetric properties, turbulence elements with finite wave numbers in the direction 
of the field being damped strongly under normal physical conditions. The turbulence 
consists of aperiodic motions as well as wave motions. An introduction of an angular 
velocity, inclined to the field, destroys the axisymmetry and modifies the damping 
effects and periodicity. The influence of the magnetic field on the damping is counter- 
acted by the Coriolis force. A linear stationary theory on the action of the field gives 
results consistent with those of the theory of decay. From the results of both theories 
an explanation is given of the observed inhibition of turbulence in mercury by a magnetic 
field. 

I. Introduction. From the large linear dimensions which one encounters in cosmical 
physics one may infer turbulence as basic for many applications in this field. In most 
applications the electrical conductivity is high enough for electrodynamic forces to 
play an essential role and for turbulence to be governed by the laws of magneto-hydro- 
dynamics. However, it has been pointed out by Chandrasekhar’ that even a small 
angular velocity of a medium of cosmical dimensions may sometimes be associated with 
a Coriolis force of the same importance as the electrodynamic force. In the general case 
a systematic magnetic field, created by constant sources, has also to be taken into account 
and turbulent fluctuations have to be superposed on the field. 

Even if the influence of the boundaries is neglected, the general turbulent state is 
neither homogeneous, nor isotropic. The variation of mean velocity with position due 
to rotation destroys homogeneity also when the mean magnetic field is homogeneous. 
In this paper, however, we shall discuss turbulence in an incompressible liquid, where 
the centrifugal potential plays no essential role and regions are considered which are 
sufficiently small to justify the assumption of a homogeneous external field. 

The problem of homogeneous, isotropic magneto-turbulence has been treated by 
Batchelor’, Chandrasekhar® and Lundquist* among others. The second of these papers 
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gives a treatment in terms of invariant theory and the third presents a decay law for 
the spectral tensors for large wave numbers. 

The purpose of this paper is mainly to discuss how the law of decay of homogeneous 
magneto-turbulence is modified by the introduction of an external magnetic field of 
the strength B and of the Coriolis force due to a constant angular velocity Q (MKSA- 
units are used in the following). Small amplitudes will be assumed and triple correlations 
will be neglected. 

In Sec. III the law of decay in a homogeneous magnetic field is derived and it is 
shown that the motion becomes axisymmetric with respect to the direction of the field. 
In Sec. IV the analogous treatment in the presence of a constant angular velocity in a 
direction making an angle with the magnetic field is given. In Sec. V a brief discussion 
of’stationary turbulence is given. The results are used in Sec. VI for an interpretation 
of the experimental results of Hartmann’’® and Lehnert’’* on turbulence in mercury. 

II. The fundamental equations. The velocity field v, in an incompressible liquid 
with constant electrical conductivity o, kinematic viscosity v, absolute permeability 
uw and density p, is supposed to be nonrelativistic and the liquid is assumed to rotate 
with a constant angular velocity Q in a homogeneous magnetic field B = «H, making 
an angle with Q and produced by external sources. The conductivity is assumed to be so 
good and the rate of change with time so slow that the displacement current can be 
neglected compared with the convection current j, which is the source of the induced 
magneto-motive force h. We start with the equations 


6 


curlh =j,  curlE = «pe (1) 
Ot 

divh = div v = 0, (2) 

j = o[E+ uv X (H + b)] (3) 


and 


; : 
| 2 +(v-V)v+2a2xv+aQx(Qx 2) | 
oO 


= pj X (H +h) + »pV'v — Vip + 1d,), (4) 


where E is the electric field, p the pressure, x = (x, , 22 , £3) the radius vector from an 
origin on the axis of rotation and ¢, the gravitation potential. The reduction of the system 
into two equations is similar to the treatment by Walén” and Lundquist* and will not 
be given here in detail. If the z-axis is chosen in the direction of the angular velocity 
Q, the equations become 


av = (W-V)v + AV'V + curl (v X V) 8 


ot 
and 


or =(W-V)V+7Vv 4+ 2v X¥ Q-— Vo—-—VX ewl V+ Vv X curly, (6) 





5Kgl. Danske Vidensk. Selskab Math.-fys. Medd. 15, No. 6 (1937). 
6J. Hartmann and F. Lazarus, Ibid. 15, No. 7 (1937). 

7Arkiv f. Fysik 5, 69 (1952). 

®Tellus 4, 63 (1952). 
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where 
V=h(u/p)'?, W=H(u/e)'’, A= 1/(uo), (7) 
= p/p+o,+¢. + WV + w'/2 (8) 
and 
b. = 30x; + 23) (9) 


is the centrifugal potential. 
For small amplitudes the terms of second order in Eqs. (5), (6) and (8) may be 


neglected. Further, letting 


J = curl V, w» = curl v, (10) 
and taking the curl of Eqs. (5) and (6) we obtain 
J _ (W-V)o + VJ (11) 
and 
22 = (W-V)J + V%o + 20 X O+ V9, (12) 


where, from well-known vector operations, 
y = 2v-Q. (13) 


III. The law of decay in an external magnetic field. The decay of magneto-turbu- 
lence in an external magnetic field is of special interest in connexion with the experi- 
mental investigations on mercury. In an exact theory we must introduce triple correla- 
tions which represent the basic mechanism of turbulent interaction. However, it is 
known that they do not play a role in the final period and in this paper we shall restrict 
ourselves to this stage. 

1. Discussion of the gradient term. We start with Eqs. (5) and (6) for a liquid, con- 
tained in a finite region and suppose that outside this region the conductivity and 
viscosity are zero while the permeability has the same value as in the inside. Since v and 
V are solenoidal vectors, we obtain on taking the divergence of Eq. (6) that 


if all second order terms are neglected. Outside the liquid we have no terms representing 
the electromagnetic and the viscous forces and the result (14) continues to hold. Further, 
no surface currents are allowed to exist due to the finite value of the electrical conduc- 
tivity. Thus, @ given by the expression (8) is continuous throughout the bounding 
surface. In the external region V tends to a dipole field at large distances from the liquid 
and # tends to a constant value, representing the balance between pressure and gravi- 
tation force. Since ¢ is constant at infinity it is also constant on the boundary and in 
the interior of the liquid due to Green’s theorm and 


Vo =0 (15) 


all over the liquid. 
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2. The correlation tensors. In terms of components Eqs. (5) and (6) can be written 


in the forms 


aV , Ov; — ' 
oes a W, Tt + NWT, (16) 
Ot OX; 
and 
Ov . OV " 
— = na “s) 7) 
ry Wig, TV» (17 


where use has been made of Eq. (15). Summation over repeated indices is to be under- 


stood. Let 
x =x-+Fr (18) 


for a point at the distance r from x. We shall distinguish the values of the various quanti- 
ties at x’ by a prime. Equations for the correlation tensors may now be formed in the 
usual manner; thus multiplying Eq. (16) by V! ; we get, 
-, OV; > oop Os 
vi7-a = FLY; — 


AViV’V (19) 
ot dz, TNIV 


and adding to this equation that obtained from Eq. (19) by interchanging 7 and j and 
the primed and the unprimed quantities, we get 
0 , dd r 0 fare. % 0 yr fy r 2 s2\ (17 ry ‘ 
—(V,Vi) = W,| — (V9) + a (Vw) | + AV + OV )V 7. V3) (20) 
ot F Ox, z OX}. ; 
In Eq. (20) (and similar equations in the sequel) one of the vectors x and x’ is kept 
constant while the other is varied; then 
8/axt = 8/ar, = —d/dx, 


and Eq. (20) can be written in the form 
| 


£ VV!) = Wi (Vol — Vo) + 2V(V.V%). (22 


In the sequel all differential operators will refer to the variable r. 


Further equations are obtained in a similar manner by multiplying Eqs. (16) and 


(17) by v; . Introduce the tensors 

M (r,t) = (V(x, 0OVi(x +4, 0), 23) 
representing the magnetic energy, 

K;(r, 2 = 40.(x, Ov(x + 4, 0), 24) 


the kinetic energy, an 
£2, ) = (V(x, ox + 4, 0), (25 


the interaction energy respectively. In forming these tensors the mean values are with 
respect to time. In homogeneous turbulence these tensors should be invariant to arbi- 
trary displacements; this invariance leads to the following geometrical properties: 


M(t) = {V.(x’ — r)V,(x’)) = M;,(—-D), (26) 


K,,@) = K;-r) 
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and 

L;(—r) = 2(V,(x)o,x — r)) = (V(x + rv(x)) = L;,(%’, x). (28) 
Further, M,;; and K,; are symmetric, whereas L,;; is skewsymmetric. M,;(0) and K;,(0) 
are the mean magnetic and kinetic energies at the point x and time ¢. The physical 
significance of L,; is easily seen from the electrodynamic force in Eq. (17). The decrease 
in magnetic energy per unit mass and time due to the work of the electrodynamic force is 


p( MH) __ OV; 
E\Y) = W, 3a, 8 (29) 


which is a positive quantity if magnetic energy is converted into kinetic energy. From 
Eqs. (21) and (28) we get the expression 


(MH) __ 1: ] AC ov; ; a7) = {- 7 <2 — = } 
Ey; = lim s W,\e i + »; az’) = lim We oy, [L.(4) — L;(—r)]?. (30) 


Returning to Eq. (22) and taking the mean value we get 





r-0 = 


2 yf, = Wi = (Lids) — L,{—9)] + 2V?Ma; « (31) 
ot Or, 
Similarly we obtain 
o Ki, = —W, = [Lelt) — LD] + 29°" Kes (32) 
Ot Or, 
and 
° L.A = Wi Me +E +0 400M. (33) 
at ar, 


Interchanging r and —r and 7 and 7 in Eq. (33) we obtain 


2 Lid) = —Wi et (Ma — Ku) + + 9V7L(-0). G4) 


Now, assume that a spectral tensor, A,;(«), exists such that 


Aj;(%) = (8x°)~ If M,,(r) exp (—ix-r) dr, dr, drs (35) 


and 


M .,(r) 


II 


[[[ A.c00 exp (inet) dis deg ds, (36) 


where « is the wave number. A spectral kinetic energy tensor, Q;;(x), and a spectral 
interaction tensor, T,;;(%), can be similarly defined in terms of K,;(r) and Z,;(r) re- 
spectively. 

Since differentiation with respect to r, in ordinary space corresponds to multiplica- 
tion with a factor zx, in the wave-number space we have from Eqs. (31), (32), (33) and (34): 


9 2 r 
(2 + are’) As — ix,W,IT;(x) — T;,(—%)] = 0, 
Cc 


(2 =f 2m?) + ix,W,[T,,(x) — T;(—*«)] = 0, (37) 


p 6 
E +A + ie |erate — T;(—«)] — 2iK,W,(A,; — 2:3) = 0. 
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These equations are also valid for arbitrary amplitudes if sufficiently high wave numbers 


are considered. 
3. The law of decay. The system of Eqs. (37) will be satisfied by solutions of the 


form exp (mt) if the determinant of the system 


| m + 2a 0 —iF | 
| 
0 m + 2b iF \, (38) 
—2iF 2iF m+t+a+t+b| 
where 
a=)’ b=’; F = x,W,, (39) 
vanishes. Hence 
(m + a+ b)[m? + 2(a + b)m + 4ab + 4F*] = 0. (40) 
The roots of this equation are 
m2 = —(a+ b) + [(a — Bb)? — 4F*]"”, m, = —(a+ bd) (41) 
and the solutions have the form 
| Ay (x, 6) ) ac? | Tyg A® | | exp (m,t) | 
Q(x, 2) | = | 25; “oy Siri | exp (met) |, (42) 
II; ;(x, ¢) Ey II; ; IIs; | | exp (mgt) | 
where 
II,;(x, 2) = T,,(x, 0 — T;(-—«, dO (43) 


has been introduced. When m has been chosen as a root of Eq. (40) the solution of the 
linear system (37) can be written in the form 


g-a-O’yi1¢ta-?”’) 11+a-"Y7%y/n-a-)7] l 


—2if1 -(1-— yt —2i[1 + (1 — 9)"")/¢ — 20/5 | 


2}; exp (m,t) | 


257 exp (met) |, (44) 


2{? exp (mst) 


where we have introduced the parameter 


¢ = 2F/(a — b) = 2,.W,/[K(A — »)] (45) 
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and have written 
M2+2a=a—b+(a— BI — 2)”, m, + 2a=a-—b (46) 


and 


—(a—b). (4%) 


m.+2b= —-(a—-bt(a—db(1- 3)”, m;, + 2b 


For small values of ¢ the form (44) reduces to 


Q5;’ exp (mf) 
| |e 4/2 1 
| = | - | Of? exp (mt) |. (48) 
| | —7¢ —4i/¢ —2i/¢ 
| 2 exp (m,) 


All quantities have to be finite at ¢ = 0, and therefore, when ¢ tends to zero, 


Ai; = 0(¢) > 0, Ay? = ASP +079, AZ? = 0(f) >0, (49) 
Of = aA +00, a? = o(f’)-0, 2; = oo) 0, (50) 
Ii; = o(¢) > 0, Ii; = o() 3 0, ni; = mi) + od), (51) 


where A‘*’ , QS and IIS are independent of ¢ and 0(¢") are terms at least of order n. 
For ¢ = 0 we get 


A,;(w, 2) = AS?(%) exp (—2d«’2), (52) 

Q,;(x, t) = Q{)(«) exp (—2vx*t) (53) 
and 

TI, ;(x, 2) = I1s}’(%) exp [—(A + »)x°E]. (54) 


These solutions represent the final period of decay of isotropic turbulence and are 
analogous to the case studied by Lundquist where the magneto-hydrodynamic inter- 
action, expressed by triple correlations, was neglected and the magnetic and kinetic 
turbulence fields were found to decay independently of each other. The introduction 
of an external magnetic field, however, changes the situation even in first order and the 
decay is mainly governed by a coupling of the form (29). 

The solutions (41) define a decay time, 7, , given by 


W/ne=KAtY) FRA-—Y1-— PY; 1/3 =K°A+>). (55) 


Thus, for ¢? < 1 we have three non-periodic solutions, all with different decay times; 
for ¢° = 1 all the spectral tensors decay with the same time constant 


tT = 73 = V[Kr’A+ 9] (56) 


and finally for values of ¢* > 1 two periodic solutions are obtained, both with the same 
real damping, 1/7; . 

4. Physical interpretation of the law of decay. The physical meaning underlying the 
results (42), (44) and (55) can be understood from a consideration of a plane state of 
motion in a liquid between two infinitely conducting planes at a distance L, as shown 
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by Fig. 1. Introduce a homogeneous magnetic field By in the z-direction, perpendicular 
to the planes. In a plane state of motion and for small amplitudes 





a/ax = a/ay=0; v=(0,x,0; V=(,V,0 (57) 
O-co 8, 
IMB, 
Pas 
+X 
J — 





=0 





MMM fy ZZ oo 


Fic. 1. Plane state of motion of an electrically conducting, viscous liquid between two infinitely con- 
ducting planes in a homogeneous, perpendicular magnetic field, Bo . 


and Eqs. (5) and (6) reduce to 


OV , OV o7V S 

ot 4 dz a ia az" (58) 
and 

Ov OV av as 

a eee A —ze DY 

ot M Oz on Oz (59) 


On eliminating V we obtain the following equation: 


Separating the variables in the form 
v = Z(z)T(t) (61) 
we get 
T"/T — WZ"/Z — (’ + (Z"/Z)(T"/T) + WZ" /Z = 0. (62) 


Since the magnetic field lines and the liquid are attached to the infinitely conducting 
walls we may assume that 


Z = sin [xkz/L] = sin (xz) (k = 1,2, --- « = wk/L) (63) 
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in which case 7 admits a solution of the form 
T = exp (}mt), (64) 
provided 
m? + 4x°W? + 2(d + v)x?m + 4r* = 0. (65) 
The roots of this equation are 
mg = —(A\ +r) + (A — "(1 — 2°), (66) 


which are the two first values given by Eq. (41). This corresponds to a motion of the 
magnetic lines of force, regarded as elastic strings with a tension given by the magnetic 
field strength and a damping due to the Joule heat and the viscous losses. For ¢” < 1 the 
equivalent strings move aperiodically and for ¢* > 1 damped waves travel along the 
strongs. Correlations may be formed by products of the solutions of Eq. (60), giving 


time factors of the form 
TY pap = exp [3(my 2 + mM, >) t] (67) 
for a given value of «x. These factors are consistent with those given by Eqs. (41) and 


(42). 
We shall now return to the discussion of the results (42). Whatever distribution we 
may start with at ¢ = 0 these results will always tend to a solution, which for the kinetic 


tensor reduces to 
Q,;(%, t) = 25;'(x) exp (—t/7,), (68) 
where 1/r, is the value of the one of the decay factors (55) having the smallest real part. 


The solution (68) will differ for turbulence elements due to their size, the properties of 


the liquid and the strength of the external field. 
For small values of ¢ and for large ratios \/v (A/v = 10° for experiments with mercury) 


the decay factors (55) become 


Ut 


Qu? + aK°7(A — vv? (69) 


1/7; 


and 
1/r. = 2x” — 4K°(A — ve’, (70) 


whereas 1/r, remains unaltered; further, 1/7, is much smaller than 1/7, and 1/r; , if ¢ 
is sufficiently small. Now, Eq. (50) shows that the corresponding factor, Q{}’ , has a 
term of zero order in ¢ and consequently ©;;(«, ¢) will be represented by the asymptotic 


law 
Q,;(x, 2) = Qf? («) exp {—2[vn? + 3W?/K7(A — v) JE} (71) 


"7 


during the largest part of the final period of decay, if the x5-axis is chosen in the direction 
of the field and Eq. (45) is used. Analogous discussions may be carried out when A < ». 
The results of this section may be summed up in the following statements: 

During the turbulent decay of a liquid with \ > v all periodic turbulence elements as well 
as the aperiodic ones with small extensions in the direction of the field (large values of x3) 
are damped out relatively rapidly. The asymptotic state of decay is two-dimensional with 
respect to the direction of the field. Only vortices with x; = 0 are left in the final state; 
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the damping is by viscosity only and since the electric field becomes independent of 


2, no induced currents flow. 
The main kinetic decay factor, 1/7; , for mercury, given by Eq. (71) is shown in 


Fig. 2, where we have put x; = «, . The curves clearly show action of the field in sup- 
pressing elements with finite wave numbers in the X,-direction. 
nem") V(r MAY (mm 
4, 
4200 = 80061257 


er 7000063] 
4000 6000167] 


(— 500 [2] 
' —__ 400125) 





ne AB BB 500%" 


gn 
ha 010 [0013] 





f /s0 MAK ie 
120+ |\\ WRN ores 
\ \ \ ae: ™~ 
\ \\ SAO 03010030 
\ \ | e 
St toate te ~~ 040100531 
a (8) A A 


a \ \ \@s0l0065] 





800 250} ,; \ \ \ 
900 1000 1200 4.0 2015 1.0 080.7 06 
“u(m") a5) [01 wemss) 
n,Wis) CB(Vs/rf)] 
Fig. 2. The factor of decay, 1/7: , for the final period of turbulence in mercury as a function of an 
external magnetic field [W = Bo(up)~* being the corresponding wave velocity] and the wave numbers 


ky = x, and x;. The values \ = 0.763 m?/s and »y = 0.116 X 10‘ m?/s at 18°C have been used. 

IV. The decay in a magnetic field in presence of the Coriolis force. It was pointed 
out at the beginning of this paper that the influence of the Coriolis force may sometimes 
be comparable to the electromagnetic force, both leading to first order effects which 
may be significant for small as well as for large turbulent intensities. The effect of these 
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forces, when they act separately on a liquid, shows many similarities. The stationary 
flow of a rotating liquid becomes two-dimensional at high angular velocities’® and the 
same situation is true in a conducting liquid in a strong, homogeneous magnetic field.’** 
Similarly both have inhibiting effects on the onset of convection."’*’’ But if both forces 
are simultaneously present a complicated situation arises, which is not the same as the 
superposition of the separate effects. This is evident for example’ from Chandrasekhar’s 
investigations of the thermal instability of a fluid layer heated below under the joint 
effects of a Coriolis acceleration and a magnetic field. 

1. The correlation tensors. We shall now extend the theory of the decay of turbulence 
described in the preceeding sections to allow for the effect of a stationary angular velocity 
in the x,-direction: thus 


Q = (0, 0, 2) = 4 curl U; U=aQ~xx, (72) 


where U is the velocity of rotation of the liquid and x is the radius vector from the origin 
on the axis of rotation. We choose the coordinate system with the magnetic field vector 
in the x, , 2;-plane so that 


W = (0, W., Ws). (73) 

As we shall see it is most convenient to discuss this problem in terms of the current 

density and the vorticity. We now start with Eqs. (11) and (12), which written in terms 
of components have the forms 


OSs _ yy, de 4 y¥92 
at We ag, FAV (74) 
and 
Ow; , Od; oy 
_—* —W, ny + vV%0a, + 2€:1n: Qn + a (75) 


where ¢;;,, is the usual alternating symbol. The equations governing the various corre- 
lations can be obtained in the same manner as in Sec. III; the only difference is that now 
we shall have terms containing Q added to the right hand sides of the various equations; 
thus 


< (ww) _= 7e + 2w/€ 51m: Ln + 2s 5€; 1m} Ln + wi =< + @; 2. (76) 
é ay’ 
2 (Feel) = 28 + Ween e + I. - (77) 
Introduce the tensors 
R,,(x, x’) = (x) J (x +) = R;.(x’, d), (78) 
T(x, x’) = (w,(x)w,(x + 2)) = T;,(x’, x), (79) 
S,,(x, x’) = (J(x)w,(x +n), (80) 
a a 
P, (x, x’) = (2 [J (x + ¥2)1) - (2 [J (x) ¥(x + ®l) (81) 


1G. I. Taylor, Proc. Roy. Soc. Al00, 114 (1921). 
118, Chandrasekhar, Phil. Mag. 43, 501 (1952). 
128, Chandrasekhar Proc. Roy. Soc. A217, 306 (1953). 
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and 


ie = , . 3 ' 
Oi, x) = (; w(x + va); Qii(x, x’) = Q&(x, x’) — Q#@’,x), (82) 
where relations (18) and (21) have been used. These notations will be used for the present. 
They are valid also in the inhomogeneous case. The condition for homogeneity will be 
applied later. From Eqs. (76) and (77) we now find the equations governing the corre- 


lations: 
] > , r CY] rQ , Y fal 2p , >« 
ay eai(®, X = W, *s S;;(x, x’) — S;(x’, x)] + 2AV R;,(x, x’), (83) 
OT; 
0 ry f ; d 0 Y , 4) ; ry , , 
FY’ 7 ,(x, x!) = —W, a [S,,(x, x’) — S;,(x’, x)] + 2°V°T;;(x, x’) 
O77; 
+ 20, [€:imlr;(X, X!) + e102 5(x, x’)] — P;,(x, x’) (84) 
and 
0 ‘al f , r a) ‘Dp *\ 4 Al f ; j 2a f , 
ay Silk, X) = W, “- [R, (x, x’) — T;,;(x, x’)] + A + 0V'S,;(x, x’) 
o7 
+ 22 n€;:imS:(X, X’) + Q(x’, x). (85 


Interchanging 7 and j and x and x’ in Eq. (85) and subtracting the resulting equation 
from Eq. (85), we obtain 
os ; ea paca 


ap Sis, &”) — Sil’, X)] 


a ee Pome org ’ ' 
= 2W, ay [R, (x, x’) — 7;,(x,x)] + A+0V [Si(x, x’) — S;(x’, x)] 
OT; 


+ 20 n[€;1mSi(X, X’) — €:rmS;1(X, x’)] — Q;;(x, x’). (86) 


2. The spectral tensors and the law of decay. Before the equations governing the 


a 


correlations are transformed into a spectral representation we shall consider the properties 
of the tensors (81) and (82). It is easily seen that 


P;; = Q:; = 0, (87) 
since J and w are solenoidal. If homogeneity is assumed, the tensor Q(x, x’) will have 


the property 


OF(, x’) = (2 [w(x + )¥2)1) = (2 (w(x) Y(x — D1) 


‘ 


= -(2 [w,(x) W(x + pI) = —Q4(x’, x). (88) 


ar; 
From Eqs. (81) and (82) we get the result 


Q;;(x, x’) — Q;:(x, x’) = P,,;(x, x’) — P;,(x, x’) = 0. (89) 
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Now introduce the Fourier transforms ®,;(%), V,;(«), Ts;(%), x<;(%), 04 (x) and d;;(x) 
of the tensors R,,(r), 7';;(r), S.;(r), Pi;(r), Q*(r) and Q;,(r) respectively; thus 


®,,(x) = (82°)™" [|| R,,(r) exp (—ix-r) dr, dr, dr; , etc. (90) 


Also we shall write 
I; = 1T,,(-—x*); or’ = 3 *(—x). (91) 


3y applying the Fourier transforms to Eqs. (83), (84), (85) and (86) we obtain the 


following equations for the spectral tensors: 


a = 3 
(2 + 2a), — iF(T,; — 1%) = 0, (92) 
Cc 
(2 + 2b bh + F(T; — Th) — 223%; + 63%) + xi; = 0, (93) 
-d — @ ob b l’;; al iF (®;; a= WV;;) om 2; 13T 51 —_ or! = 0, (94) 
\ot 
re] , 
(2. +a+ br, + iF(®;; — V,;) — 2Me;,,T, — d* = 0 (95) 
und 
(<. ~~ ¢ > b) I’; —_ | - 21K (®;; = V;;) _ 22(E;:3T = €:131'42) + v5; = 0, (96) 


where the abbreviations (39) have been introduced. Equation (95) has been obtained 
from Eq. (85) by interchange of 7 and 7 and x and x’. 

We shall now seek solutions of Eqs. (92) to (96) which have an exponential dependence 
on time of the form exp (mt). As before we shall obtain a characteristic equation for m. 
For the diagonal terms (¢ = 7) the four first equations give 


(m + 2a)%,;; — iF(T;; — T,) = 0, (97) 
m + 2b)W,;; + iF(T;; — Ti) — 29.3%; + Vi) = 0, (98) 
(no summation) 
(m a+ b)T;; — 7F(®;; — V;;) — 29.31 ;, = 0, (99) 
(m+ a+ DT, + iF(®;; — Vii) — 2Qe.3Th, = 0, (100) 


where use has been made of the relation (87). Interchanging 7 and 7 in Eq. (92) and 
subtracting the resulting equation from Eq. (92) we obtain 


(m + 2a)(%,; — %,,) — iF(T;; — T); — Ty; + Ti) = 0, @ # J). (101) 
Similarly for Eqs. (93) and (96) we obtain 
m+ 2b)\(%;; — Vj.) + tF(Ty; — Ti, — Ty, + Ti) = 0 (102) 
and 
(m+ a-+t b\(T.; — The — Tye + Th) — QF ((® — 9) — Wy — Yd] 
+ 2%!';,(T,;; + Ti, + T;;, + Ti) = 0, (¢ + j; no summation) (103) 
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where e/;; refers to off diagonal terms only and ¢{;, = — e/;; . Further, from the equa- 
tions of the diagonal terms we derive 
(m+ a+ OT. +I + Ty, +) — 20e(T, — T%, — Ti + 1) = 


(¢ # J; nosummation) (104) 


The system of Eqs. (101), (102), (103) and (104) will be satisfied if the determinant 


m + 2a 0 —iF 0 
0 m + 2b iF 0 | 
(105) 
—2iF 2iF mt+at+ob 2; 5s | 
0 0 —20! 5 m+atd| 
vanishes. Hence 
(m + a + b)*[(m + 2a)(m + 2b) + 4F*] + 407€/7,(m + 2a)(m + 2b) = O. (106) 
Letting 
M=m+a+b and c=a-— b}, (107) 
Eq. (106) can be written in the form 
M* — M°*(c? — 4F’ — 407%e%;,) — 407%e!;.c> = 0. (108) 
The possible values of m are therefore 
m, = —(a+ b) + (2) {c? — 4F”? — 407%), + [(C — 4F’)’ 
+ 807€/7,(c? + 4F? + 207e/7,)]'/7}'” (k = 4, 5, 6, 7). (109) 


3. Discussion of the law of decay. It is seen from ia (97), (98), (99), (100) and 
the result (109) that the time dependence of the spectral tensors with at least one index 
in the direction of the angular velocity vector is uninfluenced by the Coriolis force and 
the solutions have the same form as given by Eqs. (41) and (42). However, if 7 and j 
differ from 3 the results of Sec. III are modified, each tensor splitting up into four terms: 


@,,(x, 2) = >> &” (x) exp (m,t) (kK = 4, 5, 6, 7;¢ = 3,7 € 3). (110) 


We have corresponding solutions for the wont ensors. The behaviour of the solu- 
tions in the limiting cases W = 0,2 = 0 and W = (0 will not be discussed in detail 


here. 
We shall now discuss the manner in which the angular velocity modifies the spectral 


tensors when 7 ¥ 3,7 ¥ 3 and ¢/;; = 1. Let 
€ = WAKoW. + x3Ws)/[K (A — v)] (111) 


and 


meee 


n = 22/[K(\ — »v)]; (112) 


¢ and 7 are associated with the forms W,.L./(A — v) and U,.L./(A — v) respectively, 
where W, , L, and U, refer to characteristic wave velocities, lengths and rotation veloci- 
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ties. They may also be connected with parameters used in earlier works.*'’''* For small 
values of 2 the factors of decay become 


ms = —at+b)+ecfl — Pl — 7/1 -— Y)}}'” (113) 
and 


Me.2 = —(a + bd) + ten/(1 — 2°)”, (114) 


provided n°(1 + ¢°)/(1 — ¢&) « 1. From a comparison of Eqs. (109) and (41) we con- 


clude that 


Ms.5—> M2 and Ms.7 —> Mz when 2— 0. (115) 


Thus, if the aperiodic solutions (41) are perturbed by a small angular velocity, the 
damping effect of the magnetic field in the factors, m, and m, will be counteracted 
while the term in m, will be split up into two periodic solutions. For sufficiently small 
values of 2 and F the spectral vorticity tensor will have the asymptotic behaviour 


ll 


V'?(x) exp {—2[vn? + (K2We + x3W3)7[1 — 427/x*(A — v)*]/0°(A — v) JE} 
(t # 3,j ¥ 3); (116) 


WV, ;(x, 2) 


this corresponds to the earlier law (71). 
The occurence of periodicity can be discussed in terms of M, given by Eqs. (107) 


and (109). We have 
2M =c(l— Pf — afl (1+ 40/1 — 2 — 9°)]'"}. (117) 
If ¢° + »° < 1 the solutions corresponding to the roots m,,, are aperiodic while the 


solutions corresponding to m,,; are periodic. On the other hand if ¢’ + 7’ > 1 the roots 
m,., lead to periodic solutions while m,,; lead to aperiodic solutions. The critical case 


is given by 
P+, =1, (118) 
when m, has the values 

ms} = —(a + b) & (2Me)"”’, ms? = —(a + db) + i(2M)"”. (119) 

The result (109) may be confirmed by a discussion similar to that in Sec. III, 4. 

Thus, reconsidering the situation indicated in Fig. 1 we shall now assume that the system 

is partaking in rotation and for simplicity suppose that the axis of rotation is in the 

z-direction. Further we shall assume a homogeneous state and that it is possible to 
introduce a local rectangular system of coordinates in which the approximations 

d/dx = 0, 0/dy = 0 (120) 

are valid. Equations (5) and (6) for the z- and y-components of small disturbances, 


V and v, become 


d °\, a 
ae Nae a r r ‘ai 7 VU ° : 9 
(2 r aa :,V,=W ae (v, , vy) (121) 
and 
3 ca . 
(= ait (v, ,v,) = W 7 (V., V,) + 220, , —2,). (122) 





8S. Chandrasekhar, Proc. Roy. Soc. A216, 293 (1953). 
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If V, and V, are eliminated we get 


a a (2 ) _ | ‘ (2 *) - 
—_— — ——— —_— —2 a —Zz } ) = 2¢ —_—_ — -——z —w_). 12% 
(2 . 3) at dz” " dz” W. » %) N51 ae (v, , —v.) (123) 


Separating the variables in the manner 


v; = v;°” sin (xz) exp (nd), xk = ak/L (k = 1,2, ---), (124) 
we find that with the abbreviations (39) we must satisfy the equation 
[(n + a)(n + b) + F’], , v,) = 2Q(n + adv, , —2,). (125) 
The condition that Eq. (125) allows non-zero solutions is 
[((n + a)(n + b) + F*]? + 407%(n + a)’ = O. (126) 
In other words 
n = —3[(a + b) + 270] + F{[(a — b) + BQ)’ — 4F*}"”. (127) 


° ° . . I . (II ° : 
A product of two of these solutions gives exponents, n"’ + n°’, which are easily seen 


to be of the form (109); thus 


nm? +n? = —(a+ db) + BH(c + 219)’ — 4F*}'? + Bie — 210)? — 4F"}'” 

a wigt Be Ww. (128) 
where 

N’ = 3° — 40° — 4F’) + 4[C — 40° — 4F’)’ + 160°c*]"”; (129) 


in other words the same values of M’ as given by Eq. (117)."* 

The results of this section may be summed up as follows. If an angular velocity, Q, 
is introduced into the turbulent state of motion of an electrically conducting liquid in a 
magnetic field the factors of decay of the tensor components perpendicular to Q will 
be split up into four solutions. The passage from periodic to aperiodic motions is displaced 
from ¢’ = 1 to¢’ < 1, i.e., to weaker magnetic fields. For small values of 2, the damping 
effect of the magnetic field in the aperiodic solutions is counteracted by the angular velocity. 
The remaining tensor components are uninfluenced by the angular velocity. 

4. The importance of the Coriolis force. From the foregoing discussion it is seen that 
the effect of the Coriolis force will be of the same order of magnitude as the electro- 
magnetic effects if ¢ and » are of the same order, or from the expressions (111) and (112) if 


W = B/(up)'” = QL/2r, (130) 


where L represents the linear dimensions of the disturbance. For the sun typical values 
are W = 4m/s and 2 = 2.7 X 10°° s"’; it is seen that for these values the Coriolis force 
plays an essential role when the linear dimensions of the disturbance become larger than 
L = 6 X 10° m. Since the solar radius is about 7 X 10° m this will apply to a wide range 
of possible disturbances. It will probably be even more pronounced for stars with larger 
angular velocities. 

In experiments on turbulence on laboratory scale, the largest turbulence elements 

“The same solutions are also obtained in the non-dissipative state from results by Chandrasekhar, 
Astrophys. J. 119, 7 (1954). 
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will hardly exceed 10~° m; and in mercury a common value of the wave velocity is about 
5 m/s and of the angular velocity about 10 s~*. Under these circumstances the electro- 
magnetic effect will be about 50 times larger than that of the Coriolis force. 

The influence of the Coriolis force on the critical damping in liquid sodium, as given 


by 
(x,W,)? + 2 = «(A — v)?/4 (131) 
is shown in Fig. 3. 
W (m/s) [0175] [0351070] VWs”) - 
[BVsi20 40°50 Miso nm) 
OAy(emy 





wer) > 








Q 1s) 





{iw fem] 1700 80 60 40 20 0 50 100 150 200 250 
{1] [12570167] [25] [5] 

Fic. 3. Critical values of the wave velocity, W, and the angular velocity, Q, in liquid sodium for turbu- 

lence elements with velocity components perpendicular to Q. The total wave number is ® and xy is the 

component in the direction of the magnetic field, B, where W = B(up)*. \ = 0.057 m?/s at 120°C. 


V. The stationary state. In the previous sections we have assumed a certain spectral 
distribution of turbulence at an initial time, ¢ = 0, and have studied the decay within a 
region, insulated from external sources of energy. However, in most actual cases the 
turbulence will be stationary and the mean intensity in every part of the spectrum 
will be constant with time. There is an input of energy in some parts of the spectrum 
and an equilibrium state is obtained, when the same amount of energy is converted 
into heat motion in other parts of the spectrum. 

|. The stationary equations. We shall now investigate how far conclusions on the 
homogeneous, stationary state may be drawn from a linear theory that starts with the 
time-dependent Eqs. (31), (32), (33) and (34) for an electrically conducting liquid in a 


magnetic field. Introduce the tensors 


ne eit e (132) 


OM e] (K) 
5 { a sone j 6. Dp ; = 
I r, t) FY T,,(r, t E;; (f, FY 


and 


Ev" r,t) = -—W [L.(r, 0) — L,(-—r, Od] (133) 


- € 
k 
Or, 
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(MH) 


. ny : (M K 4 
and the corresponding Fourier transforms ¢;;"’(«, t), «:;’ (x, 2) and €;;"’ (x, ¢), where 


6.(x, t) = (82°) iff E;;'@, 0 exp (—ix«-r) dr, drz drs etc. (134) 


It is clear that E‘*” (0, t) and E{* (0, 2) are the total rate of magnetic and kinetic energy 


(MH) 7 ° ° ° Pe 2 e 
losses and E‘{*” (0, ¢) is the decrease in magnetic energy per unit time and unit mass 
due to the magneto-hydrodynamic interaction as shown by Eqs. (29) and (30). Equations 


(31,) (32), (33) and (34) can now be written as 


Er, 2) = Ei", 2) — 2aV’M...(z, 0), (135) 
ESP, ) = —ES""(, 2) — 2V’K;,(r, t) (136) 
and 
0 . 2 |p MH); r2 a [ es Q7 
——(\+ VV’ IE;; (, ) = —2W. = [M,,0, ) — K,,(r, 0]; (137) 
ot Or; 
and the corresponding spectral equations are 
ar (x, 2) = &( "" (x, t) + Q2rx A; (x, 0), (138) 
eS (x, 2) = —ef M(x, t) + 2x’; ;(x, 2), (139) 
9 2 > v2 
| 2 + (X\ + net fee (x, t) = 2,W;[A,,(x, 0 — 2:;(x, b] (140) 
C 
and 
O(n, t) = —ix,W,IT,;(x, 2 — T;.(—x, dD], (141) 


where the tensors of Sec. III, 2 have also been introduced. 

If no external energy sources are present, the right hand sides of the contracted 
equations (138) and (139) represent the energies lost per unit time in the range dx, dx2 dx 
due to interaction and dissipation. In a stationary state this loss is balanced by an equal 
assumption that these con- 


input of energy from external sources. We now make the 
at r = 0 but can be general- 


siderations are valid not only for the contracted expressions 
ized to all components and all values of r, i.e., the stationary problem corresponds to the 
case when there is no dependence on ¢ in Eqs. (133), (135), (136), (137) and their spectral 
equivalents. The quantities E{” (r) and E‘* (r) now represent the external magnetic 
and kinetic ‘‘inputs of energy”’ per unit time. 

There is no formal difficulty in extending the foregoing results to the general, non- 
linear case. The various correlations can be formed in the same straightforward manner 
as in Sec. III and the derivations need not be given here. The stationary equivalents of 


Eqs. (138) and (139) become 


(Rn) = (ue) + 2d" A. (x) + a; |(%) (142) 
and 

Ee) = ee (un) + Qu? O, (x) + B;,(x), (143) 
where a;; and 8;; are the non-linear terms. The sum of these equations is 


e,,(w%) = 2dx°A,;(x) + Qux° 2; (%) + a;;(%) + B:;(x). (144) 


! 
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2. The terms of interaction. The equations in the preceding section form a starting 
point for a discussion of the physical consequences of different approximations. The 
simplest situation occurs for small amplitudes and without external fields when there 
are no interaction terms in Eqs. (142) and (143). If the external sources are taken away 
the turbulent fields of the magnetic and the kinetic energies decay independently of 
each other as is evident from Eqs. (52) and (53). In the stationary state every spectral 
range requires a separate magnetic and kinetic energy source; further there will be no 
exchange also of the energy in the different wave numbers and in different directions’. 
However, if an external magnetic field is introduced the expressions (142) and (143) 
show that there is now a coupling between the magnetic and the kinetic fields in every 
spectral range, also in first order. With an input ¢{%” (xo) = «,"” of magnetic energy per 
unit time and unit range of wave number space and an input «5“’ of kinetic energy at 


the wave number x, the equations (142), (143) and (140) become: 


6 = 60 (wo) + KoA s(Ho), (145) 
6 = E(u) + Dvn Ds ( 0) (146) 

and 
(X + v)wses (x9) = WnioWE[Asi(%0) — 2:i(o)]. (147) 


The solutions of these equations are 
A. i(%0) = { [0 + Ks + Foden"? + Foes} /{2Ko(a + v)(rvKs + Fo)}, Fo = KroW, (148) 
and 
Q;.(%) = {[AMA + »)Ka + Foden”? + Foeo?}/{2Ko(A + »)(AvKo + Fo)}. (149) 
In particular for W © we have 

As:/Q, = 1. (150) 


This is to be expected, since the losses become less important for large values of the 
external field and the magneto-hydrodynamic waves approach the “‘ideal’’ state, where 


V = h(u/p)'” =v (151) 


and there is equipartition between the magnetic and the kinetic energies. Especially 
for zero input of magnetic energy, the kinetic energy given by Eq. (149) decreases with 
a ratio 


2, ;(% ’ F, = co ) ‘2; (Xo ’ Fy = 0) => v/(r + v) (152) 


when the quantity /’, increases from zero to infinity. This conversion of three-dimensional 
turbulence into two-dimensional turbulence is consistent with the decrease in the time 
of decay given by Eq. (69). 

Finally, we may conclude as in earlier discussions of turbulence that an energy 
transport through the whole spectrum is only possible in the presence of the non-linear 
terms. The sum of these terms in Eq. (144) generally differs from zero and there is an 
energy flow through the spectral range in question. 


A discussion of the properties of the interaction terms in hydrodynamics may be found in G. K. 
Batchelor, The theory of homogeneous turbulence, Cambridge Univ. Press, Chap. V, 1953. 
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VI. Applications to experiments. Experiments on turbulent flow of mercury in a 
homogeneous magnetic field have been carried out by Hartmann *** and Lehnert *** with a 
maximum field strength of about 1 Vs/m’ (10* gauss). Hartmann has shown that the 
necessary pressure difference to keep constant volume flow between the ends of a rec- 
tangular channel decreases in the turbulent state when the magnetic field increases 
(similar investigations have been carried out by Murgatroyd’). In the latter experiment 
the torque between an inner, rotating and an outer, fixed cylinder decreases for a given 
angular velocity if a homogeneous, axial magnetic field is introduced. Thus, in both 
experiments the field seems to have the effect of suppressing the intensity of turbulence. 

These results may be given a qualitative explanation in terms of the linear theories 
of Secs. III and V. The energy sinks formed by the dominating dissipation for large wave 
numbers will cause a mean transfer of energy from the largest whirls, which receive the 
greatest part of the external input of energy. It is reasonable to assume that a great part 
of this energy will reach the energy consuming region for large wave numbers before 
being dissipated. In this region we are allowed to apply the linear theory and obtain a 
connexion between input power and turbulent intensity of the small scale motion. 
Since it is unlikely that the external magnetic field will have a reverse effect on the 
large scale motion we may expect that the linear theory would not depart too far from 
the real physical situation. 

From Eq. (71) and Fig. 2, giving the result of the theory of decay of turbulence, the 
damping effect on vortices having finite wave numbers in the direction of the magnetic 
field the reduction of the total isotropic intensity by a factor of about 2/3 is apparent. 
This implies that a given intensity of the kinetic energy Q;; requires a larger input of 
energy in the presence of the field than if the field were absent. Further, the tangential 
force due to turbulent exchange of momentum to the walls of the experimental arrange- 
ments increases with the intensity of the kinetic energy. Thus, a given pressure difference 
or a given torque requires an increasing input of energy for an increasing field strength, 
i.e. an increasing volume flow and an increasing angular velocity respectively. Con- 
sequently, if a magnetic field is present, the curves for the pressure difference and the 
torque as functions of the volume flow and the angular velocity respectively must lie 
below the corresponding curves in the absence of the magnetic field. 

These same conclusions may be drawn from Eqs. (149) and (152), representing 
stationary turbulence; these equations show how the intensity of the kinetic energy 
decreases with increasing field and constant input of energy. Even if we assume an 
input of both kinds of energy the ratio 


. ‘ K V . 
Oe ee ©)/0; ("x , F 0) = ve, (A+ ve, +. ) (153) 


will be very small, since \ > v for mercury and /*’ and E‘"” will probably be of the 
same order during the transport of energy towards high wave numbers. Consequently 
the conclusions both from the theory of decay and the stationary theory of turbulence 
are in agreement with experimental observations. 

Finally, it also follows from the results of Sec. IV, 4 that the Coriolis force does not 
play any essential role in the experiment with rotating cylinders. 

VII. Concluding remarks. A linear theory gives some information about the physical 
nature of turbulence in magneto-hydrodynamics, especially in the case of linear inter- 


'6Phil. Mag. 44, 1348 (1953). 
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action forces caused by external magnetic fields and Coriolis fields. However, this does 
not show the feature of the general turbulent interaction and the process of distribution 
of energy. One way to solve this difficult problem is to continue the work on the “quasi- 
statistical” theory of correlations. Another way, which at the present stage also seems 
to be associated with great difficulties, is to make a purely statistical attack, somewhat 
along the lines of Burgers’ and Onsager’*. But such a theory, even in the absence of 
external fields, is probably not available. 
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A NEW SINGULARITY OF TRANSONIC PLANE FLOWS* 


BY 
A. R. MANWELL 


University College, Swansea 


Introduction. The limit line singularity of the hodograph method [1,2,3] cannot 
arise in transonic plane flows past a physical boundary of finite curvature [4, 5, 6, 7]. In 
this note we discuss a new singularity of the hodograph method which gives infinite 
acceleration of the fluid but avoids the difficulties of the limit line theory. 

1. On the non-occurrence of the limit line in transonic plane flow. We follow 
Morawetz and Kolodner [7] in regard to notation and write the hodograph equations as 


+ hee p Qe ’ (1.1) 
vo=p ge [¢ —c ly. (1.2) 


Here ¢(q, @) is the potential function, ¥(q, 6) the stream function, p is a known function 
of g and c’ = dp/dp where p = p(p) is the equation of state. 
We define characteristic variables a and 8 by 


2da = pq 'k dq — dd = de — dé, (1.3) 
2dB = pq '‘k dq + dé = de + dé, (1.4) 
say, where 
= peg — 0), (1.5) 
so that k = k (a + 8), and k increases with q. 
Then (1.1) and (1.2) may be written 
Ga = —ky. , (1.6) 
gs = ky. (1.7) 
The limit line is defined by the locus J = 0, where 
A(x, y) eee ci ee 37,2 
= Fy Oe We — Reg Vi) = 7K Vas - (1.8) 
By (1.6) and (1.7), ¥ satisfies 
2hpasp + kala + Hs) = 0. (1.9) 


We shall prove that no limit lines occur if: 

(a) the physical boundary has finite curvature, say K = K(@), and 

(b) the hodograph image of this boundary has a continuously turning tangent 
which is never in the characteristic direction, say 8 = B(a) with 6’(a) finite and non- 
zero for the equation of the hodograph boundary. 

We first show that ¥, and y, , and therefore J, are bounded in the supersonic region 


*Received October 27, 1953. 
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and on the sonic line. As in [7] we have 
K(0) = —(pqJ)'y, . (1.10) 


Combining Eqs. (1.3), (1.4), (1.8) and (1.10), with ¥./~s = — B’(a) we see that, given 
K(6) and 8(a), both y, and Ws, take determinate finite values at each point of the hodo- 
graph boundary, i.e. knowledge of the physical boundary and of its hodograph image 
leads to Cauchy data for the hodograph solution ¥(a, 8) which is thereby uniquely 
determined in the supersonic region [6]. Clearly y, , and yz are finite at all supersonic 
points since any singularity would be propagated along at least one characteristic to 
the boundary where we know that y.. , ¥, are finite. For sonic points we take the origin 
of a, 8 on the sonic line and note that k ~ (a + 8)'”*. (We use the sign ~ to denote the 
form of the dominant term in an expansion near a singularity and the sign > when the 
numerical factor is significant). The dominant terms in the solutions of Eq. (1.10) must 


therefore satisfy 

G(a + B)Was + va t vs = O, (1.11) 
which is homogeneous in a, 8 while the terms neglected are of a higher order. Now y, , 
Ws are finite at all supersonic points near F, so the theory of Darboux [8, p. 90] applies 


as a, 8 tend to zero. 


In our notation 
Via, B) ~ (a + 8)” ‘g(a, 8) + hia, B), (1.12) 
where g and A are regular. 
We see at once that ¥, ~ ¥s — ¥. and y, ~ (a + B)'” (Wa + Ws) are bounded on 
the sonic line. Writing U for the values of ky, on the characteristic line a = const., 


Eq. (1.9) gives 


dU/dk = —y . (1.13) 

From (1.8) and (1.13) and for values on a = const. 
J = —@k*(d/dk)(U*k”). (1.14) 
If U vanishes for / ky and W, is bounded, say | y | < M, Eq. (1.13) gives | U | < 
M(k — k.) so that U*k™* certainly vanishes at k = k, . The case ky = 0, i.e. where the 

. . . . . . : 4° r27-1° we, 

zero is on the sonic line is included in this argument. Since U"k~ is positive near k = kp 
we deduce that it is increasing near k = k, and so the derivative in (1.14) is positive. 
We see then that / changes sign from positive to negative with k and q increasing and a 
similar argument would show that if V = ky, vanishes J changes sign along the character- 


istic 8 = const. 

Now conditions (a) and (b) also ensure that J > 0 on the boundary and the reasoning 
of Manwell [6] and Morawetz and Kolodner [7] can be applied to complete the proof 
that J is positive throughout the region. 

2. A singularity induced by a corner on the hodograph boundary. We now relax 
condition (b) of Sec. 1 by admitting a single corner at point C on the hodograph boundary, 
(see Fig. 1). The solution (a, 8) is still uniquely defined by giving 8 = B(a) and K = K(@) 
but the solution has discontinuities in the first derivatives with respect to a and 8. 

Observing that the leading terms in the expansions near C must satisfy Yas = 0 we 
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find without difficulty that the solutions in regions J and JJ may be written 


vi = dla(l + &)/(1 — ¢ + 8], (2.1) 
v2 = dla(1 + &/(1 + ¢’) + 8], (2.2) 
where the tangents at C are o = «0 and « = — ¢’6 while d is a constant and the origin 
is at point C. 
B a 







A Hodogroph boundary 


Sonic line 









B G F A 


Fic. 1. Hodgraph Plane, Supersonic Region in Characteristic Coordinates a, 8. 


The discontinuity in y, , say [Wa], is, see [9, p. 54] and Eq. (1.9) above, propagated 
along the characteristic a = 0 according to 


2k[ Pals + Kelva] = 9, (2.3) 
which gives, with C another constant, 
[vo] = Ck”. 2.4) 
Since y, is continuous at a = 0 
ive] = —3Ck'”. (2.5) 


Now the arguments of Sec. 1 still apply to the solution y, which is determined by finite 
data on the smooth curve AC independently of how AC is continued beyond C. In 
particular (Y¥,),, (¥,), and J, are finite along CF. Using Eqs. (2.4), (2.5) and (1.8) we 
find without difficulty that, at F, (¥,), = (W2), and J, = J, , which is non-zero according 
to Sec. 1, but (¥.), tends to infinity near F with k~'”’. According to Eq. (1.10) the stream- 
line curvature remains finite at F but for the fluid acceleration, say f, we have f, = q 
dq/ds = q(dq/d@),y conse. (€9/ds) = — (qKyWe/,)2 which tends to infinity with yy . 

In the limit line singularity both K and f become infinite together, a difficulty which 
is avoided here. 

3. Expansion of the solution near the singularity. We write [y] for ¥2 — W,.2) 
and set [fy] = 0 in region J while y,,.2) may be any flow past a smooth profile and de- 
termined in region JJ by data on AC and AC slightly produced beyond C with no dis- 
continuity in slope at C. Taking the origin at F we construct the leading term in the 
expansion of a solution of Eq. (1.9) as a homogeneous function of a and @ satisfying 
Eq. (1.11), viz., 


[y] _— ap~'’°F(1/6, 7/6, 2, —ap"), (3.1) 
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for 0 < — af" < 1 and with F(a, b, c, z) denoting the hypergeometric series (see [10, 
Chap. 10)). 
For a = 0 we have [y] = 0 and, in conformity with Eq. (2.4), 
[y.a]l~ Bs ~k™., (3.2) 

Now [y] is the only homogeneous solution satisfying this requirement and regular in 
0 < — a8” < 1. For| a6 | < 1 Eq. (3.1) may be differentiated term by term and we 
find that [y], ~ 6’”° and [y], ~ @~””°. For almost all points in the hodograph plane near 
F, z = — af” is nearly unity and we use the identity given by Copson [10, p. 252] 
yielding 
[vy] ~ a6” [PF(1/6, 7/6, 1/3, 1 + a6’) 

+ Q(1 + af™')*“*F(11/6, 5/6, 5/3. 1 + aB™')], (3.3) 


where P and Q are certain numerical constants. 
Writing g = c, + At we choose the constant d so that o = ¢ 


3/2 


and for all finite values 


of t/0 
[vy] = o(P + 2°°Qt0-*”*). (3.4) 
We can show too that Eq. (3.4) may be differentiated giving 
[v], = (5/6)Pe’”*, (3.5) 
[y]. = 2°°Qe”. (3.6) 


From these results for the derivatives of [y] it follows that J, — J,,(2) is uniformly 
small and that in spite of the infinity in acceleration at / the physical coordinates near 
F are changed very little by adding [y] to ¥,.2, . Although there is a discontinuity in 
y. along FC the physical solutions join up smoothly because, according to (1.6) and 
(1.7), not only y but ¢ is continuous between J and JJ. 

4. The continuation of the solution across the sonic line. In general, of course, 
we may not joint up piecewise differentiable hodograph solutions, for the corresponding 
physical planes would not join up. In particular we require continuity of both y and y, 
at the sonic line (see Fig. 2). We write 


13/2 


6 = 8 cos 9, |¢ |" = ssin 9, (4.1) 


with 0 < 9 < z=. Suitable solutions are found as the real and imaginary parts of 


Vv ~ 83°” exp [7i(x — 0)/6)]F(1/6, 7/6, 2, exp [—270]), (4.2) 
and two further solutions are found by replacing F by 
G = F | ZF" — Z)""* dZ (4.3) 


where Z = exp (— 270). 

Except in the vicinity of 6 = 0 we have either 6 nearly zero or © nearly and so Z 

is nearly unity in either case. Using Eq. (3.3) for complex values of Z and Eqs. (4.2), 

(4.3) it is not difficult to show that, except for @ nearly zero, the most general complex 
W is given by 

Vv ~ 6°” exp [7i(r — 0)/6)[R + S | sin © |*”* exp [i(r — 20)/3], (4.4) 
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where # and S are arbitrary constants. 
Taking S zero we find a solution 


v, = 0° sin [7(r — 0)/6], (4.5) 
and similarly making R zero we get 
Vv, = 6° | ¢| sin [11(r — 0)/6). (4.6) 
Once more it can be shown that these expressions may be differentiated with respect 
to t or 6 and easily d6/dt = 0 for t = 0 so both derivatives vanish on 8 = 7, i.e. Vy = 
Vv, = 0 for 6 < 0. We form the combination 
[w] = 2Pv, + 2°°Qv, (4.7) 
and verify that it satisfies both ¥ = ¥, = 0 for |¢| = 0 and 6 < 0 and conditions (3.5) 
and (3.6) for |¢| = O and @ > 0. 


Line of discontinuities 





in the derivatives e 
70*0 +0 
* 
=O 
a I 
q I 
Supersonic region r i 
¥,- ¥,7 Lv] A Sonic fine (itt =O) 
@ increasing J 
F 
Subsonic region 
@=0 e-:n 





Fic. 2. Hodograph Plane Near Critical Point F. 


A simple discussion shows that the other two linearly independent solutions given 
by Eq. (4.4) cannot occur and so [¥] which was uniquely determined in region II, save 
for a constant multiplier depending on data near C, has a unique continuation across 
the sonic line, viz. [¥] of Eq. (4.7). It is easily verified that for any ¥ and near 6 = 0 we 
have V, ~ | t|"'“* and ¥, ~ | t |"'”*. (We use Eqs. (4.2) and (4.3) and prove that F(Z) 
has no roots on | Z | = 1). Noting also Eqs. (4.5) and (4.6) we see that W has finite 
derivatives at all subsonic points near F. The solution in the physical plane can be con- 
structed by adding [¥] to the continuation across | | = 0 of the function y, and perform- 
ing the usual integrations. 

The discussion of Secs. 2, 3, 4 applies with obvious changes to the propagation of 
discontinuities along the other characteristic CG and a singularity develops similarly 
at G (Fig. 1). 
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5. The occurrence of the singularity in flow past a given profile. According to 
Sec. 1, transonic flows can be constructed in the neighborhood of an arbitrary convex 
physical boundary of finite curvature along which the velocity is determined by a smooth 
hodograph image. By Sec. 2, however, there are infinitely many adjacent velocity dis- 
tributions for which the flow must break down at the sonic line, viz. those for which the 
hodograph images have any discontinuity in slope at supersonic points. 

We turn now to the possibility that the singularity described above may arise for a 
given body in a subsonic stream as the Mach number of the stream is slowly increased. 
It must be admitted that the existence and uniqueness of plane transonic flow solutions 
in this situation is by no means certain. As for uniqueness, the author has proved [11] 
that the solutions for cyclic flow round a circular cylinder are not unique at an enumer- 
able infinity of values of the circulation parameter and that these values lie close together. 
For non-circulatory flows past a precisely defined body it is possible that existence of the 
solutions will fail at certain (perhaps ‘almost all’) Mach numbers with transonic flow 
conditions (c.f. Frankl [12], although the argument of [12] is by no means complete) 

On the other hand in a number of approximations worked out by the present writer 
it seems clear that the general consequence of increasing the strength of the stream for 
an approximately determined body is to give velocity distributions whose hodograph 
boundaries develop a corner at the point of maximum velocity in the supersonic region. 

An example of this process is furnished by [13] which although, as appears on closer 
examination, merely an approximate solution of the problem implied by the title is a 
useful contribution to the problem of the behaviour of transonic flow when the boundary 
is kept fixed while the stream conditions are varied. (The solution of [13] satisfies correct 
boundary conditions on the arc BD and on the axis AB, (Fig. 3) and the flow becomes 


E 








Fic. 3. 


uniform at great distances but, owing to lack of analytic continuation the velocity vector 


along DE need not be parallel to AB except at D and £.) 
According to [13] the image of the boundary is determined in the hodograph plane 


by an equation of the form 
cos 6 = A(q)/A(q), (5.1) 
where g; is the maximum velocity on the boundary. If 9q, d2 , Where \’(q.) = 0, this 
equation may be expanded as 
g— a = —4A(q)) Na) @. (5.2) 
If however 9g; = @2 the expansion must be written as 


’"(q2)(q — G2)” = —A qr)’. (5.3) 
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For the semi-circular body which has in fact been treated in [13] it was found that 
\ & g/e — (5/24)(q/e)® — [125/3456 + 35(y — 1)/576](q/c)’, (5.4) 


where the notation is that of Courant and Friedrichs [9]. To this approximation it appears 
that \’(q) vanishes for g. = 1.03 c and that the branching streamlines in the hodograph 
plane make an angle are tan (qd@/dq) = 1.86 with the line 6 = 0. 

In this example the singularity develops for only a very small supersonic region but 
in other cases of transonic flow worked out by the author to the same degree of approxi- 
mation the supersonic region appears to develop considerably before the flow breaks 
down. 

REFERENCES 
1. G. 1. Taylor, Recent work on the flow of compressible fluids, } Lond Math. Soc. 5, 224 (1930). 
2. Th. von Karman, Compressibility effects in aerodynamics, J. Aero. Sci. 8, 337-356 (1941). 
3. H. S. Tsien, The viting line in mixed subsonic and supersonic flow of compressible fluids, T.N. 961, 
N.A.C.A. (1944). 


\. A. Nikolskii and G. I. Taganov, Gas motion in a local supersonic region and conditions of potential 


j 
ransl. from Prikl. Mat. Mekh. 10, 481 (1946) by S. Reiss, N.A.C.A. Tech. Mem. 


flou breakdown 
1213 

5. K. O. Friedrichs 
Appl. Math. 1, 287-301 (1948) 

6. A. R. Manwell, A note on the hodograph transformation, Quart. Appl. Math. 10, 177-184 (1952). 

7. C. 8S. Morawetz and I. I. Kolodner, On the non-existence of limiting lines in transonic flows, Com. 


‘ome 
Pure Appl. Math. 6, 97-102 (1953). 
8. G. Darboux, Lecons sur la théorie générale des surfaces, vol. 2, Paris, 1915. 
R. Courant and K. O. Friedrichs, Supersonic flow and shock waves, Interscience Publishers, Inc., 


ind D. Flanders, On the non-occurrence of a limiting line in transonic flow, Com. 


New York, 1948. 
10. E. T. Copson, Theory of function of a complex variable, Oxford, 1935. 


11. A. R. Manwell, The variation of compressible flows, Quart. J. Mech. Appl. Math. 7, 40-50 (1954). 

12. F. I. Frankl, On the formation of shock waves in subsonic flows with local supersonic velocities, transl. 
N.A.C.A., T.M. 1251. 

13. A. R. Manwell, Expansion in series of the exact solution for compressible flow past a circular or an 
elliptic cylinder, Phil. Mag. (7) 36, 499-510 (1945). 


























351 


CANONICAL EQUATIONS FOR NON-LINEARIZED STEADY 
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1. Introduction. A steady flow is conical if it contains a vertex, P, such that the 
velocity components, pressure, density, and entropy are constant on every half line 
emanating from P. Numerous authors [e.g. 5, 9, 11, 12, 19] have investigated linearized 
conical flows, which can be attacked with the aid of complex function theory inside the 
Mach cone of the vertex. A few [2, 16, 22] have considered second order approximations 
or perturbed axisymmetric non-linearized conical flow [20]. So far, the study of non- 
linearized conical flow has been confined to special examples [3, 4, 10, 13, 21], or the use 
of relaxation methods [14] or three-dimensional characteristic treatments [15]. It is 
natural to ask how far the discussion of the non-linearized case can be carried. Although 
no attempt has been made in this note to reach the ideal conclusion, an existence and 
uniqueness theorem, the present inquiry does lead at least to the formulation of an 
interesting boundary value problem and to a basis for numerical computation. 

Most of the discussion is concerned with the flow about a conical body completely 
surrounded by a conical shock. The change in entropy at a shock wave is of higher 
order than the changes in the other flow functions. For this reason the equations of 
supersonic flow are frequently simplified by neglecting the variation of entropy behind 
the shock, as will be done here. Now the flow field may be but is not necessarily irro- 
tational. In this note the flow will be assumed to be irrotational. This makes it possible 
to reduce the partial differential equations to simple canonical forms, to determine the 
class of transformations under which these are invariant, and eventually to map the 
conical surfaces of obstacle and shock onto known curves. Thus the difficulties which 
arise from not knowing the location of the shock wave can be avoided. Finally, it should 
be remarked that Taylor-Maccoll flow has all of the assumed properties for the com- 
pletely elliptic case. Accordingly, this special example is a fruitful source of conjectures 
as to the influence of the various parameters and of estimates for numerical computation. 

2. Fundamental equations. The velocity potential function, ¢, of a steady irro- 
tational flow satisfies 

(a°5,, — u,u,;)d°@/dx,; dx; = 0, i,j = 1, 2, 3, (2.1) 
where the velocity 
u; = dp/dz; , (2.2) 
the squared velocity of sound 
a’ = b(y sas 1)(1 o UUs), (2.3) 


*Received Dec. 11, 1953. 
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x; are cartesian coordinates, Kronecker’s 6;; = 1(0) accordingly as i = (#)j, and 
repetition of indices implies summation over their ranges. In a conical flow with vertex 
(0, 0, 0) let 

@ = 20(X), Ae = &,/2, Z2=97;3, a = 1, 2. (2.4) 


Note that points in the X,X, plane correspond to lines through the origin of the 
Z,X_%,-space, and conversely. Now (2.1) becomes 


[a°(dag + X.Xs) — (Ua — WXq)(ug — wN5)] 0°@/OX, AX, = O, (2.5) 

where 
Ug = O6/dX, , (2.6) 
w= = ©@—u,Xa. (2.7) 


If u,.(X) are functionally independent, subject (2.5) to the Legendre transformation 
(2.7) to obtain [4] 
[a°(bag + WaWs) — (Ua + WWa)(Ug + wws,)](—1)*"* d°w/dua., OUg,, = 0, (2.8) 
where a + 1 and 6 + 1 are to be reduced mod 2, and where the velocity u, , u. , w(u, , ue») 


is assigned to the point 


X, = —dw/du, = —W, (2.9) 
Now introduce more general parameters pa = Malu; , M2). Then 
Uz; = Ups » Me)- (2.10) 
Let 
q = uu; q 0q/Op, = U; OU;/OMe ; Jag = (OU;/Op,)(OU;/Opusg), 
(2.11) 
bas = (911922 — Po * det || 07x, ‘Opa Og , Ou;/Op, , Ou; /Ope 
For the special choice nu, = wu, (2.8) implies 
(a°gas — 7 9q/Apa 9g/Aus)(—1)**basips: = 0. (2.12) 


Since gas , bag , and dqg/du,. are covariant tensors, (2.12) is a tensor equation. Hence 
its form is independent of the particular choice of parameters u, . Now observe that 
(2.9) implies x,du,;/du. = 0. Thus the point u,;(u) on the hodograph surface corresponds 
to the half line 

XL; =1V;, v; 0u;/Ou, = 0, (2.13) 


where v;(u) is a vector normal to this surface at u;(u), and r is a parameter. 

So far the three functions u;(u) have been subjected to (2.12) only. Two more equa- 
tions are required to determine u; . As suggested by [6, p. 130], they can be chosen to 
simplify the coefficients of (2.12) in a way which depends on the nature of the quadratic 
form 

(a° Gas _ q 0q/Ope 9q/Ops) Lag . (2.14) 


This is definite (indefinite) if and only if the determinant of its coefficients is greater 
(less) than zero. For u, = u, and by (2.9) the determinant becomes a*[(1 + X,X.)(a° — 
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q) + (u.Xq + w)’]. In subsonic flow (2.14) is always definite. In the supersonic case 
(2.14) will be definite (indefinite) if and only if the half line (2.13) lies inside (outside) 
the Mach cone based on u;(u) and with vertex at (0, 0, 0). In the former (latter) case 
equations (2.5), (2.8) and (2.12) will be said to be of elliptic (hyperbolic) type for reasons 
which will become apparent in Sec. 3 and 4. 

3. The hyperbolic case. This section has been included to round out the discussion, 
for the sake of contrast with the material in Sec. 4. It is not essential for the sequel. 
Examples of conical flows which contain regions where x; = rv; is outside its Mach 
cone have been discussed in [10]. 

If (2.14) is indefinite, choose the parameters y, and y, so that 


A =a@u;u;, — (uu)? = 0, (3.1) 
B= Wujttio — (Uj) = 9, (3.2) 

where u;, = 0u,;/du,. . By (2.11) and (2.12) 
Uns = Corti (3.3) 


for some functions C, . Eliminate u;,. from A, = B, = 0 to obtain 
(0°gas — 9 Ga%)C3 = WaGiz2 — WasiJaa (a not summed). (3.4) 


Conversely, suppose u,;(u) is a solution of (3.3), where Cs; are defined by (3.4). Then 


A B, = 0, so 
A=A(u), B= Bw). (3.5) 


Equation (3.3) is a hyperbolic system of second order quasi-linear partial differential 
equations with the same principal parts, and yu, and yu, are characteristic variables. The 
existence and uniqueness, in the small, of solutions of non-characteristic and charac- 
teristic initial value problems for such systems have been discussed in [6, pp. 316-326] 
and there are numerous references in [1]. In conical flow the initial conditions of fluid 
mechanical origin, discussed in Sec. 4, appear to be more complicated than those ordi- 
narily considered in the literature. To these must be adjoined an initial condition such as 


A=B. (3.6) 


By (3.5) this implies A = B = constant. Then it will suffice to impose A = 0 at one 
point, to obtain A = B = 0. 

To simplify the formulation of initial value problems, it would be desirable to stand- 
ardize initial curves in the u,u.-plane as far as possible. This raises the question as to 
what class of transformations 4. = ua(u*) will preserve the forms of (3.1) and (3.2). 
The coefficient tensor for the transform of (2.12) is 


a*g*,; — q* tq’ = (a°ga5 — Y 09s) (Ope ‘Ou*) (Ou; /Ou4s). (3.7) 


By (3.1) and (3.2) its diagonal elements vanish if and only if (0u;/du%4)0m2/du4 = 
(Ou, Ou%)Ou./du% = 0. Hence either uw, = u,(u%), wo = we(ud), or else uw, = ui(u4), 
U u.(u*). That is, the most general transformation with the desired invariance consists 
of changes of scales along the axes, possibly combined with reflection with respect to uw, = wh, . 
By such transformations two initial curves can be mapped simultaneously onto straight 
lines. 

Examples such as those considered in [10, 15] make it seem likely that every conical 
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flow which contains regions where the equations are of hyperbolic type also contain 
regions of elliptic type. Hence the canonical form (3.3) cannot be used for an entire flow. 
4. The elliptic case. First it should be remarked that at all points in a Taylor- 
Maccoll flow the ray rv; is inside its Mach cones. Hence (2.14) is definite throughout. 
It seems reasonable to conjecture that this will also be the case at least for slightly 
yawing or slightly distorted circular cones. 
If (2.14) is definite, choose isothermic (relative to the tensor a’°gag — 9°qa9s) param- 


eters uw, and pp , i.e. set 


A = WU lin — UU, ,UjUj2 = O, (4.1) 
QB = Wut — Usgtin) — (us)? + (uj)? = O. (4.2) 
Then by (2.11) and (2.12) 
Usaa = Cyttip (4.3) 
for some functions C,; . Now form 
0A/du, — OB/dp. = O, 0A /dOp. + OB/dpu, = O (4.4) 
and eliminate u;,, to obtain 
(a°gas — F9a9)Cs = 3Gas(a” + Qa — Gass - (4.5) 
Since 


Conversely, suppose u,(u) is a solution of (4.3), where C, are defined by (4.5). 
this implies (4.4), A + 7B is an analytic function of uw, + tue. 

Equation (4.3) is a system of elliptic second order quasi-linear partial differential 
equations with the same principal parts. Since the coefficients Cs, are rational functions 


of their nine arguments wu; , u;,; , and u;. , the solutions of (4.3) are analytic [6, p. 339]. 


The existence and uniqueness of the solutions of Dirichlet’s problem for such a system 
for a small enough region of the yu,y.-plane have been discussed in [6, p. 287]. Unfor- 
tunately, the conical flow problem leads to the non-linear boundary conditions described 
below, rather than to the linear Dirichlet, Neumann, or mixed type usually considered. 

As a preliminary to the formulation of a boundary value problem, first determine 
= u(u*) that will preserve the forms of (4.1) and (4.2), 


the class of transformations pu, 
(3.4 becomes 


and therefore of (4.3) and (4.5). The transformed coefficient tensor 
a**g*, (9q:) |(Ou./Ou*) (Ou, /du%). Then p* will also be isothermic 
parameters if and only if (Ou,/du*)(Ou./du%) = Dé. for some function D. Hence, 
any conformal map of the u,u2-plane onto the w*y*t-plane has the desired invariance. 
Suppose that the conical body (shock) is described by f(X, , X.) = O(g(X, , X.) = 0). 
Assume that f = 0 and g = 0 are closed curves in the X,X,-plane which do not intersect 
0 is inside g = 0. Suppose that their images 


#2 kK - Ia* 
g8— QT dads 12 Ji 


themselves or each other, and that f 
F(u) = Oand G(u 0 in the y,u.-plane have the same properties. All of these conditions 
can be satisfied by Taylor-Maccoll flow, and should also be attainable in slightly per- 
turbed Taylor-Maccoll flows, at least. Then the region between F = 0 and G 0 can 
be mapped conformally onto the annulus 1 S | ut + wu%| < R for some R{7]. Assume 
that this has been done, that the body maps onto the unit circle, the shock onto the 
other boundary. Hereafter replace u* by u. . Now proceed to the boundary conditions 
for conical flow. First observe that to assure that (4.3) and (4.5) will imply the original 
equations, the analytic function A + 7B must be forced to vanish. This can be done by 
imposing 
A =@ (or B = 0) (4.6) 
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on both boundaries, and 


B=0 (or A = 0) (4.7) 
at one point. Also on | wu, + tq. | = 1 the equation of the cone 
f(X1 , Xe) = 0 (4.8) 
and the condition of tangency 
(ue. — wX,) 0f/aX, = 0 (4.9) 
must be satisfied, where XY, = 2,/z must be expressed in terms of du,;/dus by means 


of (2.13). Let (0, 0, Q) be the undisturbed velocity ahead of the shock. At the (actually 
unknown) conical shock g(X) = 0, the change in velocity must be normal to the shock, 


1.e. 
(0g/0X.)/Uq = (Xs 0g/dXs)(Q — w) (a not summed). (4.10) 
These equations imply ég/AX, = Ku, for some function K, and also 
Q= X,u, + w. (4.11) 
Then by (2.7), ® = Q, ie. the shock wave is an equipotential surface. Conversely, any 


equipotential surface = Q = constant satisfies (4.10). Finally, on g = 0 the equation 
of the shock polar 


Utely — 1+ 2Q° — (y + DQu] = (Q — w)*[y7 + DNQu — 7 + 1) (4.12) 


must be satisfied, where XY, must be determined as above. 


If the constant Q is prescribed, then R has to be found by first determining the flow. 

It is more convenient to prescribe R and to impose 
R dQ/ds = (—p. 0Q/du,; + uw, 0Q/du.) = 0, (4.13) 
on | uw; + tu, | = R, where Q is defined by (4.11), and where the constant value of Q 


is to be found after the flow has been determined. The shock wave must likewise be found 
by (2.13) on this circle after a solution has been found. This is comparable to the situation 
in numerical computation of Taylor-Maccoll flow, where for prescribed cone angle and 
surface velocity the free stream velocity and shock angle are not known until the end 
of the computation. 

To recapitulate: The problem is to find on the annulus 1 S | uw, + qu. | S R for 
given R and f(X) a solution of the elliptic system (4.3), (4.5) which satisfies (4.6), 
(4.8), and (4.9) on | wu, + tu, | = 1; (4.6), (4.12), and (4.13) on | uw, + tu. | = R; and 
B = 0 (or A = 0) at one point. Note that although the free stream’s speed is unknown, 
its direction is fixed. Hence (4.8) determines both the shape of the body and its orienta- 
tion relative to the undisturbed flow. It should be stressed that mere solution of this 
boundary value problem is not enough to guarantee the existence of irrotational conical 
flow about a given cone. In the units employed here, 0 S q S 1; also Q must be super- 
sonic. Thus a solution is not acceptable unless (y — 1)'?(y + 1)7'? < Q < 1. Fur- 
thermore, the map (2.13) from the hodograph to the physical space must yield a single 


valued velocity field. 
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5. Application to Taylor-Maccoll flow. To find the axisymmetric flow about the cone 


f = (X,X.)"” — tan & = 0, (5.1) 
of semi-vertex angle 6, , first seek a solution of the form 
u, = U(t) cos ¢, u, = U(#) sin ¢, w= Wid, (5.2) 
t = (X,X.)'”’, tand = X,/X, . (5.3) 
For the non-isothermic parameters yu; = ¢, uw. = ¢, the definitions (2.11) and (2.13) yield 
qg=U+W’, gi = U" + WwW", fio = O, Jo. = U’, 
b, =(UW" — U"W'\(U” + es. bj. = 0, b, = UW(U" + WwW”), 
Z, = —rW’ cos ¢, x, = —rW’sin 9, 22 FU’. (5.4) 
By (5.3) and (5.4) 
W’ + tU’ = 0. (5.5) 
By (2.12) 
W’ =a U/fa°(1 + #) — (U — tW)’), (5.6) 


where WV” has been eliminated by means of (5.5). At the shock (4.11) and (4.13) yield 
Q(t) = tU + W, 5.7) 
H(t) = U?(2Q? — (1 + Oy + DOW — 7 + 1)} = 0. (5.8) 


At the cone ¢ = tan @ q 
Ue = Welo (5.9) 


for U(t,) = U,, Wit = W,. 

A short digression on the boundary value problem (5.5) to (5.9) will be instructive. 
The problem can be replaced by an initial value problem by choosing W, > 0 such 
that W2(1 + #) < 1, defining U, = Wot, , and integrating (5.5) and (5.6). It can be 
shown that —U(é), U? + W’, Q(t), and —H(t)/U*(t) are increasing functions of ¢, 
that H(t,) > 0, and that H = 0 before a°U or a° U/W’ can vanish. By (5.7) and (5.8) 


0< Ut = 201 —(L+ €*)A7V/7 + DU + FQ 


where A? = 1/2(y — 1)(1 — Q’). Hence Q is supersonic. If Wo(1 + t) 2 1, QO > 1 
at H = 0. Thus some solutions of (5.5) and (5.6) are rejected by the restriction Q < 1. 
It is also known that there exists q*(¢,) such that for W,(1 + ey" = q* (to), (U™ 4 al la 


is subsonic all the way up to the shock; for g*(t.) < Wo(1 + @)'? ¢ (fy — 1779 4+1°" 
the flow is subsonic near the cone, supersonic near the shock; and for (y — 1) “(vy + 
I’? < Wil +t < 1 the flow is completely supersonic. To obtain completely 
supersonic flow and Q < 1 one must have 4; < 2/(y — 1). For y = 1.4 this yields the 
crude estimate 6 66°, by contrast with the computed bound of 57.5°. A similar 
crude bound for finite free stream Mach numbers can be obtained by replacing the upper 
bound 1 for Q by M[M* + 2/(y — 1)]”’”. In general, for a given (small enough) ¢ 


there is a minimum free stream Mach number for Taylor-Maccoll flow. For greater 
Mach numbers and the same f, there are two solutions, one of which has a strong, the 


other a weak shock. 


























‘ NON-LINEARIZED STEADY IRROTATIONAL CONICAL FLOW 357 


19 


— 


nt 
v 


a 


Isothermic parameters u* = T(t), uv = ¢ can be defined by 
[a(U"? + W”) — (UU’ + WW’)’\(dt/dTY’ = aU’. 


By (5.5) and (5.6) this becomes (d7'/dt)?> = W’/Ut’. Then choose 


at 


T() =14+ | (W/Ue)'? at. 


At the shock let ¢ = tan 6, = tl, , and T(t,) = R. 

R=1+] (w/uey'” dt. (5.10) 
In polar coordinates 7’, ¢, the cone corresponds to T’ = 1, the shock to T = R. Approxi- 
mate values for R as functions of free stream Mach number and cone semi-angle for 
y = 1.4, computed on the ENIAC and ORDVAC, are shown in Fig. 1. The curve for 
6 = 40° (not shown) practically coincides with that for 6 = 35° above M = 3. 
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Fig. 1. Outer radius for Taylor-Maccoll flow. 


6. Generalities and speculations. A. Numerical solutions of the boundary value 
problem can be attempted by relaxation methods. Estimates of the outer radius, R, 
of the annulus for prescribed free stream Mach number can be based on Fig. 1. The 
treatment of boundary conditions and the choice of coordinate lattice will be greatly 
simplified by the transformation pt + iu% = log (wu: + tu.) which maps the annulus 
1 < | yu, + w.!| S R onto the rectangle 0 S wt S log R, 0 S vt S 27. Computation 
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could be directed toward a great variety of goals, such as (i) comparisons with Kopal’s 
yawing cone calculations, based on A. H. Stone’s theory [20]; (ii) checks of Moeckel’s 
three-dimensional characteristic method for calculating the same types of flows [15]; 
(iii) use of computed results to estimate the importance of Ferri’s singularity in the 
entropy [8]; (iv) determinations of limitations on shape, such as (a) maximum inclination 
between a surface element and the incident flow; (b) maximum eccentricity for elliptic 
cones, i.e., can one pass to the limiting case of the delta wing at an angle of attack; 
(c) maximum yaw for a circular cone, i.e. can the angle of attack exceed the semicone 
angle? 

B. Since Taylor-Maccoll flow is a solution of this type of boundary value problem, 
it might be possible to apply the iterative method described in [6, p. 287], Green’s func- 
tion for an annulus being known, to show that flows can be constructed about slightly 
yawing almost circular cones 0 = f(X, , X.) = X,X. — tan’ 6 + eg(X, , X2) for 
reasonable g (not a function of X,X, only) and small enough e. 

C. M. Shiffman [17] has recently developed a proof of the existence and uniqueness 
of subsonic plane flows about fairly general shapes. This involves reformulating of the 
plane flow problem in variational form and then applying direct methods in the calculus 


of variations. It can be shown that this is impossible for the system (4.3), (4.5), i.e. 
that these are not the Euler equations for any integral 
J = I F(t , bo , U, OU/Opu, , OU/Opo) du, dus . (6.1) 


D. Can a new linearization of the conical flow equations be based on (4.3), e.g. by 
suppressing the right members? Instead of defining in advance the transformation 
from the X,X,-plane to the yu,y.-plane, as is done in the usual linearization, use (2.13) 
for this purpose. If possible, retain all of the boundary conditions in Sec. 4. 

7. Remarks on plane flow. Shiffman [17] and Bers [unpublished results] have 
recently established existence and uniqueness theorems for completely subsonic flows 
about plane profiles. Accordingly, it may be significant that their boundary value problem 
can be restated, as shown below, in a form intimately related to the problem of Sec. 4. 

Since u.(%; , 2%) are functionally independent in subsonic plane flow, apply the 


Legendre transformation 
ku, , Ue) = b — Vella ; (72) 
to (2.1) to obtain 
(a°5ag — Uqls)(—1)**? d°k/duas, OUgi, = O, (7.2) 


where a + 1 and 8 + 1 are to be reduced mod 2, and where the velocity u, is assigned 


to the point 
La = —Ok/du, . (7 3) 
Now introduce the more general parameters wa = Ma(U , U2). Then 


Ua = Ualtr » Me)- (7.4) 


Let 


Q 0Q/dusg = Uz OU,/Opsz , Gag = (du,/Ou—)(du, /Aus) (7.5) 
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and let the third order determinants 
Kap dk/Ops 
det = Bus ’ (7.6) 
| Uy, ap du, /Ous 
where k,.g and u,,.s are the second covariant derivatives of the scalars k and u, with 
respect to u. and us based on Gg . For the special choice nu. = u, (7.2) implies 
(a’Gag — Q 8Q/duq 9Q/dus)(—1)* Basi gsr = 0. (7.7) 


Since G., and 0Q/dug are tensors, while B,, is a relative tensor of weight one, the form 
of (7.7) is independent of the particular choice of parameters », . Also observe that 
(7.3) implies 


La WUa/Ipg + Ok/Ops = O. (7.8) 
For subsonic flow choose isothermic parameters defined by 
A = Waar — Uglaitglgs = 0, (7.9) 
2B = a’(Ualar — Ualar) — (Water) + (Water) = 0. (7.10) 
Now by (7.6) and (7.7) 
k/Ape One = Dz dk/dpsz , Ou, /Ope Ota = Dz du,/Oug . (7.11) 


Since (7.5), (7.9), and (7.10) can be obtained by setting w = 0 in (2.11), (4.1), and (4.2), 
then by analogy with (4.5) 


(a°Gas — Q 0Q/du. 8Q/dus)Dzg = 0.5G53 A(a” + Q”)/Oua — Gag Oa" /dys « (7.12) 


In the present case Ds, are rational functions of the six arguments u, , 0U./dus . As in 
Section 4, if k(u) and u,(u) are solutions of (7.11) where Dg are defined by (7.12), then 
A + 7B is an analytic function of uw, + tue . 

Suppose that the airfoil is described by f(x, , x2) = 0. Assume that f = 0 is trans- 
formed into a closed curve F(u; , w2) = 0, and that the flow is mapped onto the exterior 
of F = 0. Since (7.9) and (7.10) are invariant under conformal transformation, it may 
be assumed with no loss of generality that F = 0 is | wu, + qu. | = 1, and that the point 
at infinity in the x, + ix, plane maps onto the point at infinity in the u, + tu, plane. 
On the unit circle impose the boundary conditions A (or B) = 0; f(x, , x2) = 0, where 
xq are defined by (7.8); and u,df/dxz, = 0. At infinity u. = U,. = constants; and 
2, + ix, = ©. Finally, at some point impose B (or A) = 0. 
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FREE VIBRATION OF A RECTANGULAR PLATE WITH DAMPING 
CONSIDERED* 


BY 
MILOMIR M. STANISIC 


Armour Research Foundation of Illinois Institute of Technology 


Introduction. This paper presents a method for calculating the natural frequencies 
of the normal modes of free vibration of a rectangular plate, fixed along each edge, with 
arbitrary shape ratio. Viscous damping of the plate material is considered. 

Generally speaking, the objects of this paper are (a) to present a simple and practical 
solution for the stated problem, (b) to develop this solution in closed form so that a 
designer can evaluate the natural frequencies of a clamped plate of any shape ratio in a 
relatively short time, and (c) to make possible the determination of the damping co- 
efficient by comparing experimental and theoretical values. 

The theory developed is subjected to the following restrictions: (1) the plate is 
composed of a material which follows Hooke’s law; (2) the deflection of the plate is 
small compared to its thickness; (3) the thickness of the plate is small compared to its 
lateral dimensions. The method which is used to solve this problem is that of Galerkin’ 
[1] and belongs to the same general class as those of Rayleigh and Ritz. This method 
can be used (a) to determine an approximate solution of a differential equation with 
given boundary conditions admitting only the functions which satisfy the prescribed 
boundary condition exactly, and (b) to treat the problems which belong to non-con- 
servative systems. This method was described by E. P. Grossman [2] and also by W. J. 
Duncan [3], [4]. 

The degree of accuracy expected can be increased by increasing the number of 
independent functions which are used in the solution. 

W. J. Duncan has shown in his papers that when the functions are well chosen an 
excellent approximation can be obtained by use of a very small number of admissible 
functions. Since our system is non-conservative this method is most convenient for the 
solution of the problem. 

Solution of the problem. If damping forces are proportional to velocity, the motion 
of the plate is governed by the following partial differential equation: 


k 
w.(z, yy, = —— wz, y, 0, (1) 


V‘w(2, y, + rs D 


D 
where w(x, y, ¢) is the transverse deflection of the plate, k is the damping coefficient, h is 
thickness of the plate, p is the mass density of the plate material and D = Eh*/12 
(1 — »’) is the flexural rigidity of the plate. For a uniform plate vibrating harmonically 
with an amplitude ¢(z, y) we can write 


w(x, y, t) = (2, y) exp (—at) cos wl, (2) 
*Received Dec. 16, 1953. Presented at Eighth International Congress on Theoretical and Applied 
Mechanics, August 20th to 28th, 1952, Istanbul, Turkey. 
‘Numbers in square brackets refer to the Bibliography at the end of the paper. 
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where w is the angular frequency. Equations (1) and (2) lead to 


V‘o(z, y) cos wt + (2, ME (a? — ww’) —a | COS wt 


\ 


} Pa. 
oh (2 5 aw — w x) sin ut} = Q. (3) 
Since Eq. (3) must be satisfied for all values of ¢ we obtain (w/D) (2pha — k) = 0. 
Therefore 
k 
q= 2 ph’ (4) 


Now, Eq. (3) becomes 


; k ca ‘i 
VY 'o(x, y) | D 3 (2° — 0’ fore, y) = 0 
or 
V‘o(z, y) — Ad(z, y) = 0,7 (5) 
where 
= = Zao)" | , 
~ 4phD E +( k ; (6) 





& 4 
Fic. 1. Coordinate System for the Plate 


Assume the solution of Eq. (5) to be given in the form 


m n 
o(z,y =D daX,¥, (r,s = 1, 2,38, ---), (7) 
r=1 s@#1 
in which the X, and Y, are admissible functions that satisfy only geometrical boundary 
conditions, but need not satisfy any “natural boundary conditions’, and a,, are the 
amplitude coefficients. Since the function ¢(z, y) is an approximate solution of the 
problem, then the error caused in the solution is of the magnitude 


e= > doa, [V(X,Y.) — X,Y, ]. (8) 
r=] s=1 
The criterion that the approximation is best if « tends to zero requires, as stated by 
Galerkin, that the integrals 


| 
a b 


J = | | ee dx dy (p, q= i. 2, 3, a +) (9) 


“0 #0 
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be a minimum. The lateral dimensions of the plate are a and b. Equations (8) and (9) 
lead to 


pa e 


yreL | 


“0 #0) 


[7*(X,Y,) — AX,Y,]X,Y, dx dy = 0 (10) 
which can be written in the form 
> , & @:. | | [A c.asant o Y Secs l ocee Ack ccs ~~ Abel ike Oe =6. iD 


Equation (11) represents a system of linear homogeneous equations in the unknown 
coefficients a,, . The natural frequencies \, , A. , «++ are determined from the condition 
that the determinant of the system must vanish. 

Let the appropriate characteristic functions X, , Y, be given, following Young [5], 
in the form 


g g : g . g 
WV .(€) = | cosh p, . COS BL, i) ~ 8,\ sinh pz, — sin Me} 3 (12) 


such that 
if v=, ¥,=X,, £ = 2g, | = @; 


if y = 8, ¥, = /,, 


str 


= ¥, { = 0. 


The function V,(¢) has to satisfy the following boundary conditions: 


v,0) = ¥,() = ¥,,,0) = ¥,,.() = 0. (13) 
Then 
cosh uw, cos uw, — 1 = 0, (14) 
_ ©0sh u, — COS mu, 5 
~ sinh p, — sin yp,’ (15) 


The values of 8, and yu, are presented in Table I. Since the characteristic functions are 


TABLE I 
Values of B, and p, 








3 Me ue 

l 0.98250 +.73004 500. 564 
2 1.00078 7.85320 3,803 .537 
3 0.99997 10.99561 14,617.630 

j 1.00000 14.13717 39 943.799 
5 1.00000 17. 27876 89, 135.407 
6 1.00000 20. 42035 173, 881.316 
> 6 1.00000 (2v + 1) r/2 (2v + 1)4 r*/16 
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of the form of Eq. (12), we obtain 


| 
| 


X, exe = X,, Y, we» = i !Y., (16) 
a b 
ise j! if Pp =f 
| X,X, dx = aé, , where 6; = ) (17) 
| 0 if pwr 


and 


| Y.Y, dy = bé;,, where 5; = ) (18) 
0 if q # 8. 


Using the abbreviation 


. \! if q=s | 
| 


H,, =a] X,X,.2d2%, He=b| YY. dy, | 
Jo Se 
knowing that 
| Vv be dé = | VW; : — | Y, WV; .: dé 
0 “0 | 
and using boundary condition Eq. (13), we have 
| 
| UM, de = —] Wey .¢ dé. (19) 
« 0 
TABLE II 
Values o & e,eV; é lé 
l 2 3 } 5 6 
1 12. 30262 0 9 73079 0 —7.61544 0 
2 0 416.05012 0 —17.12892 0 15. 19457 
3 —Q 73079 0 98 . 90480 0 —24.34987 0 
t 0 -17.12892 0 171. 58566 0 31.27645 
5 7.61544 0 —24.34987 0 263 .99798 0 
6 0 15. 19457 0 —31.27645 0 376. 15008 


Numerical values for the integrals, Eq. (19), are presented in Table II. Hence 


aa 


H,, = —a | Dnothe.o G0; (20) 


/0 


Hye = —b | YosYou dy. (21) 


“0 





7 
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Now Eq. (11) becomes 


na 


m n f at 
Fe eet § RR T+ ERA Pin 


$+ 2X,Xi0eV Yew) dedy—dr] | XX,YY, de ay} =0 


or 
’ n -a 4 sb a ob 4 
2. } a,,‘ | ps XX, dx | Y,Y, dy + | X,X, dx | a YY, dy 
+ 2 [ Le a [ YY... dy — dX [ X,X, dx [ Fane ay} = 0. (22) 


After substitution of the corresponding values, and multiplication of each term by a” 


we obtain 


pb y a {| (4 4- Ma base +2 = HH | -- na'vase”| = @, (23) 
\L \a b ab 
where 
6° = | X,X,dx] YY, dy. (24) 
Equation (23) may be rewritten as 
— - l i) 1 rs) (rs ~ 
} > oe A| (E48 + 73 ut)as ee r H, n..| - nas} = 0, (25) 
= LL \¢ ) 
in which 
le ii ‘ 
A = dab = Fp E + ( ; |. (26) 
Let 
h 3 
Eo =  ( ).! 4 (<) us (8? +25 HH « (27) 
a b b 
Then Eq. (25) becomes 
> Dd 4..1E” — A8,”] = 0, (28) 
where 
jl for p=r and q=s 
hag = 
lo if either pwr or q#S8 
and 
os = ar « a 
ES up + (2) wu, +2 vile for p=r and q =8, 
a ) 
(29) 
Eee 2 : if either pr or q #68. 
) 
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Equation (28) is a system of mn linear homogeneous equations in mn unknowns 4a,, 


) 


yr =1,2,---ms=1,2 


n). Since a,, 


cannot all be zero, the determinant of the 


system must be zero. This leads, in general, to an algebraic equation of degree mn in A. 


From Eq. (26) it follows that 


AA 


Then for any shape ratio ¢ = b 


If k = 0, then 


- | AD _ k )'] 
* pha’b \2 ph 


a, one obtains 





7 | AD _ ( k y |" 
~ L phoa* \2ph/ _! 


w phoa 


A=) 


(30) 


(32) 


which is the eigenvalue parameter for the plate with fixed boundary conditions. From 
Eq. (30) it can be seen that the frequency decreases with an increase in the value of the 
damping factor k. The zero value of w is obtained for 


pad ap a] . 


a 


(33) 


The frequencies are calculated for the first three modes of a square plate fixed along 
each edge, and the results are presented in Table IV. 


rABLE III 
Coe ts for \ of Damped Squar 

r ; i E E E 

13 22 
1 1 1.30 R4 0 — 239 .43 0 0 
1 2 0 7.18 0 0 0 
1 3 2 i 17 1.77 0 0 
2 1 ) ( 0 137.18 0 
2 2 0 0 11,848.30 
2 2 rs) 0 0 896.21 0 27 
3 1 — 239 .43 0 189.38 0 0 
3 2 0 896.2 0 0 0 
3 189 .38 ( 1,924.84 0 0 

TABLE I\ 
The } of A!/2 of Vibration of 


Mode 
All 
D. Young 


Table III represents the coefficients £ 
Using only the first two terms of the series for E,," 


\2 


ra AD _ ( k ) 
~ L phea* 2 ph 


Ist 2nd 
36.11 73.43 
35.99 (73.41 


pa 


e Plate 
E E I i 
2 31 2 
0 — 239.43 0 189.38 
0 0 —896 .21 0 
0 189.38 0 924.54 
896 .21 0 0 0 
U 0 0 0 
7,530.32 0 0 0 
0 17,551.77 0 124.84 
0 0 27 . 530 .32 0 
0 -1,924.84 0 18,799.58 


Se uare Plate Fixe l A Ong Each Edae 
1 


3rd 
108 .85 
108 . 27) 


for the vibration of a clamped square plate. 


, we obtain approximations for 
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A’ which are compared with the D. Young solution [5] (given in parenthesis in Table 
[V). It can be seen that a good agreement between the results is obtained. 






































hhh Lhhhly hhh 
ae 4 
A 4A ty 
4 — 77 - ‘4 
y, y ae 4 
7 tLhhh 7CIIA77 7 
lst 2nd 3rd 


Fic. 2. Nodal Line for Square Plate 


Conclusion. The problem of free vibrations of a rectangular plate fixed along each 
edge and having internal viscous damping is presented and solved by means of general- 
ized Galerkin’s method. By choosing the function ¥(£) to satisfy the boundary con- 
ditions, Galerkin’s method can be applied to a plate with any type of support, even 
when damping forces are present. 

The influence of the damping factor on the natural frequency is given by Eq. (30). 
It can be seen that the natural frequency decreases with increase in the k-values. 
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ON THE DEFORMATION OF ELASTIC SHELLS OF REVOLUTION* 
BY 
P. M. NAGHDI ann C. NEVIN De SILVA 


Department of Engineering Mechanics, University of Michigan 


1. Introduction. In a recent paper, E. Reissner [1] formulated a theory for finite 
deformation of elastic isotropic shells of revolution where the theory of small deforma- 
tion (linear theory) is also discussed. In the present note, there is derived a single complex 
differential equation for small deformation of shells of revolution which is valid for 
uniform thickness, as well as for a large class of variable thickness. 


z 

















Fic. 1 


With the use of cylindrical coordinates r, 6, z, the parametric equation of the middle 
surface of the shell (see Fig. 1) may be represented by 


r= re, z = 2(é). (1.1) 


Denoting by ¢ the inclination of the tangent to the meridian of the shell, then 


, 


r’ = a cos ¢, z’ = asin ®q, (1.2) 
where 
a = [(r’)? + (2’)?]'” (1.3) 
and prime denotes differentiation with respect to &. 
We note for future reference that the principal radii of curvature r,; and r, are, re- 


teceived January 18, 1954. 
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spectively, the radius of the curvature of the curve generating the middle surface and 
the length of the normal intercepted between this curve (generating curve) and the axis 
of rotation. It follows from the geometry of the middle surface that 
r=fr,an ¢. (1.4) 
: , Ng and Q, and the stress couples M; and M, , acting on an 
element of the shell, are shown in Fig. 1. Also, as in [1], it is convenient to introduce 
‘thorizontal” and “vertical” stress resultants, H and V, given by 
aN; = rH +2’V, aQ = —2/H +7’V. (1.5) 


We now record the basic equations of the small deflection theory of elastic shells of 


The stress resultants NV; 


revolution with axisymmetric loading, as given by Reissner in [1]. 


rV = — [rap dé, 


aN, = (rH)’ + rapa ; 
rN; = (rH) cos¢@ + (rV) sin ¢, 


—(rH) sin ¢@ + (rV) cosa, 


Di. 2 1.6) 
M;=- E + y- al, : 
a [ 


© 
i] 


ws / F (N yN»o) — a | dé, 


where 8 is the negative change in @ due to deformation; u and w are the components of 
displacement in the radial and axial directions; p, and py denote the components of 
load intensity in r and z directions; h is the thickness of the shell; and 
ryp3 
, , Ith “s 
Cc = Eh, Dp = 5 ere 
IZ1 — » 
E and v being Young’s modulus and Poisson’s ratio, respectively. 
With 6 and rH as basic variables, proper elimination between Eqs. (1.6), differential 

equations of equilibrium and compatibility, leads to the following two second-order 
differential equations: 


TT) (rD/a)’ _, ; ‘a (r’D/e r] z’ yr’ , 

B T =— 5 | ) —y - Ir : rH) = - rV). 1.8) 

rD/a \T (rD/a) (TL) /a) (rD/a) 

. T al j | (r’ " / c oy | 
(rH) 4 Hy —| (2) +02 lon 
] al 4 / / al } 
— [r’e’  (2’/a)’ , pee 

— — 8 =x a +p mn HV) + yp — (rh (1.9) 

r al ) SG, (r/al T 


y , ,- 
(r/al oe 
oll a a = Po ee) — (apy) 
(r/aC) y 
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2. Normal form of the differential equations. Substitution of the quantities C 
and D from (1.7) into (1.8) and (1.9) and rearrangement of terms result in 


i) +f Clad 4g OHO Dy 4, am (I), _ am (ha)(mr¥) 
L,(6) + | (r/a) +s (rh) iF a Toho \h ont ¢, 


roto \h /\ Eh? 
1 r’ fa)’ | 9 Wh’) am ("*) 
Ly) — | clay  * ly rato \h /? 


(2.2) 
h”? (r’h’) (r/o)! h’ _ 
pe 9 ME A {a7 =a Z 
7 | h + Gh) * (Ja) h ly Z 





(2.1) 





Eh?’ 
where Z denotes the right-hand side of (1.9), hy is the value of h at some reference section 
(say & = &), and 


(2.3) 
_ mi _ neve . rt”? 
= he? m = [12(1 ee 


ae m 


rh aa Qu’ f(é), 


here » is constant and it is to be noted that f(¢) is independent of the thickness h(é). 
Then, multiplication throughout (2.1) and (2.2) by {h[hof(é)]"*} results in 


L(8) + vAB + Quy = F, 


(2.5) 
Liv) —-(wr- dy 2Qu’B = G, (2.6) 
where 


— [ag Ler ed!  30rh| 
d k re | {ue + 


rh 


~ 

lI 

i) 
a 
=| = 


re | {HE + ay CHD 4 ey Mh 


9 
LA (rh) (r/a) h (2.7) 


a ho | 
C= ER k f@) | 4. 


Introducing the complex function 


U=8+ikv; i=vV-1 (2.8) 








372 P. M. NAGHDI AND C. NEVIN DE SILVA [Vol. XII, No. 4 


where k, an arbitrary function of £ is to be determined, the differential equations (2.5) 


and (2.6) may be combined to read 


7.[ vr ) 
; E - ae -5,) ” | ) 

wena Be 
2 nA) 

au 


J ho . ‘| | y’ (r/a)’ | 
. “—_ ee 9 S PY ibn ‘ J? 
= iE se | hk’ +12 J — (r Je) +3 h k’ ey 
+ (F + 2kG). 
Taking k in the form 


f 2)1/2 
oi ental (on _ 8) ne E (nn _ ‘) (2.10) 
; 2u 2/ 5, 


Qu? Z 
with the restriction (the implication of this restriction will be discussed later) that 
k’ =k” =0, (2.11) 
(2.9) transforms into 
vr 


L(V) = aul ik att 5)U + (F + ik). (2.12) 
2u 


By putting the last complex differential equation in the form 


. vr 
2u*ha( k + iz :) 
4 (CU) — } . = = ~ = - f => 0 y" ik ; 
by hf ‘(é) ( hf'(é) (F + kG) 


and observing that the coefficient of U’, resulting from the application of the operator 


L, , defined by (2.3), is 
(r/a)’ h’ 
R= kt bam & 
(r/a) h 


exp E / R ae = en(2) : 


then, with the aid of the transformation 


W = (;-) (2) ‘U, (2.13) 


WwW" + [2° vw) + AOIW = E ‘| f®[F + ikG] 2.14) 
(y Ql 


and that 


we obtain 


which is the normal form of (2.12), and where 


v= (e+e )(O)re 


1 (r/a)’’ 1 | (r/a)’ |’ (“) 3 (r/a)’ h’ 3h” 3 (“) . 
A = --- + - i — |— — = . oes — = 3 2.1: 
Y | | vs 2 (r/a) h 2h $ \h 15) 


2 (r/a) 4 L(r/a) 
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We now return to (2.12) and observe that condition (2.11) is fulfilled only if k isa 
constant, and this may be achieved by proper choice of \ and 6. In particular, we note 
the following two cases. 

(a) For shells of variable thickness and with reference to differential equation (2.12), 
the condition (2.11) is satisfied, provided (vA — 6/2) is constant. Thus, by (2.7), 


(:), oe (z)' . | (Ee + ea = (‘)xre). (2. 


This equation is directly integrable and its solution is 


h = Ky’ [+ wonal ff (2)r at d§ + ear’ +c,” [ roa dé, (2.17) 


where ¢, , C2 are constants of integration. It may be noted that setting K = 0 corresponds 
to the vanishing value of (vA — 6/2) or, by (2.10), tok = 1. 
(b) For shells of uniform thickness, 6 vanishes identically and we have 


= pre C2], ditties a 


Clearly, \ is a function of — and its form is determined by the geometry of the middle 
surface. However, for numerous shell configurations \ and, by (2.10), k is either exactly 
or very nearly a constant. 

It should be mentioned that whenever the radius of curvature of the generating 
curve r, is a constant (r, may be a function of £), then, with proper choice of & (§ = ¢) 
and by (1.2), (2.4) and (2.18), A, and thus k, are in fact constant. The cases of conical 
shell and toroidal shell, treated recently by Clark [2], are included in this class. 

3. Remarks on the solution of equation (2.14). Particular solutions of Eq. (2.14) 
may be obtained approximately by the membrane theory of shells given by Hildebrand 
in [3] or sometimes by a more recent method developed by Clark and Reissner [4]. In 
this section, we shall discuss briefly the homogeneous solution of equation (2.14). 

If v’ and A are suitably regular over a finite interval of the axis and furthermore, 
if ©* is bounded from zero everywhere within this interval, then the classical method of 
asymptotic integration leads at once to the solution 


— 


6) 


W = v-'”*{ Ae~* + Be} (3.1) 


which, by means of well-known relations for Bessel functions, may also be written as 
W = & [vac passat + BJ,,(n)}, (3.2) 
where A and B are constants and 
n= (2°) | ¥ de. (3.3) 


According to Langer [5], if ¥° is not bounded from zero everywhere within the interval 
in question but vanishes to the degree n at some point £ within this interval, then (3.2) 
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may be generalized to 
af 1/2 
W = v ; | V(¢) dg {AJ_, n+2() + BJ yn+2(n)} (3.4) 


which is valid at & . 


In (3.4 n = (20°)'"p | W(g) dg. 


“so 


An appropriate form of this solution was recently employed for toroidal shells by Clark 
[2]. 

If, moreover, the coefficient function A of the differential equation has a pole of 
second order at & , then by a more recent method of asymptotic integration developed 
by Langer [6], a representation of W in terms of Bessel functions is again possible and is 
valid at & . It may be mentioned that application of this method yields a solution for 
ellipsoidal shells of revolution which is valid at the apex (where a pole of second order 
occurs); this will be given on another occasion. 

4. Acknowledgement. The authors are indebted to a referee for helpful criticism. 
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DISTRIBUTION OF THE EXTREME VALUES OF THE SUM OF n 
SINE WAVES PHASED AT RANDOM* 


BY 
8. O. RICE 


Bell Telephone Laboratories, Inc., New York 


1. Introduction. The statistical behavior of the sum of n sine waves phased at 
random has been studied in connection with a number of technical problems. These 
include radio wave fading and overloading in multichannel telephony.’ 

When the sine waves are of unit amplitude their sum may be written as 


n 
z= >> COs Yn (1.1) 
m=1 
where ¢, , ¢2, °** , ¢, are independent random angles, each distributed uniformly over 
the range —z to x. z cannot exceed n. When n — 2 < z < n the probability density 
of z may be expressed as a power series in (n — 2), as is shown in See. 2. 


There is a close relation between the distribution of z and the problem of the random 
walk in two dimensions, and the two are often treated together. Several equations 
connecting them are given in Sec. 3. In Sec. 4 the results of Sec. 2 are used to obtain 
the first few terms in a series for the distribution of the extreme values in the random 
walk problem. 

When n is large the central portion of the distribution for z approaches a normal law. 
In Sec. 5 an attempt is made to obtain an approximation to the distribution over the 
entire range of z by interpolating between the normal law result for small z and the 
results of Sec. 2 which hold for extreme values of z. The work is carried out first for the 
random walk distribution and then translated to the z distribution. This procedure is 
used because the random walk distribution seems to be better suited to our method of 
interpolation than does the z distribution. Figure 1 is associated with the interpolation 
between the results given by Pearson’ and Rayleigh’ for the random walk and those of 
Sec. 4. 

I wish to express my thanks for the many helpful suggestions concerning this paper 
which I have received from Mr. John Riordan and others. 

2. Series for the probability density of z when z is near n. Let q,(z) denote the 
probability density of the random variable z defined by (1.1). Then, when n — 1 < 


*Received Jan. 27, 1954. 

See, for instance, W. R. Bennett, Distribution of the sum of randomly phased components, Quarterly 
Appl. Math. 5, 385-393 (Jan. 1948). References to earlier work will be found in this paper. Mention 
should also be made of papers by F. Horner, Phil. Mag. (7) 37, 145-162 (1946), and R. D. Lord, Phil. 
Mag. (7) 39, 66-71 (1948). The second paper gives our equation (3.5). 

2Drapers’ Co. Research Memoirs Biometric Series III. Math. Contributions to the theory of evolution— 
XV. A mathematical theory of random migration, Karl Pearson assisted by John Blakeman, London 
(1906). 

tayleigh, Phil. Mag. 10, 73 (1880) and Phil. Mag. 37, 321-347 (1919). 
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an 
Qa+1(2) = qi(z — u)q,(u) du 
v¥z—-1 
] si [ 21-1/2 9 
=- il — (#@ — w) ] In(U) du. (2.1) 
FT Js-} 


The change of variables 


carries (2.1) into 


rere (x/2)'* fp’, concrete 
Qavi(n + I SS | [vl — vx/2))-°'’"g,[n — (1 — v)a] dv (2.2) 
us 0 
which holds when 0 < 2 < 2. When we place the assumed expansion 


qn — 2) = A,’ be iE + >a (x/4)*/n(n + 2) «+> n+ 2k — 2) | (s 
k=1 


in (2.2) and use 


bo 
OO 


qg(z) = q@l—2) =r [2x —x]’” (2.4 
we obtain 
A, = (Gey ""/Tin/2Z (2.5) 


and a set of recurrence relations, the lth of which is 


Onei,1 = >, Gp,1-10,/K!, l1=0,1,2,-:-, (2.6) 


where a@,. = 1, a = 1 and 
ar, = 1-3-5 --- (2k -— 1] 
The series in (2.3) converges for | x | < 2 as may be seen by substituting ¢,(n — x) = 
—l+n/2 ° l/o —-1/2 + a « ° ° ° ° 
x jAt), ii(@) = « (2 — 2) in (2.2). An integral is obtained which may be 
used to show in succession that f.(x), f;(x), --- f,() are analytic functions of x inside 
the circle | x | = R, R < 2, in the complex z-plane. 
Equations (2.6) and a,, = a;/k! lead to 
2 Gus = 1 D>, als’ ‘| ‘ (2.7) 
& Li 


This is merely a formal result because the series on the right does not converge. 
It is interesting to note that (2.7) fits in with some heuristic manipulation of the 
integral 


Rk aneki ; 
q(2)=s- | e7**[Jo(O]" dt. (2.8) 
Thus if we raise the asymptotic expression 


x a 1’? 2 a y it 
Am " ‘ 
Jt) ~ ) e am | 7cpaAmei7a T DAMtI7a (2.9) 
(22t) (—720) 


a= 6omi4 
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(where —r < argt < mand arg (—7) = —7/2) to the nth power we obtain the sum 
of terms of the form exp i(n — 2l)t times a series in 1/t with 1 = 0, 1, --- n. Let the 
i-plane be cut along the negative real ¢ axis and let the limits of integration in (2.8) be 
—i + o instead of +. When we substitute (2.9) in (2.8), assume n — 2 < z <n, 


and use 


: iat 294 
~-t=-e— > 
| i-e (it)” P'(v) 
0, z<0 


only the terms multiplied by exp [7(n — z)t] (corresponding to 1 = 0, x = n — 2) contri- 
bute to the value of the integral. Furthermore, they lead to the same series in x for 
q,(n — x) as does (2.7). If this procedure could be justified and generalized it might 
lead to an expression for (2.8) which would supplement the one obtained by the method 
used by W. R. Bennett.’ 

The coefficients a,, may be expressed in terms of Bell’s Y polynomials.* Thus, 
multiplying both sides of (2.7) by t"/n! and summing n from 0 to © shows that a,, is 


the coefficient of t"s*/n! 


© i @ k 
F 2 § t § , 2 2 
[exp al exp >» (1a? = )] = ¢ zz. kl Y,(ta; , tag, -*-, tai) 


1 j! k=0 


in 


where the Y,(y, , yo +++ y) are Bell’s polynomials: 
Y, = 1, Y,=%1, Y,=yt+yi, Y; = ys + 3yoy1 +i, °°? 


Consequently, 


! 
— Y,(tai , tas, +--+ , taj), (2.10) 


k! 


aa = 


where, after writing the right hand side as a polynomial in é, ¢” is replaced by 1/(n — m)!. 
We obtain in this way 


, (2/20)? [ 42 (n+ 8)x? . (n? + 24n + 200)2° | (2.11) 
0,n — 2X) = vn ; = eee é- 
‘ 2nT(4n) 1° 32(n + 2) " 384(n + 2)(n + 4) 


where the series converges when 0 < x < 2. 
= 


The probability V(Z) that | z| > E is given by 


an an—E 
W(E)=2] g@daz=2] aln—a)ae (2.12) 
“E “0 
and when n — 2 < E < n this leads to 
7 2 n-Er ¥ n(n — E) , nn+ 8)(r — E)’ | . 
FA) = r[a(n + 2)] | Qr | E + 4(n + 2) + 32(n + 2)(n + 4) - (2.13) 


*Actually what is used is the slightly more general version of these polynomials given by John 
tiordan, Derivatives of composite functions, Bull. Amer. Math. Soc. 52, 664-667 (1946). 
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3. Relation between z and the problem of the random walk. The random variable z, 
defined as the sum of n cosines by (1.1), may be regarded as the projection of the resultant 
r (of n unit random vectors) on the z-axis. Hence, we may write z = r cos 6 where @ is 
a random. angle distributed uniformly over the interval (0, 27). The probability p,(r) dr 
that the length of the resultant lies between r and r + dr is given by the random walk 
distribution when the n elementary linear walks are of unit length each. We shall use 
®,(r) to denote the probability that the resultant equals or exceeds r. Then the connec- 
tion between z and r leads to the following relations between the probability functions: 





9 an 
V,(E) = - | p,(r) are cos (E/r) dr, (3.1) 
TWJE 
9 a an : - 
¥,(E) = = | r'@? — E)-8,0) ar, (3.2) 
WT JE 
ar 3 ere : 
(7) = —r>— | (@ - )'*¥,(£) dE, (3.3) 
dr J, 
@(r) =r° | E”(E? — 1)? (BE 'W,(£) — Vi(E)] dE, (3.4) 
, gee en a ae — 
qn\2%) = — | oe. =e PAY) ar, (3.2) 
pr) = —2r | ( — r’)~"/*1dq,(2) dz| dz. (3.6) 


In these equations FE, 7 and z are assumed to be less than n and p,(r) is to be inter- 
preted as an impulse function. In (3.6) n must exceed two but this causes no difficulty 


since it is known that 


> /d 2) —1/2 i. ‘ 27) 
PAT) = (2 m)(4 — r’) P —2<r< 2. (3.7) 
In going from (3.1) to (3.2) we have integrated by parts. Settingn” — ro =i,n°-K° = 2 


in (3.2) converts it to a special case of Abel’s integral equation* whose solution gives 
(3.3). The remaining equations are obtained by the same kind of analysis. 

4. Random walk distributions when r is nearn. Whenn > 2andn —2<r< vn, 
substitution of the expression (2.11) for q,(z) in the integral (3.6) for the probability 


density p,(r) of r gives 


1/2 \ (n-3) /2 
n n-—-r?r (n — 1)(n — r) 
palt) = oT [3(n — 1)] (' Qn ) [ paces 4n 


2 j —s —_—~ 2 
ptm mes] ay 
32n 
When n = 2 the method fails, but it is not difficult to show from (8.7) that (4.1) also 
holds in this case. Expression (4.1) may also be obtained from the recurrence relation 
for p,(r) by a method similar to that used in Sec. 2 to obtain g,(n — x). Pearson” has 


*See, for example, Whittaker and Watson, Modern analysis, 4th ed., Cambridge, 1927, p. 229. 
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given, essentially, the leading term in (4.1) and has given one or two more terms for 


n = 8, 4, 5, 6. 
Integrating (4.1) termwise gives 


b (7) [ y - |" - ry" aE 4 (n — 1)*(n — 1) 
AD] Pe OP Tha + DIL Qe in(n + 1) 
(n — 1)(n? + 4n — 9) . | 
+ —~,—____—“ (n — r)° vee 2 
: 32n7(n + 3) " a + (4.2) 
which holds forn — 2 <r < n. 


5. Approximations for #,(r) and ¥,(£). Numerical values of the various proba- 
bility densities and distributions have been given by Pearson’, Slack*, Bennett’ and 
others for values of n up to 10 (and somewhat beyond in certain cases.) The values of 
®,(r) and V,(2) given by Bennett were computed from series which converge for all 
values of r and E between zero and n. In this section we shall consider a method of 
estimating values of ®,(r) and ¥,(/) which involves only a small amount of calculation 
but which, of course, lacks the accuracy of the computations mentioned above. 


When n is large, but E and r of 
®,(r) 
W,(E) 


erf (x) = 


moderate size, it is known that 
(5.1) 


(5.2) 


exp (—r’/n), 


1 — erf (E/n'”), 


ios [ exp (—2x) dx. 


“0 


Our method of estimation is roughly equivalent to interpolating between values given 
by these formulas and those given by the formulas of Sec. 2 and 4. We shall first work 
with the random walk distribution. 

Equation (5.1) suggests that we introduce a function y of r and n defined by 


—ny 


=e 


®,(r) ( 5 3) 


or 


1 . 
y= — log ®,(r). (5.4) 
Comparison of (5.1) and (5.3) shows that, for small values of r/n, y = (r/n)* and hence 
in this case y depends ‘‘much more” on r/n than on n. The same is true when r/n is 
nearly unity since Eq. (4.2) for ®,(r) leads to (assuming n large so that the series may be 
approximated by exp [(n — r)/4]) 
n— 1 r 1 r 1 T 
y = ———— lo (1 -2)-3(1-4)+ log — 5.5 
2n 5 n a n 2 Be? (5.5) 
and this again tends to become a function of r/n only as n > o~. 


In order to test the dependence of y on r/n, Bennett’s values of ®,(r) for n = 6 and 
10 were used to compute y from (5.4). The results are plotted as the “exact” values of y 


‘Margaret Slack, Zhe probability distributions of sinusoidal oscillations combined in random phase, 
J.I.E.E., Pt. III, 93, 76-86 (1946). 
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(indicated by the small triangles and circles) shown in Fig. 1. It is seen that the two 
sets of values tend to follow a common curve. 
The dashed curves in Fig. 1 were computed from the approximation (5.5) forn = 6, 10 

















and «. Although (5.5) is valid only for r/n = 1, it yields values of y which are fairly 
' 8 4 
» / 
____ EQUATION (5.5) j 
VALID FOR F/N=! / ‘d 
1.7 } j } 4 } ig 4 
EQUATION (5.6 y Fd 
aaa 4 
VALID FOR r/N SMALL y ie 
1.6 avy Z 7 r 
a 
1.3 
1.2 
3 
0.9 
ao 
0.4 
0.98 
Fic. 1. The triangles (n = 6) and the circles (n = 10) show exact values of y = (—1/n) log, ®, (r) 
®, (r) is the probability that the resultant of n random two-dimensional unit vectors is longer than r. 
The curves show the approximations (5.5) and (5.6) for n = 6, 10, and ~. 


close to the exact values even for values of r/n as small as 0.4. The solid curves were 
computed from 


FP YI EPHOL SO 6 
= sa + - a a = —_-_—_- = —_ — _— s ‘ns 5.6 
y (£) ' 4\n i 36 \n nL2 \n ” 3 \n 12n? \n (5.6) 


which holds only for small values of r/n and large values of n. Expression (5.6) is obtained 
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from 
FE 2f. (6n — 11)f; 
0 «11-2. 24 St 5 ... 5.7 
| 2n 3n* 8n* ; (5.7) 
where x r-/n and 
— m > m(m— 1) 3 
f aes fe ee —m ~ pe A = aE 
te x,F,(—m; 2; x) x+ Tike 913! x + : 
i’',() being a confluent hypergeometric function. Expression (5.7) may be obtained 


by integrating a result given by Pearson’ and, later, by Rayleigh*. Pearson’s formula 
suggests that the next two terms in (5.7) are + (50n — 57) f,/15n* — (270n’ — 2125n + 
1892) f,/144n’. 

Thus, a rough idea of how ®,(r) behaves for all values of r/n and a particular value 
of n may be obtained by (i) computing approximations to y from (5.5) and (5.6),* 
(ii) plotting them on semi-log paper as shown in Fig. 1, (iii) joining the two portions 
by a smooth curve, and (iv) using these approximate values of y to compute ®,(r) 
from (5.3). 

When we turn to ¥,(/), the procedure used in dealing with ®,(r) suggests that a 


new function y’ of E and n be defined by 
V,(E) = 1 — erf [(ny’)'”’]. (5.8) 


However, it is found that the expression for y’ when E/n = 1 does not have the sim- 
plicity of its analogue (5.5). Instead of following this line of thought further, we note 
that the analogy between y and y’ suggests that y’ should not differ greatly from the 
function obtained by replacing r by E in y. That the difference is small may be verified 
by comparing Bennett’s exact values of ¥,o(2) with approximate ones obtained from 
(5.8). In using (5.8), the values of y’ are taken to be those of y as computed from (5.4) 
and Bennett’s exact values of ®,.(r). The exact and approximate values are shown in 


the follow ing table for three values of E. 


Ki exact V,.(/) approx. V,0(£) 
l .6604 .6624 

5 02337 02440 
8 .00009 00011 


Thus, once we have obtained approximate values of y from curves of the type shown 
in Fig. 1, approximate values of #,(r) and ¥,(2) may be obtained readily from (5.3) 
and (5.8) (with y in place of y’). 

*Instead of (5.6) one may use (5.7) and (5.4). This gives more accurate values of y at the cost of 


more computation. 
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THE MECHANICS OF THE RIJKE TUBE* 


BY 
G. F. CARRIER 
Harvard University 


1. Introduction. An interesting phenomenon, the analysis of which has received only 
scattered attention, is the thermal-acoustic oscillation discovered by Rijke in 1859. 
It was noted, under certain circumstances, that superimposed on the anticipated steady 
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Fic. 1.1. Open ended cylindrical tube containing heated ribbon at x2, = cL. 


convective flow in the apparatus indicated in Fig. (1.1) was an acoustic oscillation. 
The frequency of the phenomenon is essentially that of the first free oscillation mode 
of the pipe (i.e. the wave length is 2L’)," but the occurrence (or lack thereof) and the 
intensity depend on a number of parameters. 

The analysis of the effect of these parameters which is presented here proceeds from 
the fundamental conservation laws with an investigation of the response of the heater 
to a fluctuating velocity, an analysis of the non-isentropic wave propagation in the pipe 
with its attendant viscous losses and end radiation, and, from the interaction of these, 
the deduction of the complex eigenfrequency of the tube. In a more qualitative manner, 
the consistent appearance of harmonics is explained and the role these waves play in 
determining the sustained intensity of the sound is discussed. 

The results are in excellent agreement with the observed physical facts, and it is 
believed that they, as well as the details of the analysis, may be useful in the future 
investigation of various combustion oscillation phenomena. 

2. The wave propagation in the tube. It is convenient to divide the analysis of this 
composite problem into various pieces. The description of the waves in the tube for 
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1’ is the length of the pipe corrected for end effects and non-uniform temperatures. 
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x < o and those for x > o are treated individually.” The boundary conditions at ¢ 
“joining” these waves are treated as though a disk shaped heat source, whose response 
to local velocity fluctuations is known, were located at o. However, in the detailed 
investigation of the heater response, a more realistic description of the flow past the 
heater will be adopted. 

In analyzing the oscillating flow downstream of the heater we must anticipate that 
it is non-isentropic and, in particular, that important viscous losses occur at the wall 
and that decaying temperature fluctuations are convected from the heater toward the 
exit. We assume that the steady flow upon which the oscillation is to be superimposed 
is known and introduce a small perturbation (in velocity, pressure, density, etc.) to 
represent the oscillatory phenomenon. The conservation laws (mass, momentum, and 
energy) are 


(put); + p’, = O, 2.4) 
Oui, + Uji.) Die = BUG + Gi) + Fe, 2.2) 
p’c(T’, + uiT’;) — (p’/p’)(pt + ulp’;) — kAT’ = 0, 268 
and we adopt the state equation 
p’ = p'R’'T’. 2.4 


Here, p’, p’, 7’, ui are the usual thermodynamic variables and the velocity, and /’; is 
the gravitational body force. It is convenient and appropriate in this analysis to treat 


Cp , C, , u, k, as though they were independent of the thermodynamic state. It is also 


convenient to introduce the notation” 
ui = v(x, 7r) + algrad g(x, r, 7) + curl ky¥(z, r, 7)], 
p’ = p(x, r)[1 + pla, r, 7)], 

p’ poz, Til + 7(2z,7, r) |, 


T’ = T(z, r{1 + Oz, 7, 7)I, 

z= 7,/L, r = (x. + 23)/L’, a” = ypo/po , r =at/L, 25 
where k is the unit vector perpendicular to both the x and r directions. The tensor 
notation differentiations of Eqs. (2.1) to (2.3) are taken with regard to the physical 
coordinates (the z;), but the vector notation differentiations of Eq. (2.5) pertain to the 
variables x, 7. In the definitions of x, r, 7, the constant length parameter L may be thought 
of as the length of the tube. It will be defined more precisely later. 

Equations (2.5) are now substituted into Eqs. (2.1) to (2.4), the terms containing 
no perturbation contributions are removed (since the steady terms are themselves 
solutions of the conservation equations), and the remaining equations for the pertur- 
bation quantities are linearized. Furthermore, the following simplification is adopted: 





*See Fig. (1.1). 
8This notation is appropriate for either the region upstream of o or that downstream of o. We shall 
distinguish these in Sec. 4 by subscripts 2 and 1 respectively. 
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+= 
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% , Po, Lo , Vo are replaced by appropriate constant ‘‘average’”’ values (in particular 
7 Po ; J 


v, = Ma) and », , v; are taken to vanish. The resulting perturbation equations are: 
n,, + Ag = 0, (2.6) 
o.t+y p= $e dg, (2.7) 
6,—(y — 1)n.. = fre AO, (2.8) 
¥., =€Ay, (2.9) 
p=nt 0. (2.10) 


Here we have taken the Prandtl number to be 3/4 for algebraic simplicity; € is defined 
by u/pal and y = c,/c, . The operator denoted by (_ _),, implies @( )/dr + Ma(_)/dx. 
That is, it is the convective time derivative.‘ 

These equations must be solved subject to the boundary conditions that, at the wall, 
the temperature and velocity fluctuations vanish, and at the plane x = a, the velocity, 
temperature, etc. are consistent with the heater behavior. At the exit, and at the inlet, 
a reflection coefficient must be adopted but we shall discuss that later. Equations (2.6) 
to (2.10) may be integrated in detail subject to such boundary conditions. The simplest 
presentation of the solution is obtained by eliminating 7, p, and 6, and by anticipating 
that the form of the solutions is exponential in x and 7. For example, ¥(z, r, r) = 
v(r) exp [i(ar — kx)]. Actually we must anticipate that there will be two modes of 
propagation; one which is essentially an acoustic wave and one which drifts with the 
stream. The former should consist principally of a contribution from ¢ where k* ~ a? 
(this implies that it propagates with a phase velocity which is essentially the acoustic 
speed), whereas the latter requires a less elementary description. The acoustic wave 
will consist of two parts; one which propagates in the downstream direction and a 
reflection from the tube end which propagates back upstream. 

Using the same letters to denote the functions of interest with the exponential 
dependence factored off, we obtain the following ordinary differential equations for 


yg and y after eliminating 7, p, and @. 
Lido) = 0, L(y) = 0, (2.11) 


where L, = A + 8°/(1 + 4iBye/3), L. = A — 3i8/4e, L; = A — iB/e, A = r-*0*/dr? + 
r'a/or — k’, and 6 = a — kM. 

For the upstream section analysis, it follows from our earlier remarks that 
8 = va — kM. The propagation modes are found by applying the homogeneous boundary 
conditions at the cylindrical wall and finding the eigenvalues k. If we anticipate that 
o ¢, + g with L,(¢,) = L.(¢.) = 0 for the acoustic mode, and further anticipate 
that in this mode ¢, and y will each display a boundary layer character’ whereas ¢, 
will essentially be a plane acoustic wave, the motion is described to terms of order ¢ by 


g(r) = AJ, [(8? — k’)'?r] + B exp [—(r. — 7378/40)" ] (2.12) 
‘When we consider the waves in the upstream region, it will again be convenient to define r = at/L 
giving a downstream value. The definition of ( ), , then will become »d/dr + Md/dxr where v = a;/az 
and M is evaluated for the upstream quantities. 
’The term boundary layer is used here in the more general sense as in [1]. 
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and 
vir) = C exp [—(r, — r)(iB/6)"”’), (2.13) 
where r, = R/L. Since @(r) = p — n = [—7tvB + (78) ‘Alg(r), the characteristic equation 
for k becomes fusing ¢.("7,) + ¥2("1) = ¥-("1) — ¢:(71) = Or,) = 0], 
—ik —ik (18/«)*’* 
k? — p’)/2 (378 /46)'”° ik = 0, (2.14) 
uy — 1)8 —3/de 0 
or (again to order e 
k? — B°)(iB/e)'”” = 2[k? + (y — 1)(4/3)'’"B"] 
and 
; (y — 1)(4/3)'” | ,.n1/2 
k = +a 1 + - . “ (é 18)' : ‘ (2.15) 
ry 


The imaginary part of this number defines the decay rate of the wave which is associated 
with the dissipation and heat conduction near the wall. Its order of magnitude for many 
experiments (in particular, ours), renders it an important source of energy loss. 

The propagation of the other “drifting’”’ wave cannot be expressed in as elementary 
a form. However, we may obtain a reasonably simple representation of this wave if we 
anticipate the nature of the result. It is clear that the wave would be a plane wave drifting 
with the stream if it were not for the cylindrical wall boundary conditions. The effect 
is an encroachment on this wave of such a nature that one cannot 


yi 


of the walls, however, 
expect a simple product type representation corresponding to the acoustic wave. If 
one attempts the product type representation, an infinite series of such products is 
required. For an efficient representation, one should formulate an initial value problem 
in x such that the temperature has a uniform value (no dependence on r) at x = o but 


such that v and 6 vanish at the wall. The solution to such a problem to the order of accu- 


racy we require is given by y 0 and 


af ; 
ear =—] ~ e"K[t(a? + S*)'”] dS, (2.16) 


3M /8e and ¢ = (¢ + 3ia/4e)'’~. A more accurate but no more useful (cor- 


where i 
solution would again involve boundary layer solutions in y and ¢ 


rections of order ¢'”” 
in addition to a slightly modified ¢“”’. 

Equation (2.16) is valid only for small enough x so that ¢ * (x, 0) takes essentially 
; However, for those x violating this condition, the wave has 


the value ¢°’(z, —@). 
decayed so far as to be of no further interest. The critical observation concerning these 
eigenmodes of the tube is that the mode associated with yg provides velocity and 
pressure contributions which are very small compared with its temperature and density 
contributions. Conversely, the k, , k, , modes give contributions of the same order in 
each of these four (dimensionless) state variables. This implies that the velocity and 
pressure fluctuations are almost entirely contributed by the acoustic modes while the 
final mode merely picks up any discrepancy in the temperature condition at x = o+. 


We shall see this more clearly in Sec. 4. 
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The wave in the upstream portion of the tube contains only the acoustic modes. 
The other possible mode is again a downstream drifting wave which has zero “input 
in the upstream section of the tube and is therefore not present. The formulas and k 
values for this section are identical with those of the downstream section except that the 
acoustic velocity, Mach number, M, etc. must be the values for the upstream state 
instead of the downstream state. In particular a is replaced everywhere by va except in 
the exponential exp (ar). 

Before closing this section, we must discuss the relation between the outgoing (away 
from the heater) waves and the reflected returning waves. The reflection coefficient for 
an acoustic wave in an unflanged tube has been investigated by Levine and Schwinger 
= 0. We have shown, and shall report in detail elsewhere, that the 
following is true. If the plane axially directed wave incident on the exit end of a tube 
(from inside) has acoustic pressure p; = P exp (twf) at some point 2, the reflected wave 
has pressure p, PN’ exp (iwt) where N’(wR/a, M) = N[wR/a(1 — M”’)'”’]. The 
formula for the inlet end is p,/p; = (1 — M/1 + M)N[wR/a(1 — M”’)'”]. Here N 
is the reflection coefficient computed in [2] for 1 = 0. However, for S << 1, N(S) = 1 
S°/2 + , and arg N is such that the wave is apparently reflected with reflection 
coefficient | N | from a point .613R beyond the end of the tube. Since M is of order 10™ 
for our problem, the moving gas correction is negligible and we use the above results 
with S wr/a (but with differing values of ‘‘a” in the two sections of the tube). We 
now define L more precisely. In fact x, = 0 is the point .613R below the upstream end 
of the tube and L is the point .613R beyond the other end. We shall return to, and use 
these results when we have determined the conditions at x = o which are imposed by 
the heater. To facilitate this, we record, respectively, the downstream (of oc), and up- 
stream, acoustic wave forms for use in Sec. 4, replacing J) of argument whose order is 


{2] for the case 4 


10°° by unity. 


gy’? = A, exp [iar — tk,(x — 1)] — |! N, | exp [iar + ik,(x — 1)], (2.17) 
gy = A, exp (iar + tk,x) — | N, | exp (tar — ik.2), (2.18) 


with 


+. (y — 1)(4/3)'”" ” , 7 
— of Pe Be ee Ie io)? a ke — §kit. 


"s 
ke. : va ss + — 16/3) _ iva)” = ky — thy’, 
is 
N, = 1 — (a7,)’/2 + - 
N, = 1 — (var,)’/2 + - 
vy = a,/a, 


The k{ , ki’, separation is such that k{ is real for real a. 
To the order of accuracy to which it is sensible to work, the downstream fluctuating 
velocity and pressure are given by 
) ( qa) ° (1) ‘ 
u'” = ay”, p’? = —iayy" (2.19) 


‘In our problem S is of order 0.1. 
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and, correspondingly, 
uu’ = ay, , pp” = waye"’. (2.20) 

As we shall see, the explicit formulas for the contributions of g“’ will not be needed. 
3. The heater ‘‘response”. In this section’ we shall estimate the fluctuating heat 
release from a heated ribbon to a stream of fluid whose velocity is fluctuating in a known 
manner. We say “estimate” because certain simplifying approximations will be adopted 
to render the problem tractable. However, the errors introduced should be of the order 
of a few per cent. In particular, we shall consider the problem whose simplified geometry 
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Fic. 3.1. The geometry of the simplified heater response problem. The ribbon thickness is oriented 


the y direction; Od is its breadth. V represents the fluctuating flow. 


inh 


is described in Fig. 3.1. For a known flow of an incompressible fluid with constant thermal! 
conductivity and specific heat, the energy conservation equation has the form 


pe,(T, + v-grad 1’) — k AT’ = 0. (3.1) 


Since the velocity distribution near the plate is rather complicated, this equation would 
be difficult to treat in its present form. However, the analysis of a similar problem [3], 
that of determining the viscous flow past an obstacle at low Reynolds number’, has 
indicated that an excellent approximation to the result is obtained when one replaces 
v-grad 7’ by cUT! in Eq. (3.1). Here c is a number, 0 < c < 1, (.43 is the optimum 
value in the viscous flow problem), U the free stream velocity, x the coordinate oriented 
in the flow direction and T7’(z, y, t) the fluid temperature. We introduce the notation 
cU = up + we’, T’ Q(z, y) + A(x, ye’, and T, (the plate temperature) = JT + ée’’. 
We shall take 7 and 6 to be constant rather than functions of distance along the plate. 
This is not correct in detail but two items justify this choice. A precise analysis of the 
surface temperature of the wire is very messy and the correction to the heat release 


7Unfortunately, the number of symbols needed in this paper is so large that the quantities a, 8, 
and other notations of this section are not related to those labeled with the same letters in other sections, 
This should lead to no confusion since the equations are never directly combined. 

®The nature of this approximation is discussed in [3], [4], and its verification noted in [4]. The present 
problem is very closely related to these and it is to be anticipated that the results obtained for the macro- 
scopic quantities should be highly accurate. The reader should note, in particular, that the predicted 
steady heat release is in agreement with the known information concerning the heat release from plates 


to fluids. 


20) 
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rate associated with this approximation is not an important contribution to the total 
release’. 

The differential equation may now be “separated” into two equations governing, 
respectively, the steady and fluctuating stream temperatures provided we assume that 
each fluctuating quantity is small compared to the corresponding steady quantity 


(e.g. 6/8 < 1). These equations take the form 
A6® — ad, = 0, (3.2) 


Aé — ab, — 186 = (aw/uw)9, , (3.3) 


II 


where the second order term in w@, has been omitted. In these equations a = pc,t)/k 
and B= pCc,W k. 
The boundary conditions are 


6— 0 as gm (x + iy)'? 3 @, O(z, 0) = T on the plate, 


and 


@—0 as gm (x + iy)'? 3 @, (x, 0) = 6 on the plate. 


We consider first the problem of determining @(z, y) when d(Fig. 3.1) is indefinitely 

large. This problem is readily solved by invoking Fourier Transforms and the Wiener- 

Hopf technique. Omitting the details of the analysis, the reader can readily verify that 
@(£, y) is the Fourier Transform *° with regard to x of the required O(z, y). 

@/T = - — exp {—| y | [#(a — 7é)]'”’}. 3.4 

it(a = i£) I ( | y | [2é( £)] ) ( ) 

The inversion integral is indented to pass beneath the origin (in the £ plane). In particular, 


the transform of @,/7' [(which will be needed to solve Eq. (3.3)] is 
@,/T = [a/(a — it)] exp {—| y | [t&(a — 7€)]'"}. (3.5) 


Other useful results are: 


0,(x, 0)/T = —(a/rx)'”, x>0 (3.6) 
A(x, 0)/T = erfe {(—az)'”’], x <0 (3.7) 

ane 
Q, = 2kl | 0,(x, 0) de = 4kIT(ad/x)"”. (3.8) 


This last quantity'’ will be associated with the steady heat output once we identify d 
with the true wire “breadth” and show that the contributions to the heat release rate 
of the finite plate is essentially that of the semi-infinite plate in the region 0 < x < d. 
To see this we note that the heat release rate of the problem just completed, which in 
the physical problem should be zero for x > d, is given by Eq. (3.6). We can now solve 
the problem where we postulate that 0//T on x > d, y = Ois|y| (a /mx)'’*/y and 0’ 

°This has been analyzed crudely but will not be reported here in detail. 

10@ is defined as |*., exp (—7&r) O(z, y) dz. 

11] is the length of the wire perpendicular to the plane of Fig. 3.1. 
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vanishes for x < d. This is not quite the boundary value problem the superposition of 
which on the preceding one gives the correct answer, but if it turns out in this problem 
that @/ is extremely small on x < 0, then the error in using the superposition of these 
two problems will be negligible. Actually, we anticipate that the change in heat release 
rate on 0 < x < d, as well as that on x < 0, will be negligibly small and so we replace 
the boundary condition on 6/ in our superposition problem by the requirement that 
6’/T be | y | (aw/d)'’’/y. If the heat release on x < d associated with the solution to 
this problem is small, that associated with the properly formulated problem will be 
smaller. The solution of this problem may be treated like the previous one and the 
essential part of the result is: 


(rd/a)'’O/(x, 0)/T erfe [a(d — x)]'* — [ra(d — x)]"“~ exp [—a(d — 2)], r<d. 


In our problem, a is of the order of 50 em~' and this analysis for the finite plate is justified 


provided d is of the order .1 cm or so. The heat release from the plate associated with 
0’ is implied by 
QO) = 2ki 6/(x, 0) de = IkT /(xad)'’” 
That is, for d of order .1 em or more, Qj is of order .05Q, and can be safely neglected. 
( .9.¢ 
We now turn to Eq. (3.3) and the boundary conditions on @(7, y), again for d in- 


definitely large. Another application of the same technique gives 


Thus Eq. (3.8) defines the steady component of the heat release rate. 


6(é, y)/T = {a/(a 1£) | exp j—| ¥ l7t(a — i€)]’ oie 
— [a/(a, — it)]'’? exp {-—| y | [(@E + a.)(a, — 7€)]'* }iaw/Bu 
+ [éa,/it(a, — 7€)T] exp {—| y | [((E + a.)(a, — it)”, 2.10 
where a, , —a. , are the roots of 2 — az — i8 = O(a, — a, as 8 ~ 0). In particular, 
6,(a, 0)/T [a(a’ /wx w/2 up B\[1 — exp (—a,2)] 
+ [(a,)'’* 6/T][(rx)~“’* exp (—ay.x) + ay” erf (a,x)'’*}, 
and 
e 2ia”’ “klTw vee oe a 
2kl | 6,(x, 0) da es [(a,)'’? erf (a,d)'”? + (ad)'*(e™*** — 1)] 
Jo Bu 
| S f / 1/27/ \1/2 4 f \1/2_ —a2d ‘ 
+ §kl(a,/az)'’*[(a,d + d) erf (a.d)'’* + (a.d/m)'”"e |. 3.11) 


Again, we can estimate the correction required by the fact that the semi-infinite plate 
treatment is only approximate but again this correction is smaller than is consistent 
with the accuracy of our basic model. 

In order to deduce the surface temperature fluctuation 6 we must analyze the flow 
of heat within the wire. A very crude analysis of this item will demonstrate that the 6 
contribution is negligible. This rough analysis is conducted as follows. We consider an 
infinitely long slab of thickness 2b with surface temperature ée'’. The differential equation 
to be solved (conservation of energy) is 


K.0,.’ = p.C.%, = 0 
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and the boundary conditions are 
o°"'(+b, 7) = be" 
Here, (w) merely denotes wire. 
The solution can readily be obtained and, in particular, 


6," '(b, r) = 26 tanh (QDd)e"’, (3.12) 
where 
Q = (iwp,C,/K,)'”’. 


The heat release rate for a wire of length / and breadth d is given by Q’ = 2K,1d6;"? (b), 
omitting the e’’. This must be equated to the heat release to the fluid as given by Eq. 
(3.12) in order to determine 6. For ribbons of thickness of order .03 cm and breadth 
.1 em (within rather large factors), 6 is so small that “‘the 6 term” of Eq. (3.11) is of 
order 10~° or less times “the w term.” 

Thus the fluctuating heat release from the ribbon is essentially given by 


—2ia® *kITw 


o ((a)'”? erf (ad)? + (d)'(e7*** — 1)]. (3.13) 


2kl | 6,(x, 0) dx = 
It should now be noted that (a,/a) becomes imaginary and small when 6/a* — 0 but 
has angle 7/4 and becomes large when a*/8 — 0. This implies that for large B/a’ the 
bracketed quantity in Eq. (3.11) is nearly a pure imaginary and the phase lag between 
the heat release rate fluctuation and the velocity fluctuation is small.’ Conversely, 
when 8/a° is small, the phase lag approaches 37/8. For the conventional Rijke tube 
experiment, the phase lag is of order 32/8. 

Before we abandon the discussion of the heater behavior, we should discuss quali- 
tatively the large amplitude behavior. To do this, consider first the extreme case where 
the free stream velocity is zero. In this case the response of the heater to a positive 
velocity fluctuation must be identical with its response to a negative velocity fluctuation. 
This is a direct consequence of the symmetry of the situation. Thus while the velocity 
fluctuation traverses one period, the heater is aware only of the absolute value of the 
velocity fluctuation and thus executes two cycles of heat release fluctuation. In other 
words, an input velocity fluctuation of frequency w excites a heat release fluctuation 
of frequency 2w. If the phenomenon is not linear, the higher harmonics will also be 
present. If one now estimates the expected behavior where w/w is not small and where 
u,) is not zero, he must conclude that the input signal at frequency w gives rise to a heat 
release fluctuation with both w, 2, and probably higher order components. The ratio 
of the 2w content to the w content of the output should be an increasing function of 
w/ uy . This information will be of use when we try to explain the harmonic content and 
sound level of the Rijke tube output. 

4. The matching conditions at the heater. The problem of combining the results of 
the foregoing sections is now a relatively straightforward matter. We merely write down 
the requirement that mass, momentum and energy be conserved across the heater. 
If we denote the quantities of interest at the heater inlet by the usual symbols with 


This corresponds to the situation conventionally encountered in hot-wire instrumentation problems. 
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subscript 2 and those at the exit with subscript 1, these laws take the form 
aX 
Pilly — Pollo = —| p, dx, t.1) 
+ X 
piu Pp Polls + Pp» — D/rR* — | (pu), dx, (4.2) 
+, 
puc,1, — pottec,T, = Q/rR* — h (oc-T); dz, 1.3) 


where Q is the heat input rate to the gas from the ribbon and D is the drag of the ribbon. 
It is tacitly assumed in these equations that 0 and X are in those regions where conduc- 
tion is not important, and it can readily be verified that the omitted kinetic energy 
terms are of higher order than we have consistently retained. In limits of the integrals 
in these equations, 0 and X are appropriately chosen points just upstream and down- 
stream of the heater. The integral term of Eq. (4.1) can be expected to be of the same 
order of magnitude as the other terms in that equation and a careful estimate of its 
value is required. As we shall see, the value of X is most readily chosen after we obtain 
the information which allows us to circumvent the necessity of finding p, inO <x < X. 
Before proceeding with this, however, we note that the integral terms in the other 
equations are of higher order in .V/ than the dominating terms and a few simple estimates 
lead to the results 

Pi =p, 


ur, = “(1 + A), (4.5) 


i 


with X = (y — 1)¢/yrR' po , when q is defined by Eq. (3.13). The primes indicate the 
time-dependent components of the state variables and velocities. To obtain a more 
useful form for Eq. (4.1), we write the conservation equations for a heat conducting 
gas in a passage assuming a given heat input rate’* Q/rR* and a one-dimensional flow. 
For brevity, we accept the implication of Eq. (4.4) that p = p(¢) only, in the x range of 
interest (i.e. 0 < x < X), rather than obtain this result again from the momentum 
equation in differential form. The other conservation equations (mass and energy) are 


[with p = p(d)] 
(pu), + p, = O, t .(5) 


puh, + ph, — (k/c,)h,, = Q/aR’, (4.7) 


where hf is the enthalpy 
Writing p = po(x) + p’(a, t), u = U(x) + wu'(az, 0), etc. and noting that the steady- 
state pressure is a constant, p) , we have 


pou const = A, Polo = const = [y/(y — 1) |po 


and 


Aho.2 — (ke/€:)ho.2z = Qo/aR?. 


13) is the heat input per unit time per unit distance in the flow direction. 
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From this, ho(x) can readily be obtained when Q, and h(— ©) are known. The pertur- 
bation terms in the sum: [Eq. (4.7) plus h times (4.6)] lead to the equation 


| aoe (pout + pt) — (k/e,)ht. = Q’/aR’, (4.8) 
Lf 
where the primes denote the fluctuating components. 

It follows that 

’ p= tatk 2) —-1 , 7 | 
on nt + (eR?) [Qt dx — ~~ pi |. (4.9) 
YPo c Jo 7-4 

Thus, us — uj is given by Eq. (4.5) only when X is large enough so that h! vanishes to 
our order of accuracy. That is, kh!/c,po < u{ . There are, of course, two contributions 
to h! in the downstream region; that associated with the acoustic wave and that of the 
entropy wave. If the contribution of each wave to each side of this inequality is esti- 
mated, it is seen that only when the entropy wave has decayed essentially to extinction 
(this occurs in a very small fraction of an acoustic wave wave length) does the inequality 
hold. Thus, we must choose XY downstream of the heater at such a distance that only 
the ‘isentropic’? waves are still of appreciable magnitude and then Eqs. (4.6) and (4.7) 
constitute the boundary conditions on the upstream and downstream acoustic waves. 
Having determined the acoustic waves, one could, if it were of interest, use Eqs. (4.6) 
and (4.9) to find w’ just at the rear of the heater and use this value together with u; 
to find the contribution of the entropy wave ¢“’. With ¢“’ so determined, the decaying 
density, temperature, and velocity, fields behind the heater are known. 

Equations (4.4) and (4.5), applied to the contributions of g“” and g”’ at x = ¢ 
(to obtain algebraic simplicity without essential loss of accuracy), lead to the charac- 


p 


teristic equation 


sin [ki(1 — o@)] — iz, exp [—7ki(1 — o)], —ylsin ki — iz. exp (—ikio)] 


| = 0 
cos [ki — ¢ z, exp [—7ki(1 — o)], (1 + A)[eos kic — 2, exp (—ik3e)] (4.10) 


where 


1 — | N, | exp (—kj’o) 


II 


Z - V, | exp [—ky’(1 — o)], 22 


This is a rather messy transcendental equation for a but it can be treated by noting 
that \ and 2, are small compared to unity (but larger than the contributions we have 


consistently omitted). For the case \ = z; = 0, the characteristic equation contains 
only real contributions (for real a) and has a real eigenvalue, a. If k, = k, , it reduces 
to sin k, = 0 and the solution of major interest is k, = 7. It is easy to compute (numeri- 


cally) the eigenvalues for a,/a, = k,/k, ¥ 1, and such a solution can be used as the basis 
of a perturbation calculation for a. Denoting the zero order values of k{ and k; by ko 
and ks») we let k’ ko + k;, . The characteristic equation can now be linearized in the 
k;, noting that k,,/k, a,/a, = v, and the resulting equation has in general, a complex 


solution. 
If we denote by yu, , the imaginary part of k,, , we obtain 


(1 — o) + v’¢)] cos [koopa] cos [kyo(1 — o)] — vsin [kee] sin [kyo(1 — o) Ju 
gm (X) sin [kio(1 — o)] cos [koe] — (22 + 2,) sin [kyo(1 — o)] sin [hoe] 
— (z, + 2) cos [k,o(1 — o)] cos [koe], 
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and a negative value of u, implies a negative imaginary contribution to a. Thus, the 
tube will sing when yp, < 0. It is consistent, of course, that gm (A) is always negative 
and the z, are positive functions of real a. u, will be less than zero for the lowest mode 
only when the bracket multiplying it is negative. This corresponds to having the heater 
in a prescribed lower portion of the tube which would be ¢ < 1/2 for a, = a, but which 
is ¢ < K where K < 1/2 for a, > a, . Together with this position requirement, a 
negative u, also requires that the \ contribution to the right side of the equation be 
larger than the other terms. It is not a desirable task to tabulate the dependence of 
this criterion on the many parameters involved, but the relation of certain predictions 
to the observed facts can readily be discussed. For a velocity, wv) , of 1.7 ft/sec and a 
heater temperature of the order of 800°F, the prediction of the foregoing analysis is 
that Jm (A) is large enough so that the tube should sing for ribbon breadths of .08 em 
or greater, but that a wire of less than this size should require a higher temperature. 
These facts were verified by experiments conducted by J. J. Bailey at the Harvard 
University Combustion Laboratory. He observed that ribbons of breadths 3/16’, 3/32”, 
1/4’, etc., play well at such a temperature (in an apparatus with L = 75 cm, 1 = 90 cm, 
R = 5 cm), but that a 1/16” wire needed such a large temperature that several wires 
were burned out during the experiments. 

Many investigators (see the bibliography) have noted that the tube will sing at the 
fundamental only when o < L’/2 (L’ denotes a correction for temperature). Bailey, 
however, has successfully ‘‘played” the second harmonic by placing the ribbon in the 
range 1/2 < o < 3/4. The eigenvalues of the characteristic equation predict a growing 
wave at essentially the second harmonic frequency under these circumstances. 

Another experimental fact not reported by other investigators is the following. For 
a given heater position and various ribbon temperatures (and, hence, various sound 
levels), the ratio of intensity of second and first harmonics has been recorded. This 
ratio increases (monotonically and in an experimentally repeatable manner) with increase 
in sound level. This is consistent with the qualitative discussion of the large amplitude 
heater response at the end of Sec. 3. There, it was noted that an increasing fraction of 
the heat release induced by the fluctuating velocity could be expected to appear at 
frequency 2w as w/u, increased. This is indeed the case and may well be the major 
influence in deciding the sound level of the tube. The conjecture is this: as the sound 
level increases, more and more of the energy associated with 9m X is stolen by the second 
harmonic until the first harmonic contribution of 9m (A) is balanced by the z; . This is a 
crude picture which could readily be made more precise by a detailed knowledge of the 
large w/w. heater response. 

We have not presented a quantitative account of the Harvard experiments here 
since the instrumentation was carried only to an accuracy consistent with the foregoing 
description. However, this singing tube analysis would not have evolved to the present 
stage without the enthusiastic experimental accompaniment of J. J. Bailey. 
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—NOTES— 


THE POST BUCKLING BEHAVIOR OF A CLAMPED CIRCULAR PLATE* 
By SOL R. BODNER (Polytechnic Institute of Brooklyn) 


1. Introduction. The von Kaérman equations for large deflections of elastic plates 
have been completely solved by K. O. Friedrichs and J. J. Stoker for the case of a circular 
plate simply supported around the circumference where it is subjected to a uniform 
radial thrust in the plane of the plate, Refs. [1] and [2]. Their solution describes the plate 
behavior from the initial buckled state to the condition of the ratio of the edge thrust 
to initial buckling load becoming infinitely large. 

The purpose of this note is to apply the methods developed in [1] to determine the 
post buckling behavior of a clamped circular plate. A modification of one of the methods 
is made which reduces the computation necessary to analyze certain aspects of the plate 
behavior. 

2. Mathematical formulation. For assumed radially symmetric deformations, the 
von Karman equations for a circular plate are a pair of non-linear ordinary differential 
equations, each of second order. These equations are 


Vo = w,,(1l/nw,, , (1) 
(yh)? Vw + (1/n¢.,w... + (/ne.,-w., = 0, (2) 


where ¢ is the stress function of the membrane stresses, w is the deflection of the middle 
surface of the plate, r is the radial coordinate from the plate center, A is the plate thick- 
ness and y” = 1/12(1 — »’) where » is Poisson’s ratio. The radius of the plate is R and 
the quantities obtained by differentiating @ are stresses per modulus of elasticity, F. 
Equations (1) and (2) can be simplified by introducing new dependent variables p and 
q where 


p = (1/r)¢,, , the compressive radial membrane (3) 


stress divided by Young’s modulus, 


q = —(R/r)w,, . (4) 
Equations (1) and (2) then become 
{r*[Gp — (1/2)q°]}., = 9, (5) 
{r°[n°Gq + paq]}.. = 9, (6) 
where 
Gq = (r°q..),-(R’/r’) (7) 
and 
n= Vh?/R’. (8) 


*Received Dec. 9, 1953. 
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The integration of (5) and (6) vields 


Gp — (1 2)q° sx (). (9) 


nGq+pq = 0. (10) 


The constants arising from the integration are zero by continuity considerations. Equa- 
tions (9) and (10) with appropriate boundary conditions completely define the problem. 
The boundary conditions at the edge of a clamped plate subjected to a radial edge 
thrust p* are 

p = p* at r= Rk, (11) 

q=0 at r= fh. (12) 
The two other boundary conditions needed for a complete solution are at the center 
of the plate and are due to symmetry requirements. They are 


Pp, = 0 at r= 0, (13) 


= 0 at r= Q. (14) 


3. Methods of solution. The boundary value problem considered here depends essen- 
tially upon one parameter N: the ratio of the edge pressure p* to the lowest critical 
pressure p,, . In [1] and [2] a perturbation method is used to solve (9) and (10) for low 
values of N, 1 < N < 2.5. This perturbation method consists of expanding p, q and 
p* as power series in a quantity e which vanishes at the onset of buckling, i.e. at NV = 1. 
The value of p., is obtained as one of the steps in the perturbation method and was 
found to be 

Der = 14.687. 
This method becomes increasingly tedious as N increases and it is convenient to use a 
power series method for N > 2.5. 

The power series method is mathematically simpler than the perturbation method 
but necessitates an estimate of the membrane stress and the curvature at the center of 
the plate. This estimate can be obtained from the results of the application of the per- 
turbation method in the region of lower NV. The power series method can be applied, 
however, without previous use of the perturbation method by utilizing an important 
aspect of the plate behavior. 

One of the notable results of [1] and [2] was that at large values of N the membrane 
stresses become tensions in the interior of the simply supported plate and change abruptly 
to compressions in a narrow “boundary layer’ at the plate edge. Since the same general 
behavior can be expected to occur for a clamped plate, the power series method can be 
modified to obtain the value of N at which the radial membrane stress at the center 
of the plate, p, , changes from compression to tension. The solution of this problem can 
serve to give the estimate of the physical quantities needed for the further direct applica- 
tion of the power series method. 

An asymptotic solution is used to determine the limit situation as N becomes in- 
finitely large. From this solution the limiting value of the membrane stress in the interior 
of the plate can be determined. 

4. Power series method. The introduction of a new independent variable a and 
new dependent variables z and x into (9) and (10) permits a solution of those equations 


by power series expansions. 
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Hence, if 


a = Ar/R, 0O<a<A, 
(15) 
x = p/A’n, k= q/A*n, 
where A is an arbitrary parameter, (9) and (10) become 
(1/a*)(a*m,«).q = (1/2)x’, (16) 
(J ‘e’)(a Ka) 0 +nrx«= 0, (17) 
with boundary conditions 
T 1* p*/n A®* at a= A, (18) rr. =0 at a =Q, (20) 
xg = @ at a= A, (19) ka =0 at a= 0. (21) 


The variables 7 and «x can be represented in the form 
2k » 
r= 7 Wp ,y (22) 
hk 0 


c= } Ka, (23) 

k=0 
The first term in each of these expansions is a function of physical quantities at the 
center of the plate. The remaining coefficients of each series can be expressed as functions 
of the first terms, 7) and x , by means of the recursion relations obtained from sub- 
stituting (22) and (23) into the differential equations (16) and (17). These recursion 


relations are 


2h(2k + 2)r, = (1/2) > KmKn 5 (24) 
m+ne=k-1 
Q2h2k + 2)mq=— Yo mak, . (25) 
m+n=k—1 


The power series method of [1] and [2] consists of prescribing 7) and «x» and solving 
for A and x* from the equations obtained from boundary conditions (19) and (18). 
The corresponding values of p, qg, p* and N can then be calculated from the previous 
formulas. 

The value of A is the lowest root of the equation 


> «A™ = 0, (26) 


k=0 


which is boundary condition (19) at the plate edge. This equation can be solved for A 
with a minimum of labor if A ~ 1. For A = 1, my = po/n’ and ko = Qo/n where py and 
do are the values of p and qg at the center of the plate. These values can be estimated for 
a particular N from the results obtained by the application of the perturbation method 
for lower values of N. 

From the results of [1] and [2] it is reasonable to expect that, for the case of a clamped 
plate, an N exists for which p, = 0. Therefore, for this special case, 7 = 0 and all the 
coefficients in the expansions (22) and (23) can be written in terms of one parameter, ko . 
The value of x) can be determined from the solution of (26) in which A can be prescribed 
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to be unity. That is, xo is the lowest root of 


> x, = 0, (27) 


which was found to be x) = 27.11 using coefficients up to x.) . The corresponding value 
of x* was found from (22) and (18) to be 29.11. The value of N for the condition po = 0 
was found from the relations N = p*/p,, = 2*/14.68 to be 1.98. 


With this information the first section of the curve for the variation of the radial 
membrane stress at the plate center, po , with N can be sketched as shown in Fig. 1. 
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A similar curve could be drawn for gp . The general behavior of these curves, from physical 
considerations, should not differ from those obtained in [2] for simply supported plates. 
These curves can then serve to provide estimates of po) and q for use in the further 
application of the power series method. When the limiting case of N — @ is solved, an 
estimate of p, and q can be obtained for the complete range of N, (Fig. 1). In this 
manner the power series method can be conveniently used for the range 1 < N < 135, 
thereby being independent of the results of the perturbation method. For N about 
equal to and greater than 15, the power series method becomes unwieldly and an asymp- 
totic development is preferable. 

5. Asymptotic solution. The solution of (9) and (10) for N tending to infinity for 
the clamped plate follows, in general, the method described in [1] for the simply supported 
plate. The computation necessary for the clamped plate was greater than that required 
for the simply supported plate. This solution shows that as N — © the interior of the 
plate acts as a membrane subjected to a uniform tensile stress equal to .131 times the 
compressive thrust at the plate edge. That is, 

lim p = —.131p*. 


V 


This is also the asymptotic value of p, and is indicated on the graph, (Fig. 1). 
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AN APPROXIMATE SOLUTION TO THE NAVIER-STOKES EQUATIONS* 
3y MORTON MITCHNER (Harvard University) 


The purpose of this note is to show how a new approximate solution of the Navier- 
Stokes equations may be constructed from any given exact solution having a certain 
specified form. We shall suppose that we are given an exact solution of the Navier-Stokes 


equations for an incompressible viscous fluid having a velocity field q’ = (q{ , 93, q5) = 
(u’, v’, w’) specified in the form 
yu’ = all(y, s. 
yp’ = aVi(y, 2, 0, (1) 
w’ = aW(y, z, @). 
U, V, and W denote three functions of the position vector r = (x, , ®2 , 23) = (a, y, 2), 


and the time coordinate ¢; a denotes a dimensionless constant. For consistency with the 
equations of motion (upon taking the divergence of the Navier-Stokes equation, and 
employing the continuity condition), the pressure p’ (and density p) must satisfy 


| is 0q' dq) 
—! vip = 5 Sai oat 
p ik OL, OX; 
Hence, it is sufficient to assume that the pressure field has the form 
p’ = a Ply, z, t). (2) 
Assuming the existence of the above exact solution, we can construct a new approxi- 
mate solution [q = (u, v, w), p] of the equations of motion for an incompressible viscous 
fluid, and this solution is given by 


u = ula + by) + u’ — (t — by) du’, 


, 


v=v', 
(3) 

w= w’, 

p = constant, 





*Received Dec. 26, 1952. Revised manuscript received Jan. 7, 1954. 








402 NOTES [Vol. XII, No. 4 


(uw , a, b, and ¢, are constants). This solution is valid to first order in a for ¢ = é under 


the hypotheses that 


(a) the spatial derivatives of q’ are bounded for t = t, | 
(b) a <1, (4) 
(c) a buy (t — th) K 1. 


The velocity field q(y, z, t) may be regarded as the subsequent time development of an 

initial disturbance q’(y, 2, é) superimposed upon a Couette-type shear flow at time ¢, . 
The validity of the preceding statement may be checked by direct substitution into 

the Navier-Stokes equations. Thus, for the u component 

Ou 1 Op 


Ou Ou 2 
+ u + y ay + u _e “= ra + vV~u 


Ou 


ot 
and hence, employing (3), 


ou’ ou’ Ou’ . | ov’ 2 
— +’ — + w’ — — ult — &)| — +0 — + wv’ — 
ot Oy 7 Oz 0) Ot + oy ” Oz 


= vpVu’ — but — t)vV’. (5) 


In virtue of the fact that wu’ and y’ are exact solutions of the Navier-Stokes equations, 
Eq. (5) states that 


+. 9. (6) 


bu,(t — ty) 
p oy 


But p’ = a’P(y, z, t) and for t = t, , P/dy is bounded. The term on the left side of 
Eq. (6) will therefore be of second order in @ provided ¢ also satisfies the condition 


a’ bu,(t — ty) Ka, or abut — t)) «1. 


In a similar fashion, it may also be shown by direct substitution that v and w satisfy 
the Navier-Stokes equations. 

A particular example of the above general result is provided by the known velocity 
field describing the decay and diffusion of an infinitely long vortex filament initially 


concentrated on the z axis. 
ue = 0, 


_«,u- exp (—r* /4vt)] 
Qa ~ r ‘ (7) 


K_ [1 — exp (—7’/4r)] : : : 
w= a —— — r=Yy . 

or Y Yr ’ y + 
K denotes the initial circulation or strength of the vortex filament, and is to be associated 
in its dimensionless form, Kb/uo , with the parameter of smallness a. For any t, ¥ 0, 
the solution (7) satisfies the hypotheses (4), and consequently Eqs. (3) and (7) describe 
the behavior [for ¢ — t, 2 0, but (¢ — &) <K 1/buoa] of a vortex filament superimposed 
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on a Couette-type shear flow, the vortex filament being aligned with the direction of 
flow. 

The particular solution stated above exhibits some interesting properties as regards 
the exchange of vorticity between components. We note that 


, _ Ow av _ dw’ dv’ i 
Oy 0z Oy Oz 7 
Ou ow Ou ov’ 
= — =—= —(t — 4)bu — 
” Oz Ox Oz 0) Duto dz’ 
Ov Ou _ Ou _ 


Dn ee ile eh ee ee 

Ox Oy OY oy 

The vorticity associated with the vortex filament disturbance, ¢’, remains completely 
unaffected by the presence of the shear flow and proceeds to decay as if the shear flow 
were absent. However, the 7 component of the vorticity, initially zero, begins to grow 
at the expense of the vorticity of the shear flow. The action of the component is thus 
analogous to that of a chemical catalyst; while é itself remains unaffected by the shear 
flow, it causes a production of y, drawing upon the infinite field of oriented vorticity in 
the shear flow. It may be shown that this phenomenon (that ~ is unaffected by the 
presence of the shear flow) is actually independent of the specific form of the shear flow. 


APPENDI X* 

In connection with the preceding remarks, the editors of the Quarterly have brought 
to the author’s attention an investigation by Berker.** Although there is no direct con- 
nection between the present work and that of Berker, there does exist a superficial 
similarity which may be worthy of clarification. 

Berker assumes that an exact solution of the Navier-Stokes equations for an in- 
compressible viscous fluid with respect to an inertial frame of reference oryz is pro- 
vided by the velocity field q’(x, y, z, 0). Using this given vector point function, Berker 
then defines a vector field qr(X, Y, Z, t) = q’(X, Y, Z, t) with respect to a moving frame 
of reference OX YZ. Corresponding to the motion qr(X, Y, Z, ¢) with respect to OX YZ, 
there will exist a motion q,(z, y, z, t) with respect to oryz. Berker then determines 
the conditions under which q,(z, y, 2, t) will be an exact solution of the Navier-Stokes 
equations for an incompressible viscous fluid (his Eq. 6.14) and indicates the construction 
of qz(2, y, 2, t) in terms of q’(z, y, z, t) (his Eq. 6.13). 

For the particular form of the assumed initial exact solution provided by Eq. (1), 

3erker’s new exact solution has the form 


uz(y,z,t) = a(t + w'(Y, Z, 0d, 
vay, 2, t) = b() — (2 —A)2+0v'(Y, Z, td cos Qt — w'(Y, Z, d sin 2, (8) 


w,(y, 2, t) = c(t) + (y — 124+ v'(Y, Z, sin Qt + w’(Y, Z, Bd cos AQ, 
*Received Jan. 7, 1954. 


**\ RR. Berker, Sur quelques cas d’intégration des équations du mouvement d’un fluide visqueux in- 
compressible, Institut de Mécanique des Fluides de L’Université de Lille, 1936. 
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where 2 is a constant, where a(t), b(t), c(t) are arbitrary functions of t, and where 
Y = (y — b) cos Qt + (2 — ©) sin 1, 


Z = —(y — b) sin Qt + ( — ©) cos NE. 


Comparison of (3) with (8) makes quite evident that these two solutions are essentially 


different. Whereas (8) is exact, (3) is approximate. Furthermore, (3) cannot be derived 


as an approximation from (8). 


A MEASURE OF THE AREA OF A HOMOGENEOUS RANDOM 
SURFACE IN SPACE* 


By STANLEY CORRSIN (Aeronautics Department, The Johns Hopkins University 


Introduction. We are given an indefinitely large space containing random surface 
or surfaces homogeneously located in the mean. The problem is to relate the average 
area of surface per unit volume of space to a simpler statistical quantity, in particular 
the average number of cuts per unit length made by a straight randomly directed 
sampling line with the surface. 

The plane case will be studied first. After the three dimensional case, illustrative 
application will be made to the problem of extending to two and three dimensional 
variables a theorem of 8. O. Rice on the average rate of occurrence of any particular 
value of a one dimensional random variable. Possible use in describing fluid mixing is 


also indicated. 


Two dimensions. Given a plane “homogeneously” inscribed with contour or con- 
tours of arbitrary shape. The homogeneity is statistical, i.e. any statistical function 
associated with the contours is invariant to a translation of coordinate system in the plane. 
Let £ be the average contour length enclosed in unit area and let n be the average number 
of cuts per unit length made by an arbitrary straight traverse line crossing the plane. 
For a non-isotropic field n is averaged over all traverse directions with uniform weighting; 
for an isotropic field, any single line will do. 

Draw a “very large’ square in the plane, L on a side, and subdivide it into “very 
narrow” traverse strips parallel to one pair of sides. 
averages over L or L’ are satisfactorily 


“Very large’ here denotes L so large that 
it requires that each traverse strip cross 


close to their asymptotic values. For example 
the contours a very large number of times and that the length of contour in L’ divided by 
L’ be as close as we like to £. “Very narrow” denotes 6 so small that virtually all of the 
intercepted contour segments in a strip can be approximated by secants. This gives 
restriction on the permissible number of corners and contour intersections. 

We imagine the square and strip structure rotated through 180° for averaging 
purposes in case the field is not isotropic. Then the average number of crossings in one 


*Received Jan. 18, 1954. Supported by the Office of Naval Research. 
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strip is nL. The number of strips is L/é. Hence the expected total contour length in 
the square is 
nL* 
£.. = * 0p, (1) 
where (l) is the average contour segment per crossing, so (1) can be written 


n 


£= (i). (2) 





To obtain (1), consider a typical crossing (Fig. 1). For any one crossing 
6 
| = ———_.- (3) 
| Sin 6 | 
However, for each traverse orientation, the probability of intersection at angle 6 


with a contour element is proportional to the projected length of the element on a line 
perpendicular to the traverse direction, i.e. to | sin @ |. With normalization this becomes 


the probability density of 


: sin 6 —“<e<* 
B.(0) = >" nd 2 (4) 
\O elsewhere. 
But 
=f wiya=2]  10B,(6) de. (5) 
Substituting (4) into (5) we find 
T 
() = 5 6, (6) 
whence 
2 - 
n=-&. (7)* 
T 


This result is easily seen to be consistent with the solution of Buffon’s “needle problem”’t 


*Arrived at independently by P. V. Danckwerts of Cambridge University (private communication). 
tSee, for example, Uspensky: Introduction to mathematical probability, McGraw-Hill Book Co., 1937. 
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in the special case when the contours are equidistant parallel straight lines. 

For an isotropic field £ can thus be determined by the intersection rate along a 
single sampling traverse. 

Three dimensions. Now we have a three dimensional space homogeneously inscribed 
with surface or surfaces of arbitrary shape. Define @ as the average surface area per 
unit volume of space and n as the average number of cuts per unit length made by a 
random straight line traverse across the space. For a nonisotropic field n is averaged 
over all directions with equal weight; for an isotropic field a single traverse suffices. 

Take a ‘‘very large” cube (Z on a side) and subdivide into “very narrow’’ square 
traverse tubes (6 X 6 X L) with all faces parallel to those of the cube. 

The average number of surface crossings in a single tube is nL. The number of tubes 
in L’ is L’?/5°’. Hence the expected surface area in the cube is 


L? 
G,. = nL — (a), (8) 
0 
where (a) is the average area segment per crossing. @;; = @-L*, so (8) can be written 
7 . 
@ = 3 (@). (9) 
5 


To obtain (a), consider a typical crossing (Fig. 2), with g the angle between tube 
axis and normal to surface element. 











i a m aca heganaiaantrtii as 
ae 8 
ath Ss “_ 
= r ; 
—— ar es a =— 
. * 
YY ah a ea ee 
~\ \ 8 
te _a 
ee —_——<— ppm ahi as or “s 
normal to surface element 
Fic. 2 
5 
a= —- (10) 
COS ¢g 


For each traverse orientation (i.e. cube orientation) the probability of intersection 
at angle ¢ with a surface element is proportional to the projected area of the surface on 
a plane perpendicular to the tube axis, i.e. to | cos g | times a measure of relative solid 
angle giving ¢, i.e. | sin g |. With normalization constant, this gives the probability 
density of ¢ 





F 
- T 
12 | cos¢ sin ¢ | os ¢S- 
B.(e) = 4 . (11) 
0 elsewhere; 
therefore 
i= | aB,(a) da = | a(y)B,(¢) de. (12) 


/0 “0 
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Substituting (10) and (11) into (12), we find 

(a) = 28°. (13) 
Whence 
(14) 


n= 


19 


Again, for an isotropic field, @ can be determined by the intersection rate along a single 
sampling traverse. 


Application I: Rice’s theorem in two and three dimensions. 
Given a two dimensional stationary (= homogeneous) isotropic random variable 
u(x, y). The iso-value lines corresponding to u = u, will be an isotropic field of contours 
in the plane, and we can seek the average occurrence length per unit area of the value u, : 


£.. =n. , (15) 


where n,, is the average number of u, occurrences per unit length along the linear samp- 
ling traverse in any direction. The traverse yields a one dimensional randam variable 
to which we can apply a theorem of Rice:* 


n.. = | | yu. ,u’) | u’ | du’, (16) 


d = 


where u’ = du/ds, the derivative in the direction of our sampling traverse and y(u, w’) 
is the joint probability density of u and its slope in this sample. From (15) and (16) 


£., = 5 yu. , u’) | u’ | du’. (17) 


In the special case of jointly Gaussian u and wu’, we use Rice’s simplified result to 
obtain 
l ro! —4?/2¥(0) 
a = “a i se l 
~ 7 vos ° (18) 


where V(c) is the auto-correlation, (u(s)u(s + @)). 
For three dimensions, the corresponding equations are 


« 


@,. = 2 | yu. , u’) | wu’ | du’ (19) 
and, for jointly Gaussian u and w’, 
; l eo} —u*/2¥(0) 
&. = =-\~-Sa te . 20 
«eS { vO) (20) 


For non-isotropic, homogeneous u(x, y), y(u, u’) is a function of traverse direction, 
and we must include the operation of averaging over all directions. 


*S. O. Rice, Mathematical analysis of random noise, Bell System Tech. J., 23 (3) and 24 (1), July 


1944 and January 1945. 
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Application II: Fluid mixing. 

The average interfacial area per unit volume is a significant measure of the “degree 
of mixedness”’ in an isotropic field of two molecularly immiscible liquids. It is, however, 
inaccessible to straightforward experimental determination. Equation (14) permits 
its calculation from the simpler process of interfacial-intersection counting along a 
linear traverse through the mixture. 

In the homogeneous mixing of two gases or molecularly miscible liquids (e.g. tur- 
bulent mixing) it is possible that the notion of interfacial area can be replaced either 
by the surface area on which the concentration fluctuation is zero, or by the surface on 
which the concentration gradient magnitude has a local maximum. 


EQUILIBRIUM OF MEMBRANES ELASTICALLY SUPPORTED AT THE EDGES* 


By V. G. HART (Dublin Institute for Advanced Studies 


Abstract. The problem considered is that of finding the statical deflection of a 
stretched membrane, subjected to a uniform pressure on one side and elastically supported 
at the edges. The deflection of the membrane is supposed small and the problem reduces 
to solving Laplace’s equation with mixed boundary conditions. Solutions are given for 
the cases where the bounding curve of the membrane is (i) an equilateral triangle, and 
(ii) a rectangle. 

1. Introduction. We consider the problem of finding the statical deflection of a 
membrane originally lying in a plane (the neutral plane), when subjected to a uniform 
pressure on one side, its edge being elastically supported. This means that the edge can 
move in a direction perpendicular to the neutral plane, but is restrained at any point 
by a force proportional to the deflection at that point. Small deflections only being 
considered, the tension is a constant, and the problem reduces to solving the boundary- 
value problem (3) for an edge of arbitrary shape; solutions are given for (i) a membrane 
in the form of an equilateral triangle and (ii) a rectangular membrane. 

Imagine a membrane stretched to a uniform tension 7’ and bounded by any plane 
curve B. A uniform pressure P now acts on one side of the membrane which takes up a 
statical deflected position with deflection w. At any point on the edge the tension gives 
a component of force per unit length in the direction perpendicular to the neutral plane 
of amount — 7 dw/dn, (dn being the outward normal element to B in the neutral plane), 
and this is balanced by the elastic force of constraint, which we write w/k, where k is a 
constant. 

The appropriate partial differential equation and boundary condition for w are 
therefore 


Aw —-P/Tin S, (w)p = —kT(Ow/dn), , (1) 


S being the domain of the membrane. 
We now make the transformation 


») 


*Received Feb. 1, 1954. 
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where r is the distance from any fixed point in the plane of the problem, and from (1) 
we find that u must satisfy Laplace’s equation with mixed boundary conditions, viz., 


‘ ° a e 
Au = Oin S, (. +c ) = : r+e . (: ’), (3) 
On/ p 2 on \2 ° 


where c = kT. Our procedure is to solve (3) for u in the case of the particular boundary 
B chosen, and then to find the deflection w from (2). 

We note that on putting ¢ = 0 in (3) the boundary-value problem reduces to the 
torsion problem, and this fact affords a check on the results obtained. 

The Dirichlet. integral of w, which is intimately connected with Laplace’s equation, 
represents in this problem the excess potential energy stored in the deflected membrane. 
The total energy V, stored in the membrane is 


we ey (22) ga 
V,==7 I | (2 +(5,) [as +e. (4) 


where C is the energy stored in the undeflected membrane. On taking up its deflected 
position, work is also done on the elastic support and the potential energy V, contained 
in this is 
; se © 
Vv, = w db, (5) 
B 


—_ > 


where db is an element of length of the boundary B. Using the boundary-condition 


on win (1) and also Green’s theorem we find 


rs | ‘ se | dw\’ dw\* ; 
3! | ised I, (22) r (=) | - @) 


: P| wdS—V,+6, 


- . S 


S. Thus the total potential energy V stored in membrane and 


since Aw —P/T i 


support is 


as is indeed obvious. 
If B is a circle, the solution of (3) is very simple, viz., 


u = 3a’ + a, (8) 


where a is the radius of the circle, and the origin is taken at its centre. The corresponding 


Pa’ r 7 ' 
= == > y ° 9) 
’ 17 [ ay +2 | \ 


deflection is 
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2. Membrane in the form of an equilateral triangle. We now consider the case 
where the membrane is an equilateral triangle of height 3a. We take the origin and axes 
as in Fig. 1. To solve the boundary-value problem (3) for this domain, we seek a har- 


y 
4 





—» X 





ge =e (] —p> 


a. oe 








Fic. 1 


monic function possessing the symmetry properties of the triangle and also satisfying 
the boundary conditions. It is obvious that a function with one of the symmetries of 
the triangle—invariance under a rotation of axes through 27/3—will satisfy the boundary 
conditions on all three sides of the triangle if it satisfies them on one—say on sr = —a. 
Bearing in mind the other symmetry of the triangle (reflection in Ox), we accordingly 
choose 





u=A+ Be + 2*°), z=a2+1y, 2* =x — iy, (10) 
which is harmonic and has both the symmetries. We adjust the constants A and B to 
satisfy the boundary condition on x = —a, which reads: 

Ou Bove 2 
u—-c=— =-(a +y) + ca. (11) 
4 ne 2 


On substituting from (10) we find A and B by comparing coefficients of y 


a(2a" + 6ac+3¢) 7, 1 (12) 
3(a + c) 12(a +c) = 
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The problem (3) is now solved for u, and by (2) the deflection 
Pa? [3c?/a? + 6c/a +2, 2° — 32y’ es ] 
w= La . 1; 
; 2T | 3(1 + c/a) 6a°(1 + c/a) 2a @+y) (13) 


The level lines are approximately circles for large values of the dimensionless parameter 
c/a, i.e. for a very weak elastic support on the edge. 
The total potential energy stored in the membrane and support is by (7) 


3/2p2.4 2/2 ee, 
oS 2 6 | 1o¢/a + ele +8) (14) 


v= C+ or 1+c/a 








3. Rectangular membrane. We consider now the case of a rectangular membrane 
(length 2a, breadth 2b). The axes are taken so that the sides have the equations x = a, 
y = +b. 

To solve the boundary-value problem (3), we introduce an auxiliary function v by 
the transformation 


u = (2? — yy’) + BW? + Wb +2, (15) 
and (3) becomes the following boundary-value problem for v 
Ov 
(A) vtc—=0 on y = +b, 
*y 
Av=0 in S, (16) 
Ov 9 9 
* vt = =y —b —2cb on x = +4; 
L d 
note the homogeneous boundary conditions on y = +b. 


We propose to solve (16) for v and then find the deflection w through the trans- 
formations (15) and (2). 
Since v is harmonic we take as the solution of (16) 


v= > A, cosh (A,2) cos (A,4), (17) 


where A, and A, are to be determined. Application of the boundary condition (16A) 
gives the following equation for X, 


cot (br,,) = cA, « (18) 
This equation determines an infinite sequence of values. We shall use only the positive 


values, numbering them in order of increasing magnitude. For large values of n we have 


NE (19) 


cnr ~ 
The functions cos \,y form an orthogonal set; direct calculation gives by use of (18): 


ab 0 for m # Nn, 
| cos (A,,7) cos (A,y) dy = (20) 
oe lb +csin’ (r,b) for m=n. 
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We expand the right-hand side of the boundary condition (16B) in terms of these ortho- 


gonal functions, writing 


y — bb — 2be = } B,, cos X,Y, (—b<y< b); (21) 


the coefficients are found by using (20) and are 
B, = —4sin (A,b)/[AR(b + ¢ sin’ \,))]. (22 
The boundary condition (16B) is now applied to v as defined by (17) and this requires 


that A, should satisfy 


y’ — b° — 2b = > A,(cosh X,a + cd, sinh X,a) cos X,Y, (—b<y <b). 23) 


Comparing (23) with (21) we see that 

A, = B,/(cosh X,a + eX, sinh X,a), (2-4) 
where B, is defined by (22). The coefficients A, in (17) are thus determined and so the 
solution v of (16) is found. 
Tsing (15) and (2) we find the deflection to be 


Pl? 2 y\? 
w= > E +2 — (¥) 
vi | b bh 


= 


— 


: sin” X,,b cos (X,Y) cosh (A,r) 
<; (A,b)*[1 + (sin 2d,,b/2X,)](tan X,.6 + tanh d,a) cos (A,,b) cosh (A,a) _}’ 


where we have simplified the coefficients A, in (24) by the use of (18). 

As already remarked we get the torsion problem on putting c = 0; Eq. (25) is then 
easily seen to reduce to the known solution for the torsion of a beam of rectangular 
section, A, now being (2n + 1)x/2b. 

Due to the presence of the term cosh (A,x)/cosh (A,a) which has the asymptotic 
— a)] when n is large [A, being given by (19)], the series in (25) is 


e)<a@,iwi <b. 


value exp [A,(| x 
uniformly and absolutely convergent in both z and y for the intervals 
Also if we differentiate twice with respect to z or y, the resulting series is still uniformly 
and absolutely convergent for the same reason. Thus the operation of term-by-term 


differentiation is legitimate, and so w, as in (25), is the solution of the boundary-value 
problem (1) for a rectangular boundary. 
The potential energy stored in the membrane and support is 
; 2P°b* | a ( a, 
V == |7(14+ 35) 
37 D A) 
 ® sin v0 __tan (An) tanh (A,q) 4 ¢ 26) 
r,,0)°|1 + (sin 2X,,6 2X,,5) | (tan X,.b + tanh A,a) 


| 
This statical problem for the rectangular membrane as discussed above is mathe- 
matically similar to a dynamical problem treated by Rayleigh (Theory of sound, Vol. I, 
pp. 200-204) viz., that of a vibrating string with ends elastically supported; this involves 
the equation (18) for the computation of proper frequencies in a certain limiting case. 
In conclusion, I would like to express my gratitude to Professor J. L. Synge for his 


aid and criticism in the preparation of this paper. 
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THE SOLUTION OF EIGENVALUE PROBLEMS OF 
INTEGRAL EQUATIONS BY POWER SERIES* 


By J. R. M. RADOK (Dept. of Supply, Aeronautical Research Laboratories, Melbourne, Australia) 


Summary. The characteristic equations and eigensolutions of Fredholm integral 
equations of the second kind with symmetrical kernels of the homogeneous polynomial 
type are obtained by use of power series. The method is applied to the particular case 
of the kernel | x — y | and it is indicated that it may readily be extended to more general 
types of equations. 

1. Introduction. One of the standard methods of solution of ordinary linear differ- 
ential equations makes use of power series. Although certain types of integral equations 
correspond to such differential equations, to the author’s knowledge no equivalent 
method for the direct solution of these equations has been developed up to date. It is 
the object of this report to fill this gap. 

In [1] a method of solution of Fredholm integral equations of the second kind with 
symmetrical kernels of the polynomial type has been developed by Heller and Radok. 
This method was conceived as an iteration process. The only difference from the normal 
iteration procedure arose from the fact that the initial function, which in the commonly 
used process is arbitrary, was determined by the kernel of the integral equation. The 
choice of the initial function was thus based on the observation that the integral operator 
in the equation transformed an arbitrary power of the integration variable into certain 
“dominant” terms with exponents independent of the exponent of the power originally 
substituted and into one term of higher order. In other words, if the integral operator 


is given by 


Le = F K(x, yely) dy, (1.1) 
then 
ty") = x"**c,, + e eS 2 (1.2) 
ln 
where yu, v, k only depend on the kernel. Substituting for the initial function gy, a poly- 


nomial with arbitrary coefficients of the type represented by the “dominant” term in (1.2), 
applying to this initial function an arbitrary number of iterations and comparing then 
coefficients of powers of x on both sides of the integral equation 


¢ = AL(y), (1.3) 


the corresponding number of exact terms of the power expansions of the eigensolutions 
and of the characteristic equation may be obtained. 

In this report the above process, which still contains formally some of the disad- 
vantages of iteration processes, will be generalized into a straightforward power series 
method. In order to simplify the analysis, consideration will here be restricted to kernels 
of the homogeneous polynomial type. However it should be stressed that the method is 
easily extended to much wider classes of kernels. In fact, it will become obvious that 


*Received Feb. 17, 1954. 
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it is applicable to kernels involving a finite number of linearly independent functions, 
in which case the desired eigensolutions have to be expanded in infinite series of the 
complete system of functions occurring in the kernel. In particular, this will be the case 
for kernels involving trigonometric functions, when the eigensolutions will be obtained 
in the form of Fourier series. 

The general method of solution will be developed in Sec. 2, while one of the simplest 
kernels of the special type under consideration will be investigated in Sec. 3, in order to 
demonstrate the ease with which the method may be applied to actual problems. In 
Ref. [1], two more complicated types of kernels have been studied for which the exact 
solutions are known and which arise from the problems of natural vibrations of a uniform 
beam and of a wedge, clamped at one end. 

2. The general method. As indicated in the introduction, consideration will be 
restricted to Fredholm integral equations of the second kind with symmetrical kernels 
of the homogeneous polynomial type, i.e., to equations of the type: 


g(x) = AL(y) = al | P,_:(x, yoy) dy + | P,_i(y, x)oly) ay |, (2.1) 
where 
n—-1 
P,.(z,y) = D> Aaz'y’, i+tk=n-—-1. (2.2) 
(i,k) 


If it is assumed that g(x) is continuous and differentiable, it is easily seen that ¢ 
satisfies a differential equation of the form 


n—1 
ge” =r Dine”, (a) 
where p; are polynomials in x; hence the solutions g must be integral functions, i.e., 
\ - l , 
g(z) = az", (2.3) 


and the use of power series is justified. 
Consider the operator 


L(y) = [ P,_,(x, yol(y) dy + [ P,_i(y, x)oly) dy. 2.4) 
In particular, 
L(y’) = e,2"*' + > 2 As, (2.5) 
where 
endl eT 2 
aia p> (+i : Ey +1) (rz ) ) Pea a ae 
2 ? An _ of 1 1=0,1,-:: _ 
ma rees i es eee ee 8 r 
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Applying the operator (2.5) to (2.3), substituting in (2.1) and comparing coefficients 
of the powers of x, one finds 


Qna1 = ACA; , l>n (2.8) 


and hence, forO0 < k < n — 1, 


v~-l 
Apaivn = ’a, Il Crinu ’ v > ie (2.9) 


p=0 


<l<n-— 1 give 


n—1 @ 
a; _ A ze p> Y | ewe ee 


Finally, the powers of x for 0 


n-1 re) v—1 
mAD, 2», APNE Lh Gir s 1=0,1,-::,n—1. (2.10) 


In order that the system of equations (2.10) will have non-trivial solutions a, , one has 


5:, — AG,,(A) | = 0, lLu=0,1,---,n-1, (2.11) 
where 
co v-1 
G06 = 2 he ol Been « (2.12) 


The determinantal equation (2.11) determines the eigenvalues of the integral equation 
2.1), while the corresponding eigensolutions are given by 


(2 


n 1 oe] vy-1 
alt) = >> dS n’x"*"a, TT cine « (2.13) 
p=0 v=0 «=0 


It is easily verified from (2.6) that (2.13) converges with » for ail x, the rate of con- 
vergence depending on yu and d. The same is true for the function G,,,(A) of (2.12), 
i.e., for (2.10). 

3. Example. As an example of the application of the method of solution of Sec. 2 
the following integral equation will be investigated. It is one of the simplest of the type 
under consideration and in many text books, e.g. Biickner [2], is chosen for the demon- 
stration of methods of solution: 


ax 1 
g(x) = a | (a — y)e(y) dy + [ (y — x)e(y) in|; (3.1) 
hence, in the notation of Eq. (2.2), 
P,_i(z, 9) =z— y, i+tk=1, n= 2, (3.2) 
Ao = Rw = 0, Ass = —l, Bis => * 


The integral operator in the present case is 


z 1 
L(y) = [ (x — yoy) dy + [ (y — x)g(y) dy 
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and, in correspondence with (2.5), one has 


i x l 
L(y) = _ - ~ : (3.4) 
(1 + 1)(l + 2) 1+ 1 l+2 
since by (2.6) 
by Ay 
) l+2 ia l+1?’ 
) 
2 ve 
C= ° (3.5) 
(1+ 1)(1 + 2) 
Proceeding as indicated in Sec. 2, comparison of coefficients of powers of «x gives 
2da; : 
G12 = —— i> s (3.6) 
(l + 1)(1 + 2) . 
é y ‘ ¥ 
(2d) "avy _ (2d) a, ma 2 7) 
a !? Q+a = 1? y » fi, (5.4 
(2p)! (Qn + 1)! 
(2d)’ (2A)” 
a = Na +r a, - 
Xo, + 2)(2r)! Xo + 3)(2v + 1)! 
3 8a) 
— (2d) (2d)” 
a = Al Ao + a, > 
; La +1! 1 2 +d)! 
The system of Eqs. (3.8a) which corresponds to (2.10) is easily seen to lead to the deter- 
minantal equation 
° sinh ® 
cosh ®o—- WwW sinh 23) a cosh Ww 
Ww ‘ \ 
= Q (3.9a) 
w sinh w 1 + cosh w 
or, finally, to 
2 + 2 coshw — wsinhw = 0, (3.9b) 
where 
Q2=w 3.10) 


The elgensolutions are cviven by 


=~ (2X) "x . (2A) x2 ‘ 
p A» = eg 3 L 1)! 


(2y)! » (2v 4 


3.11) 


a 
+ sinh we. 


= dy cosh WU 
@ 


4. Conclusions. A method has been developed by which the eigenvalues and eigen- 
in types of Fredholm integral equations of the second kind may be 


¢ 


functions of cert 
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determined. The method which essentially corresponds to the method of power series, 
used for ordinary differential equations, may be applied to more general kernels than 
those of the homogeneous polynomial type considered here. It may also be extended to 
other systems of linearly independent functions, such as trigonometric functions, depend- 
ing on the construction of the kernel which in that case would have to involve a finite 
number of trigonometric terms. 

The method should prove of great value in the treatment of vibration problems of 
systems with polynomial mass distributions and concentrated masses. Two of the 
simplest problems of this type, whose exact solutions have been known for a large 
number of years, have been solved by a closely related method in [1]. 

5. Acknowledgment. The author is indebted to Prof. A. Pfluger at the Eidgenéssische 
Technische Hochschule, Ziirich, Switzerland, for some valuable suggestions which 
greatly assisted the further development of the method. 
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ASYMMETRICAL FINITE DIFFERENCE NETWORK FOR 
TENSOR CONDUCTIVITIES* 


By L. TASN Y-TSCHIASSNY (University of Sydney, Australia) 


In a paper recently published in this journal [1] the problem dealt with was described 
by its author in the following way: “Given a region in which Eq. (4) holds and a large 
number of points in the region chosen at random, in what way should the points be 
interconnected with ‘physically realizable’ electrical resistors in order that the voltages 
at the nodes shall be as nearly as possible the correct solutions of the boundary value 
problem characterized by Eq. (4) and appropriate boundary conditions?” 

The Eq. (4) quoted is the special case of the differential equation 


V-{ole]- VS} +7 =0 (1) 


for [e] being the unit tensor, i.e., for isotropically conducting material. In Eq. (1) ® is 
the electrical potential in a two-dimensional continuum, 7 the current density of a 
distributed transverse external source, o the conductivity of the continuum in a certain 
direction, and [e] the non-uniformity tensor of the conductivity. The scalars 7, o, and 
the tensor [e] may be functions of the position. 

The author of the present Note had been interested in a practically identical problem 
for which he had coined the expression “triangulation of a two-dimensional continuum” 
and which included the general case of tensor conductivities [2]. He replaced the interior 


*Received March 3, 1954. 
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of a triangle of any suitable shape and small size by a delta of resistors (see Fig. 1) 

so that, with the exception of the boundary, two resistors in parallel always connected 

a node pair in the final equivalent network. If the principal (maximum or minimum) 
fa 

y| “AN As 


) 








Az 


% 
A (4) 





i (a) 


Fig. 1. Replacement of the interior of a triangle by a delta of resistors 


conductivities are denoted by o = a, in the x-direction and o, = ke, in the y-direction, 


so that 


le] = 


the formula 

k l —_ k Cos (o D>) 

‘She _ T os (2) 
“ SIN Qs 


bo |+ 


2Y; O- = 


(and similarly for Y, and Y,) resulted, which for k = 1, i.e., isotropic materials, is identical 
with MacNeal’s Eq. (10). With a proper lay-out of the nodes it can, for tensor conduc- 
tivities also, be easily arranged that sides of triangles coincide with the boundaries and 
that all resulting conductances connecting two nodes are positive, i.e., physically real- 
izable. 

MacNeal’s paper embodies two remarkable advances in the solution of the problem 
in question. One is finding that an infinite multitude of simulating networks exists for 
given node points, and the other is the precise definition of the part of the continuum to 
be lumped at a node as far as the distributed external current sources are concerned. 
These advances can be extended to tensor conductivities, quite frequently met with 
in engineering, along the line of MacNeal’s method. 

In Fig. 2, which corresponds. to a part of Fig. 4 in MacNeal’s paper, let A and B 
be two nodes of the simulating network and 1 and 2 the corresponding points of the 
dual network so that the current J,, crossing the line 12 within the two-dimensional 
continuum divided by the potential difference V,, between the points A and B deter- 
mines the conductance Y,, of the resistor joining A and B. Let the direction 12 be 
found from the direction AB in the following way: Draw the line 34 perpendicular to AB. 
Multiply the distance N4 from the y-axis of a suitable point 4 of this line by 1/k. Make 
Nz’ parallel to Oz and N2 equal to 1/k-N4. Using MacNeal’s notation we obtain for 
I 4p 


a2 > 


Lay = —| odrn-{[]¥d} = —| odr{nlel}-Vo & —o-re | Vb | cosa (3) 


“1 “1 
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so that 


(4) 








Fic. 2. Portion of the asymmetrical network of triangles. 


This relation is general, whether 0 is the mid-point of AB or not. Equation 4 shows the 
splitting-up of Y,, into two partial conductances corresponding to the interior of the 
two triangles ABC and ABD. 

A closer investigation for 0 as the mid-point of AB is carried out in the notations 
of Fig. 3. The bisectors H,L, H.L and H;L of the sides of the triangle A,A,A; , with 


) 
ae Ax(*s,y3 
AS 












Az (~ , ys) 

















| A,(24,, y,) 
| od 


Fia. 3. Conditions for 0, 0’, and 0” of Fig. 2 being the mid-points of AB, BC and CA respectively. 
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the Cartesian co-ordinates of A, , A, and A, being (2, , y:), (x2 , y2) and (x; , ys) respec- 
tively, intersect at L. Let the lines H,M and H,M be the lines derived from the lines 
H,L and H,L by the procedure by which in Fig. 2 the line 12 was derived from the 
line 34. Let them intersect at the point M(x, , ys). Using basic methods of analytical 


geometry we obtain for yy 


- k(x, — ao)(%o — 2,)(a, — 23) — [y3(a2 — 2.) + yo(x, — as) + yi(t3 — 22)] (5 
uy = ; ‘ 
sid 3(Yo — Yi) + LAY. — Ys) + Ui(Ys — Yo) 


and a similar formula for x, . These formulae are symmetrical in the co-ordinates of 

A, , A, and A, , so that the line H.M derived from the line H,L by the same procedure 

passes through M. In other words, if the points 0, 0’ and 0” in Fig. 2 are the mid-points 

of the sides of the triangle A BC, the line 02 for the side AB and the corresponding lines 

0’2 and 0’’2 for the sides BC and CA intersect in one point 2. 
If (see Fig. 3) the ratio 


ol Yu — (Yr + y)/2 ; 
o — ‘ (0) 


yee Le — 2; 


r. = 


that corresponds to the ratio o-ro,/laz of Fig. 2 is computed, we obtain after some 
manipulating 
Vv o ko ~~ 
ay 5 = -+- Re (4) 
cot g. — cot ¢g, tan g, — tan ¢g 


Equation (7) is identical with Eq. (2), as can easily be shown. Hence Eq. (2) represents 
the special case of Eq. (4), if the point 0 (Fig. 2) and the corresponding other two points 
0’ and 0” are the mid-points of the sides of the triangle ABC. 

Tensor conductivities can be dealt with also by a change of variables [3]. This method 


leads to the same results. 
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IMPEDANCE SYNTHESIS WITHOUT MUTUAL COUPLING* 


By AARON FIALKOW (Polytechnic Institute of Brooklyn) 
AND IRVING GERST (Control Instrument Co.) 


In a fundamental contribution to network theory, Bott and Duffin [1] have given a 
method for the synthesis of an impedance which obviates the use of any mutual coupling 
*Received Feb. 16, 1954; revised manuscript received April 12, 1954. After submitting the paper for 
publication, we have learned that independent papers embodying substantially equivalent results but 
using different approaches from our own and each other were presented by F. Reza and R. H. Pantell 
at the January 1954 meeting of the A.I.E.E. 
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as was formerly required by the Brune realization technique. The purpose of this note 
is to present a modification of their procedure which results in fewer elements. 

The Bott-Duffin method depends upon the result [3] that if Z(p) is a positive real 
rational function (p.r.f.) of degree n, then 


_ p&p) — kZ(k) 


R(p) = pZ(k) — kZ(p) ’ k>0 (1) 


is again a p.r.f. of degree not greater than n. Using (1), they show that Z(p) is realized 
by a balanced bridge whose opposite pairs of arms are Z(k)R(p), Z(k)/R(p) and Z(k)p/k, 
kZ(k)/p respectively. If now it is assumed (after well-known preliminary reductions) 
that at p iwy , &) > 0, Z(iw) = iw, L > 0,' then by choosing k > 0 such that 


Z(k)/k L they obtain an R(p) which has a pole at p = +7 , so that 
R(p) = Rp) + P a : a> 0, (2) 


where F,(p) is a p.r.f. two degrees lower than R(p). Thus their algorithm requires six 
elements and two reduced functions. 

Our alternative synthesis of the above case utilizes an unbalanced bridge network 
lr whose pairs of opposite arms are Z, , Z, and Z, , Z; respectively, and whose bridging 
arm Z. connects the common node of Z, and Z, and the common node of Z, and Z, . 
We shall assume that Z, = W + S, Z, = 1/W where W and S are both p.r.f. If Z’(p) 
is the impedance of I, then (see [2, pp. 284-285]) 


Z'(p {Z,(Z. + Z;)W + [Z,S(Z, + Z:) + Z(Z2Z3 + 1) + 22+ 23 + Z5|W 
+ S(Z_ + Z3 + Zs) + Z(Zs + Zs)}/{(Z2 + Z3)W* + [S(Z. + Zs) 


+ ZZ; + 2.45 + 2:45 + 1\W+ 8S + Z;, + Zs} (3) 
We will now show how to choose Z, , Z; , Z; , S and a reduced function W so that 
Z'(p) = Z(p). It is desirable to impose the further condition that the right member of 


(3) reduce to a bilinear function of W. For this to be so, the resultant’ of the two quadratic 
polynomials which are the numerator and denominator of the fraction in (3) must be 
zero. After some calculation this leads to the equation 


— 2,2;)(4: Z; 
Sa See ae (4) 
Z;(Z2 + Zs) 

The simplest reactive choices of the Z’s in (4) which cause the resistors in W and 1/W 
to be eliminated at the frequency p = +iw, and which are consistent with Z(iw))/i > 0 
are 

, 1a ‘ ; 

Z,=—, Z,=bp, Z;,= ap, a> 0,b>0; 

P 
which gives 
g _ d= abw)(a + b)p 
’ ab(p* + ws) 

‘If L < 0, one considers either 1/Z(p) or Z(w?/p) instead of Z(p). 
The resultant of aoW? + aW + a and af{W? + afW + af is (afa: — a af) (af a, — aga{) — 


(aya, — afaz)*. 


’ 
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where a and b must be chosen so that 1 — abw, > 0. Substituting these values of the 
Z’s and S in (3), dividing out the common factor 
, y Ae + Zs 
W+5374 
* Ea. + 5 
and simplifying, we obtain 
abp(p* + wo)W + (a + b — ab’wo)p* + aw | (5) 
a(p’ + w.)W + plabp* + 1) . 
Now it follows from (1) and (2) that Z(p) may be written in terms of the reduced 
function R,(p) as 





Z"(p) = 


, i p(p° + o)R, + (a+ k)p” + kw ; 
Zp) = Ak) °F ree ie (¢ 
? (*) k(p + w)R, + pip + ak + wo) 6) 
A comparison of (5) and (6) shows that if we let W = R, ,b = 1/k, a = k/(ak + w), 
we obtain Z’(p) = Z(p)/Z(k). Making this impedance level change in Z’(p), we cenclude 
that Z(p) may be realized as an unbalanced bridge with 
_ Ak — abwo)(a + b)p. 


Tk: ( sS= nate 2 
Z(k)Ri(p), ab(p’ + wo) 


Z. = awZ(k)/p, Z; = bZ(k)p, Z, = Z(k)/R,(p), Z; = aZ(k)p, 


Z=W+S, W 


a = k/(ak + wi), b = I/k. 


Thus our algorithm requires five elements and two reduced functions. 
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ON TWO METHODS OF GENERATING SOLUTIONS OF LINEAR PARTIAL 
DIFFERENTIAL EQUATIONS BY MEANS OF DEFINITE INTEGRALS} 


By J. B. DIAZ* anv G. S. 8S. LUDFORD** 


(Institute for Fluid Dynamics and Applied Mathematics, University of Maryland) 


1. Definite integrals with variable limits of integration. Consider the canonical form 
of the linear hyperbolic differential equation in two independent variables 


L(u) = uy + a(x, yu. + b(a, yu, + c(z, yu = 0. (1) 
tReceived April 13, 1954. 
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It has been well known since the beginnings of the subject (see the detailed discussion 
in Goursat [4]) that if U(2, y, a) is any’ solution of (1) containing a parameter'a, for 
a, < @ < a, with a, and a, fixed constants, then 


u(x, y) = | fla)U(e, y, a) day, (2) 
where f(a) is an arbitrary function, is also a solution of (1). In this way solutions of (1) 
can be generated from functions f(a) of a single variable a. 

In 1895, Le Roux [1] showed that one or both of the fixed limits of integration can 
be replaced by a characteristic variable [say, x or y, in the case of (1)], and that (2) will 
still be a solution of (1), provided only that U(x, y, a) is chosen so as to satisfy certain 
additional conditions. Thus, if either limit of integration is replaced by x, then U(z, y, a) 
must satisfy the characteristic condition 


al + ala, y)U =0 On Z, = a, (3) 
oy 


which, incidentally, is one of the conditions satisfied by the Riemann function. Similarly, 
if either limit of integration is replaced by y, then U(z, y, a) must satisfy 


". + D(x, a)U = 0, on y = a. (4) 
As a matter of fact, the characteristic variables x and y, or an arbitrary function of 
just one of them, are the only variable limits which may be used in the definite integral 
in (2), if u(x, y) is to be again a solution of (1). This is proved by Le Roux [1] for the 
case in which the functions occurring in (2) are such that the definite integral may be 
differentiated under the integral sign. 

For instance, if R(é, n; x, y) is the Riemann function’ of (1) [so that, for fixed (&, 7), 
the function R(é, n; x, y) of the “pole variables” (x, y) is a solution of (1)] then, for 
n fixed, R(a, n; x, y) will be a solution of (1) satisfying (3); while, for & fixed, R(é, a; x, y) 
will be a solution of (1) satisfying (4); and, finally, R(a@, a; x, y) will be a solution of 
(1) satisfying both (3) and (4). Accordingly, for arbitrary f the definite integrals with 
variable limits 


az 


| fla)R(a, n; x, y) da, | f(@R(E, a; x, y) da, | S(a@)R(a, a; x, y) da, 


“ay ay 


will each be a solution of (1). 

2. The integral operator. In several papers Bergman’® has developed a theory of 
integral operators for obtaining solutions of partial differential equations. As applied 
to (1), the method is, for the most part (see the concluding paragraph of Sec. 3), based 


on the following theorem: 


'Throughout this paper precise statements concerning differentiability, continuity, and regions of 
definition of functions will be omitted, except in those cases where they are of prime importance. 
2For the connection between the integral operator of the next section and Riemann’s function, 


see [10]. 
’For a summary and bibliography see in particular [9] and [11]. 
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Theorem. Let E(x, y, t) be a solution of 

vr . l op , ; 7” 

(1 — f)(E,, + aE,) - , (E, + aE) + 22tL(E) = 0, (5) 


such that, for z ¥ 0, 
(1 — #)'’(E, + aE) z 
a ¥ (6) 

xt 
is continuous for ¢ = 0, and tends to zero for each (z, y) as ¢t approaches +1. Then, 
if f is an arbitrary once continuously differentiable function, the function u defined by 


a+1 

' : . . it 
u(x, y) = | E(x, y, Of[$e0 — ¢)] a — yz (7) 
J-1 i ) 


is a solution of (1). 
Proof. Writing, for brevity, f instead of f[42(1 — ¢’)], one has, upon differentiating (7), 


‘és . : ; Oe dt . 
L(u) = I {LE) -f + (LE, + aE) Lh ne A? . (8) 
Now, 
I 2, Of ya 
5 ( ~t +e eo = 0, 


and when zt + 0 this equation may be solved for df/dx. Substituting for df/dzx in (8) 
and integrating by parts (cf. (6) at ¢ = 0), one obtains 


w=] 47 — E Sy (E, + at) | by dt 


i- fr... . POPS te 
es ! Ort (E, + abs i 


The desired conclusion now follows from the hypothesis made about (6) as ¢ —~ +1, 


and the fact that 
L(E j=. fy" : 

firs -+- ! = -(E, + at) | = 0, 
ae or 


by virtue of (5). 

This theorem, as stated and proved above for Eq. (1), does not seem to be stated 
explicitly in full in the literature. However, for an equation occurring in gas dynamics 
[when a = b = 0, c¢ = c(x — y)], see Bergman [6]. Moreover, Ghaffari (8) considers the 
general equation (1), but does not state the essential conditions on the behavior of (6) 
at ¢ = 0. 

3. Definite integrals with variable limits of integration and the integral operator. 
The object of the present note is to show that the theorem of Sec. 2 and its obvious 
modification when x and y are complex (which constitute the kernel of Bergman’s 
operator method), are merely a restatement of part of Le Roux’s results mentioned in 
Sec. 1, and that Le Roux’s formulation is simpler in the fluid dynamical application. 
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First of all, notice that E(x, y, t) can be considered—without any loss of generality— 
to be an even function of ¢, for otherwise H(z, y, t) may be replaced by its even part 
1{H(x, y, t) + E(x, y, —8)] in (7), without altering the value of the definite integral. 
Hence, instead of (7), consider 

1 
u(x, y) = 2 f E(x, y, Of(420 — &)] a- }” , (9) 


Then, for z # 0, the change of integration variable 
1/2 
. ; . io 
a = x(1 — f), 1.e. {= (2=) ; (10) 


transforms Eq. (9) into 


; ez x —_ 1/2 ] 


Now, if U(x, y, a) is defined as 





x 
U(z, Y, a) 7 : ee ee 
’ [a(z — a)]'” 
so that 
ol’ 1 a 1 
= a \E. + =a aa Ee - 5 e}, 
OX la(a — | + on*"(x — a)’ ‘ 2(x — a) 
ou’ l 4 
, ri yi72 E, 9 
OY was = a) | 


ou 1 “a 
7 3 \Esy + 5-379 —-—sya En — 5 
Or OY alr — a) 2 { y 2° (x we a)'”? vt Hx am a) ufo 


then 
] f : a " 1 . 
L(U) =; a7a (L(E) + s-a78 iva (Ey. + aE.) — 5——s (E, + a€)?, 
eae a a) | y A (x = a) (x ame a) 
l f : a/x . . 
73 VL(E) +: 7s (Ey, + aE,) 
1a re — a) \ 2x[(x = a) x] 
l ; 
<a we hg E, LE \ 
2x|(x — a)/x] (E, + ab) 
= @. 


by (5), since the ¢ argument of EF is just [(r — a)/z]'”’. Hence U(x, y, a) is a solution 
of (1) containing a parameter a. Further 

au ; I ; a 

(oe + at ) a — 

\ Oy la(x — a)]} 
so that 


(a — 5 Ris » n  & (au ) 
ad (EZ, + ak) = x \ay + al}. 
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Now, the theorem uses the continuity, at ¢ = 0, of the function on the left hand side of 
the last equation. At first glance this appears to be a weaker requirement than the 
vanishing of the right hand side at a = zx, which is needed by Le Roux. However, as 
has already been remarked, EF can be assumed to be an even function of ¢, so that 
(1 — #)'?(E, + aE)/zt is odd in t; hence the continuity of the latter function at t = 0 
implies that it must of necessity vanish there. 

Thus the function F(z, y, t) used by Bergman may be written as the sum of an even 
part [E(z, y, t) + E(a, y, —t)]/2, which satisfies Le Roux’s condition (3) after the trans- 
formation (10) is performed; and an odd part [E(z, y, t) — E(x, y, —t)]/2, which plays 
no essential role, since it contributes nothing to the definite integral. This remark about 
the odd part of # has apparently been overlooked in working out some particular ex- 
amples. 

If the independent variables x and y in (1) are both complex, and the coefficients 
a, b, c are analytic functions of these two complex variables, then there is no essential 
distinction between equations of hyperbolic and elliptic type. In this case, which is 
the one usually considered by Bergman, the function /Z is taken to be analytic in the 
two complex variables x and y, and the function f (since it is differentiated with respect 
to z) is taken to be an analytic function of its (complex) argument. The formal calculation 
leading to the proof of the theorem is exactly the same as in the real case considered 
above. By this passage to complex z and y, however, the class of non-analytic solutions 
of a hyperbolic differential equation is excluded from direct consideration. It should 
also be remarked, see [5], that for certain purposes it is convenient, in the case when x 
and y are both complex, to consider ¢ as complex also, and to integrate from t = —1 
to ¢ = +1 along an arbitrary curve in the ¢-plane. If this curve does not pass through 
the origin, the Le Roux condition (3) need not be fulfilled, and hence, in these circum- 
stances, Bergman’s method is not equivalent to Le Roux’s. However, in concrete appli- 
cations the path of integration from ¢ = —1 tot = +1 seems to be invariably taken 
along the real axis, as in the theorem of Sec. 2. 

4. The application of the integral operator method to fluid dynamics. 
major applications of Bergman’s integral operator method is to the equations (in a 
“modified hodograph” plane) governing the plane steady irrotational flow of a perfect 
compressible fluid. Under a suitable change of independent variables the equation for, 
say, the stream function y, is of the type (1), with coefficients a, b, c depending on the 


One of the 


difference x — y alone. 

The method of Le 
which satisfies either the condition (3) on the single characteristic x = 0, or the condition 
(4) on the single characteristic y = 0, or both. Once this solution U*(z, y) has been 


toux reduces in this case to finding a solution U*(z, y) of (1) 


found, one sets 


U(z, y, a) = U*(x — a, y — a), (11) 


and then a family of solutions depending on one arbitrary function is given at once by 
(2), with a, and/or a, suitably replaced by x and/or y, as the case may be. For, since 
the partial differential equation is invariant when the variables x and y are replaced by 
x — aand y — a respectively, the function U(z, y, a) defined by (11) is again a solution, 
and will satisfy the corresponding Le Roux conditions on x = a and/or y = a. 

A special case of Eq. (1) with coefficients depending on x — y alone is the Euler- 


Poisson equation (see Darboux [2]) 
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tty — B 2+ B 4, = 0, (12) 


u v 
z—-y z—-y 





which for 8 = 8’ = 1/6 serves as a first approximation (Tricomi) to the equation for 
the stream function of compressible flow near the sonic line. Darboux uses the procedure 
just outlined in the previous paragraph, and an argument equivalent to that of Le Roux, 
to construct the general solution of Eq. (12); and from this, Tricomi [3], has constructed 
a simple definite integral representation for the solution of the Cauchy problem for 
(12) with 8 = #’ = 1/6, the function u and its normal derivative being given on the 
sonic line x = y. 

Besides the general compressibility equation, Bergman considers (in his terminology) 
a “simplified’’ compressibility equation [7]. In effect this “simplified” equation is that 
of Tricomi, under change of variables. In dealing with this “simplified” equation, he 
first finds solutions of the corresponding equation (5) for E in terms of hypergeometric 
functions, and is led to a complicated representation of a class of solutions of this equation. 

Thus, even for this “simplified” equation, it appears that Bergman’s method is 
more complicated than that of Le Roux. Indeed, Bergman’s method consists in finding 
particular solutions of an associated equation for EZ, Eq. (5), in three independent 
variables, whereas the method of Le Roux consists in determining particular solutions 
of the originally given differential equation in two independent variables. 

The complications involved are apparent from the following concluding remark: 
Bergman looks for particular functions E(x, y, t) which involve series whose general 
term is a power of ¢’x times a function of x — y. This is easily seen to be equivalent, in 
view of the change of variable (10), to just looking for particular solutions of the original 
differential equation using series whose general term is a power of (x — a) times a func- 
tion of  — y, or (because of the invariance of the differential equation when z and y 
are replaced by x — a and y — a respectively), just a power of x times a function of 


r— ¥. 
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AN INTEGRAL EQUATION GOVERNING ELECTROMAGNETIC WAVES* 
By P. R. GARABEDIAN (Stanford University) 


1. Preliminaries. We shall consider in this paper the Helmholtz equation 
Au+ ku =0 (1) 


in two independent variables £ and 7. In electromagnetic theory the problem arises of 
determining, in the exterior D of a simple closed curve C, a solution wu of (1) with pre- 
scribed boundary values or prescribed normal derivatives on C. The function u is supposed 
in addition to fulfil at infinity the radiation condition 


, {Ou . 
lim r' (2 — itu) == (), (2) 
or 


r-@ 


where r = (’ + n’)'” is the distance from the origin. Thus u will be complex-valued 
and u(é, n)e““’ will be the actual wave-function for the problem considered. We can 
think of « as representing a two-dimensional wave scattered by the finite object C, 
and the type of boundary condition to be imposed along C then depends on whether 
we treat the electric or the magnetic field. 

We shall investigate in detail the case when the boundary values of u are given 
along C. However, all our statements carry over almost verbatim when it is the normal 
derivatives which are prescribed instead. The uniqueness of the solution of (1) and (2) 
follows easily from the work of Rellich [2], but the existence of the solution as presented 
by Weyl [3] and Miiller [1] depends on a careful discussion of the eigenfunctions of (1) 
for the interior of C. The object of the present note is to reduce the boundary value 
problem for (1) in the exterior domain D to a new Fredholm integral equation whose 
study is independent of the interior domain and its eigenvalues. This will be accom- 
plished by introducing a suitable parametrix for the problem, constructed by means of 
conformal mapping. Our integral equation has, furthermore, the advantage that it can 
be solved numerically by iteration in certain important cases. 

The middle section of the paper is devoted to the development of this integral 
equation. In the concluding section, we discuss the uniqueness question from a new 
point of view and indicate a variational formula for the estimation of the scattering 
cross section of C. 

We shall assume that the frequency k = 1 in the following without loss of generality, 
since this reduction corresponds merely to a change of scale. 

2. The parametrix method. Our analysis of the equation (1) is based on the notion 
of a parametrix S(¢, 7). Such a parametrix S is not a solution of (1), but it satisfies the 


following requirements. It is regular as a function of ¢ in D, except at ¢ = 17, where 


S(¢, 7) + log |¢— 7 
remains continuous together with its first partial derivatives. At infinity S satisfies 


; 12 OS i 
lim r § _ is) = () 
or 


r-© 


the radiation condition 


*Received March 1, 1954. This work was done under a contract with the Office of Naval Research. 
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as a function of ¢, and along C we have S = 0 for each fixed 7 in D. 

Suppose that wu is any solution of (1) and (2) in D with k = 1, as was agreed above. 
Then from the radiation condition we obtain 


; aS , Ou 
lim [ (. — — S- “) ds = 0, 
~e Jigar on on 


r 1 


where s and n represent arc length and inner normal along the path of integration. 
Hence by Green’s theorem and by (1), 


; IS(t, 7 ; / 
| u(t) - u 27) ds = 2ru(r) — [| u(AS + S) dé dn. (3) 
J on Je 
D 


There are many elementary ways in which S can be constructed, and in each case (3) 
provides an integral equation for the determination of w. 
We choose S in a special and particularly useful manner. Let 
f=2z+ f@ (4) 
be the conformal transformation of the exterior | z | > FR of a circle of radius R onto D, 


normalized so that 


Am = 
f(z) => > mm (5) 
m=0 «© 
is regular at infinity. Furthermore, denote by G(z, w) the Green’s function for (1) in 
the region | z | > R. To be precise, G satisfies the equation 
AG +G=0 (6) 
as a function of z, except at z = w, where 


G(z, w) + log |z -— w 


remains continuous; G@ satisfies the radiation condition at ©; and G = 0 for|z{| = R. 
This particular Green’s function has a variety of explicit representations in terms of 


Bessel functions, all of which are obtained by the standard procedure of separation of 
variables in polar coordinates. 
With ¢ given by (4) and with r = w + f(w), we set 
S(¢, 7) = Gee, w). (7) 


The function S so constructed by application of a conformal transformation on the 
arguments of the known Green’s function G is evidently a parametrix of the type defined 
above. Indeed, S is regular in D except for a logarithmic infinity at ¢ = +, S vanishes 
on C, and S even satisfies the radiation condition, since the conformal transformation 
(4) has a derivative 1 at infinity and therefore leaves this condition invariant. 

We substitute (7) into (3) and evaluate the integrals in the z-plane. The area element 
dé dn and the Laplacian AS are altered in the transformation by multiplication and 
division, respectively, with the Jacobian 


0(é, n) zi ' 
= ae 
Page 1+ f"(2) | 








430 NOTES (Vol. XII, No. 4 


In particular, the expression AS dé dy is a conformal invariant, as is also (@S/dn) ds. 
Therefore we obtain 
G oes ae a 
[ _ Ue) mee Das = 2U(w) — [[ V@l~aG+a@|14+f’ Pldzdy, 


where U(w) stands for u(r) and likewise 

U(z) = u(f). (8) 
Application of (6) yields 
r aG Ze P > 19 
| U(z) fee ds = 2rU(w) — | U(z)G(z, w){2 Re f’(z) + | f’(2) |} dx dy. (9) 
¥\izi=R « 

(zi>R 


In order to derive from (9) a Fredholm integral equation for U with a symmetric 
kernel, we introduce the notations 


plz) = {2 Re f(z) + | f’@ |)" (10) 
V(z) = p(z)U(2), (11) 
K(z, w) = = (z)p(w)G(z, w) (12) 


We notice that when u(f) is prescribed on C, the expression 


g(w) = mf u@ 2) g 8 (13) 


will be known. Multiplying (9) on both sides by p(w) /2z, we arrive at the final Fredholm 
integral equation 
g(w) = V(w) — || V(z)K(z, w) dx dy (14) 
[2z|>R 


for V, or in other words for U, and this is the basis of our entire discussion. 
With w = a + ib, the symmetric kernel K of (14) satisfies the square-integrability 
condition 


I/ I K(z, w) |? dx dy da db 


|1>R z R 
= = I/ /| G(z, w) |? | p(z)p(w) |? dx dy da db < @, (15) 
T J 
Iwl>R |2|>R 
since | z | p(z) is bounded at infinity, by (5) and (10), and since | z |'”*| w |'”G(z, w) is 
bounded outside a fixed selabdveteoed of z = w, by standard vahediee of the Bessel 


functions. Hence the Fredholm theory is A applicable to (14), and we deduce 
that for any given g we can find a unique square-integrable solution V, provided that 
the homogeneous equation 


I VK dz dy (16) 
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has no non-trivial eigenfunction V. But an eigenfunction V transforms by (11) and (8) 
into a solution u of (1) which vanishes on C and satisfies the radiation condition (2), 
as can easily be verified by direct calculation from (16). Rellich’s uniqueness theorem [2] 
then shows that u = 0, whence V = 0 and there are, indeed, no eigenfunctions. 

Thus, given u on C, we can calculate g, substitute into (14), solve for V, and find 
u as a solution of (1) and (2) by using (8) and (11). This completes our proof of the 
existence of the solution of the boundary value problem. The hypotheses on the smooth- 
ness of C which are required in the argument are only those needed for the demonstration 
of the uniqueness theorem. In particular, our integral equation yields the solution of the 
problem for curves C with a finite number of corners, and the behavior of u at the corners 
can be derived from (14). Furthermore, the same proof is valid when it is the normal 
derivatives of u which are assigned on C. For this case, we have only to replace the Green’s 
function G throughout the argument by the corresponding Neumann’s function N(z, w) 
for (1) and (2) in the region | z | > R. The function N is characterized by the same 
requirements as we imposed upon G, except that we replace the boundary condition 
G = 0 by the boundary condition dN/dn = 0 along | z | = R. For the circle, the Neu- 
mann’s function has an explicit expansion in terms of Bessel functions, obtained by 
separation of variables. 

From the practical standpoint, the integral equation (14) has some interest. If the 
domain D differs sufficiently little from the exterior of a circle, the function p depending 
on the conformal transformation (4) will be so small that 


[/ | K(z, w) | dx dy Se <1. (17) 


l2z|>R 


Hence the smallest eigenvalue of (14) will have a modulus exceeding 1 and the resolvent 
kernel can be used to solve the equation in a Neumann-Liouville series. In fact, the 
successive approximations defined by 


Ve=gt I V,.iK dx dy (18) 
l2zI>R 


converge geometrically when (17) holds. The a priori bound 


» | | 
Lat ie oe (19) 





which follows immediately in this case, is also significant, since the maximum principle 
fails for the partial differential equation (1) and more elementary estimates of u are 
therefore not available. 

The case (17) of convergence of the successive approximations (18) to the solution 
V of (14) has an implication for the existence of solutions of the boundary value problem 
for (1), since no hypothesis of smoothness on the boundary curve C is required when 
(17) holds. Finally, we point out that in the development of our method, the conformal 
transformation (4) and the explicit Green’s function G could be replaced by the corre- 
sponding quantities for any region in which one knows how to solve the equation (1). 
3. Remarks about uniqueness and the scattering cross section. For solutions u of 
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the Helmholtz equation in D not necessarily satisfying the radiation condition, one can 
define a norm by the formula 
2 ° i ) 2 | ( 9 
u = lim ) ul + i=! pd. (20) 
If wu. is the solution of (1) and (2) with boundary values h on C and if w_ is the solution 
of (1) and (2) with boundary values h on C, it then follows that the function 


U, + U_ (21) 


9) 


Uo 


is characterized by the property that among all solutions u of (1) with the boundary 
values A on C it has the smallest norm, 


Uo = min || u ||. (22) 
Indeed, writing v = u — wu, we find 
u ||° = || wo Il? + lle II’, 
since by the radiation condition 
° ( P a, » | P aia ad 
J. OUy OD l ) (ou, Ou _. @ 
lim | up) + ——;ds = = a _--— — (u, — @) = ds 
eee ~~ on anj 21 Je on on On 
= 0. 


The radiation condition also yields for the minimum value of the norm the relation 


+ 


u, ||? = lim | u, |° ds 


. i 9 
= lim | lu_ |’ ds. 


| 


The extremal characterization (22) of uw) is a refinement of the uniqueness theorem 
for the solution u of (1) and (2) with given boundary values on C. For if h = 0 then 

Up = (0 and therefore u, = 0, whence if wu satisfies (2) and vanishes on C it must 
vanish identically, since we can put both u, = u_ = wand u. = u_ = twin (21). 

A case of special interest for (22) occurs when h = —e‘*. The function u, then repre- 
sents the scattered field due to a wave of the form e’’ incident on C, while u_ represents 
the scattered field resulting from an incident wave e '*. The estimate (22) shows that 
the total scattering cross section o of C is given by 


o = min U \\” (23) 


among all solutions u of (1) with the boundary values —e’*. This characterization of 
indicates clearly the dependence of the scattering cross section on both u, and u_ , 
in view of the form (21) of the extremal function. 

Another formula exhibiting the dependence of o on both wu, and w_ is obtained when 
we shift C by an infinitesimal amount én along its inner normal and attempt to express 
the scattering cross section o* of the shifted curve C* in terms of quantities associated 
with C. We continue to take h —e'* and we introduce the total fields 


O =e +t; 


GQ. = C . + U_ , 
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which vanish on C. We denote by ¢* and ¢* , and so forth, the corresponding quantities 
associated with C*, and we obtain by an easy application of Green’s theorem and the 


radiation condition 








fugull [ [lut ? — | uy [] ds 
ree J\i¢el=r 
_. out _ Ou 
= gm [ u* ds — gm] a, —ds 
Joe on c on 


dg* do_ 
= jm [ u* 2+ ds — gm / Us “P= ds 
Joe on c on 


dg*t du* do aus 
= J * —* _ ot —— — — : 
- il ( on ?+ “an oar 7? an ds 


iz ! x ais 
= gm [ (ue Ge — e** du* — U. de ~_ + a dus) ds 








- an an an on 
dy. de* 
= gm [ (6 sy — 9, ae*) ds. 


In the last integral we are free to choose C as the path of integration, and since yg, = 0 
on C and ¢* = 0 on (*, we find 


o* —a 


de. F Og, ) 
g [——a=— | c—— - 
mifie an das | an os | 


= —3m I {VotVe. — gte.} dx dy, 
D-D* 
where D* is the exterior of C*. We now consider only terms of the first order in the in- 
finitesimal shift 5n and derive, with the notation o* — o = éc, the Hadamard variational 


formula 
a dp- dg. ‘ 9 
ba = —Im [ a bn ds (24) 


for the scattering cross section o. The result is valid for small shifts én of either sign. 

Formula (24) has an obvious importance for the estimation of the scattering cross 
section o of figures near those for which o has been studied. As a special application of 
(24), we remark that it substantiates the conjecture that a vertical segment has the 
least scattering cross section among all simple closed curves which enclose its end-points. 
In fact, according to (24) we check that by symmetry 6¢ = 0 for any normal shift of 
the vertical segment carrying it into a neighboring curve with the same end-points. 

Finally, we call attention to the fact that only the most elementary alterations are 
required in order to extend all the remarks of this section to the case in which the normal 
derivatives of u are prescribed in the boundary conditions along C. A generalization to 
space of three dimensions is likewise possible. 
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STABILITY OF TWO-DIMENSIONAL PARALLEL FLOWS 
FOR THREE-DIMENSIONAL DISTURBANCES* 


By CHIA-SHUN YIH (State University of Iowa) 


The object of this note is to establish a relationship between the stability of two- 
dimensional parallel flows for three-dimensional disturbances and that for two-dimen- 
sional ones. The special case of confined flow of a homogeneous fluid has been considered 
by Squire’. In the present note neither is the upper surface of the fluid necessarily assumed 
to be fixed, nor are the gravitational force and variations in density and viscosity neg- 
lected. The variations in density and viscosity, which for two-dimensional flow can 
occur only in the direction normal to the plane boundary along which the fluid flows, 
may be continuous or discontinuous. 

For the disturbance, a stream function of the type 


vy’ = oly) exp i(mz + nz — met) (1) 
can be taken, in which z is measured in the direction of the primary flow, y is measured 
in the direction normal to the plane boundary, z is taken along an axis normal to the 
x, y plane, and m, n, and ¢ are constants. If a rotation about the y-axis is performed 
so that the z’-axis has the direction numbers (m, 0, n) with respect to the original 


coordinate system, then 
mx + nz = m’x’ (2) 


in which m’ = 
written as 


(m? + n’)'”. If, furthermore, c’ is defined by m’c’ = mc, Eq. (1) can be 


vy’ = g(y) exp im’(x’ — c’l) (3) 


which represents a two-dimensional disturbance progressing in the 2z’-direction with 
wave number m’ and celerity c’. Similarly, the stream function 


y’’ = o(y) exp i(mx — nz — met) = ¢g{y) exp im’(x’’ — c’t) (4) 
represents the same disturbance progressing in the x’’-direction with direction numbers 


(m, 0, —n). 
Now the stream function 


y=’ +” = 2(y) cos (nz) exp im(x — ct) (5) 


represents a three-dimensional disturbance progressing in the z-direction with celerity 
c and having wave numbers m and n in the z- and 2-directions, respectively. Since by 
symmetry y’ and y” are physically identical, it follows that the flow is stable or unstable 
for y according as it is stable or unstable for y’, because the differential system governing 
stability is linear and homogeneous in the quantity representing the stream function 


of the disturbance, so that the rule of superposition applies. 
Now if the Reynolds number, slope, and pressure gradient of the primary flow are 


*Received February 10, 1954. 
1H. B. Squire, On the stability for three-dimensional disturbances of viscous fluid flow between parallel 


walls, Proc. Roy. Soc. London A142, 621-628 (1933). 
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denoted respectively by R, s, and dp/dx, those for that component of the primary flow 
pertinent to y’, denoted by R’, s’, and dp’/dx’, are determined by 

m'R’ = mR, m’s’ = ms, csc B’ = 2 esc 6 (6) 

Ox Ox 

in which 6 and #’ are the angles of inclination of the boundary to the horizontal in the 
x- and 2’-directions, respectively, so that s = tan 8 and s’ = tan @’. As can be easily 
seen, the cross flow (in the z’-direction) for the case of ¥’ makes no consequential con- 
tribution to either the equations of motion or the equation of continuity, and does not 
affect in any way the satisfaction of the boundary conditions. Hence, the primary flow 
is stable or unstable for a three-dimensional disturbance according as it is stable or 
unstable for a two-dimensional one at a lower Reynolds number, a milder slope, and a 
reduced pressure gradient: the laws of reduction for the three quantities being given by 
Eqs. (6). 

The writer arrived at the foregoing conclusion by using Squire’s approach. Whereas 
many details were obtained as by-products, the rather cumbersome calculations involved 
are unnecessary if only the principal conclusion is desired. The general approach used 
here was suggested by the reviewer of the original manuscript of this paper, who attri- 
buted it to Professor C. C. Lin of the Massachusetts Institute of Technology. 


ON DIFFUSION IN AN EXTERNAL FIELD AND THE ADJOINT SOURCE 
PROBLEM* 


By JULIAN KEILSON (Lincoln Laboratory, M.I.T.) 


Abstract. If diffusion in an external field is described by dp/8t = DV’*p — p/r — 
V -(F(r)p), the function y(ro) describing the probability that a particle at r, will reach 
a collector surface before decaying or being absorbed by other surfaces satisfies the 
equation DV*y — y/r + F(t)-Vy = 0. This equation has no singularity to disturb 
any geometric symmetry available. Boundary conditions on y(r) at the collector surface 
and other influencing surfaces are derived and shown to be independent of the external 
field. The boundary conditions at the secondary surfaces are homogeneous. The collector 
surface boundary condition is inhomogeneous. 

1. Introduction. If particles with lifetime 7 diffuse in the presence of an external 
field F(r), changes in density p (r, ¢) are described by the continuity equation’ 


Op _ «PP 
ot T V-i. (1) 


Here j, the current density, is given by 
j= —DVp + Fae (2) 


*Received April 17, 1953; revised manuscript received Dec. 31, 1953. The research in this document 
was supported jointly by the Army, Navy, and Air Force under contract with the Massachusetts Insti- 
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1J. Keilson, On the diffusion of decaying particles in a radial electric field, Journ. Appl. Phys., 24, 1397- 
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436 
so the equation of motion is 


Op 2 p 
— = DV'p —-— — V-(Fir)p). (3) 
ot T 

Suppose one is interested in the probability y(ro) that a particle at ro will reach the 
collector surface S before decaying or being absorbed by surface S’ (Fig. 1). 


Ss! 
Fic. 1 
To obtain y(r.) one could solve the steady state diffusion equation for a source 


I, at To 


2 p 5 
DV*p —-~— V-(F@)p) = —L,6r — Pr) (4) 
T 
with given boundary conditions at surfaces S and S’, and obtain y(ro) as 
lf, \. Pe 
Y("o) = Tl. | ; —DV p oa F(r) p} -de. (5d) 
Os 


It is often simpler to obtain y(r,) as the solution of the homogeneous equation adjoint 


to (4): 
DV*y — — + Fr): Vy = 0. (6) 
The singularity at r, no longer interferes with any geometric symmetry available, 
and the boundary conditions on y, it will be shown, no longer contain the field F(r). 
The boundary conditions at the collector surface, however, are now inhomogeneous. 
That 7 obeys Eq. (6) may be seen in the following way. Let P(r,/r; t) be the solution 
of the time dependent diffusion equation (3), ¢ seconds after the particles were known 
to be at r, , i.e., that solution for which 
P(r,/r; 0) = ot — rp). (7) 
Continuity of probability demands that, for all ¢, 


t 


y(ro) = | P(t,/r; ty(x) dr + | ds | Jas , 8) do. (8) 
/0 ; 


“Ss 


jn(fs , 8) is the normal component of the current density at the point rs on our collector 
surface S at time s. At t = 0, the collector current density hasn’t had a chance to build 


up and j,(rs , 0) = 0. 
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If then one differentiates Eq. (8) with respect to ¢ at ¢ = 0, it is seen that 


Op ; 
/ ry (r,/r; 0)y(r) dr = 0, 


rf 6(r — fo) 


| { DY é6(r — 1%) -— - — V-(F@ dr — rd) i dr = 0. 


When this is integrated by parts, Eq. (6) is obtained. 
2. Boundary conditions. Suppose the boundary conditions on p are 


i, = —D 2% + Fa)p = ap at S, (9) 
on 
and 


a : ‘ 
= + F.(x)p = a’pat S’. (10) 


The corresponding boundary conditions on y(r) at S, and S’ are wanted. To find 
them one multiplies (3) by y(r), integrates over space and time, and obtains 


nt fd 9P t e P 
| ds | y(r - (r,/r;s) dr = [ ds / (| DV'P ---— v-@r)} dr. 
“0 e Cc /0 ¢ 
If one uses (6) and (7) and carries out integration by parts for both space and time 
one finds 


at 


. 


ds | 4 


5 


- >/ . ( » OY oP "id 
y(%o) = | P(r,/r; tyr) dr + | {1 _- Dy an +. Palt)P} do 


“0 


at 


os | ds [ 


“0 7 §* 


ay aP } 
2 ~ — —_ J or 
{pI re Dy = + yF(ts))P? de. 
But Eq. (8) tells us that 
. at . a) ) 
v(f,) = | P(to/rjit)y(e) dr + ds | {pe + PP} de. 


Comparing our last two equations it is seen that one must have for surface S’: 


Oy oP ’ 

2 aoe ne a G Do 
DI an Dy a7 + yF,P = 0, 
and hence 


D st ro ery (11) 


This means that if, at S’, one had perfect reflection (a’ = 0), the condition on y there 
would be dy/dn = 0. If S’ is perfectly absorbing (a’ = ©), our condition on y is that 
+ = 0. This is what one would expect since a particle near such a surface has little hope 


of reaching the collector surface. 
At the collector surface, to agree with the continuity of probability equation one must 


have: 


_ poe — pp — py wear p 
od on +P = on DY on + YF ub 
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and hence 
a=D A + ay. (12) 
on 

This condition states that at a perfectly absorbent collector surface (a2 = ©), y = 1. 
If the surface is perfectly reflecting dy/dn = 0, which is consistent with the more natural 
condition y = 0. 

A typical application of the formalism above would be to the diffusion of decaying 
particles in the presence of a radial electric field. In the earlier paper’ the. formalism is 
applied to the very simple case of a perfectly absorbing sphere about the center of 
symmetry. 

I would like to thank Dr. 8. F. Neustadter for his interest and encouragement. 


SOME PROPERTIES OF OPTIMAL LINEAR FILTERS' 
By HERBERT A. SIMON 


(Carnegie Institute of Technology and Cowles Commission for Research in Economics) 


The purpose of this note is to discuss two considerations that arise in designing a 
filter that is optimal in the sense of minimizing the integral over time of some norm, 
under the further condition that we do not wish to employ information about the statisti- 
cal structure of the signal. This is the “classical” problem of filter design prior to Wiener’s 
work on stationary time series, although, in general, design was and is guided by figures 
of merit rather than by an explicit minimization procedure. It is probable that the 
results set forth here are well known, at least intuitively, to designers of servomechanisms, 
but I have not seen them set forth. 

Since the norm to be miniminized is an integral over time, the problem may be 
viewed as one in the calculus of variations. Suppose that the norm, F, is: 


ats 


F=] $y, yl, +++ Ym} at, (1) 


to 


where the integrand, ¢, is a polynomial in some variable, y(¢), and its first n derivatives, 
Yun) (t) = d*/dt' y(t), (k = 1, ---, n); so that we may write: 


() =C+ Do dbyw + DD aiyoys + VO, (2) 


i=0 i=0 j=0 
where y/(¢) is a polynomial composed of terms of third and higher degree in y(¢) and its 


derivatives. 
To obtain a necessary condition for a minimum, we form Euler’s equation: 


d , 
Po = ‘The + tae + (—1) dt" Pn — 0, (3) 


where ¢, = 0¢6/0Y) . 


1Received April 28, 1954. The results presented in this paper were obtained in the course of research 
undertaken by the Cowles Commission for Research in Economics under contract Nonr-358(01), NR 
047-006 with the office of Naval Research. I am indebted to my colleagues A. Charnes, C. C. Holt, and 
F. Modigliani with whom I have had numerous valuable discussions on the topic of this note. 
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Equation (3) is a differential equation of order not greater than 2n in y(t). We now 


make the assumption: 
(A). Equation (3) is a linear differential equation with constant coefficients. From (2), 


we obtain: 


& =b+ p nyo + > Ai Yoiy + Vld) (k =0,--- ,n), (4) 
at n n a‘ 
qe Pe = b QiYcirry + > Qi Yoirey + dt" v(t) (kK = 1, +++ ,n). (5) 


Substituting (4) and (5) in (3), we have, under assumption (A): 


b, + + AiYuy + p ® Q;Y) + Wolb) 


n j m n d* 
- 2e (—1) { be Sam + 2, atu > Ti wco} 


i=0 


= AY + AY) + -+- + AenY(2n) = 0. (6) 


Rearranging the terms of (6), we obtain: 


bo + z, = } (QiYoise + auiteon)} 


t=O s=0 


+ vol) + (=I Sa WO — Dat = 0. 7) 


Now, since all the terms of y are of third or higher degree in y and its derivatives, 
all the non-vanishing terms of y, will be of second or higher degree, as will also all the 
terms of d‘y,/dt" . Hence, the linearity of the differential equation implies that: 


ol) + ("Se lO) = 0. 8) 


A necessary and sufficient condition for (8) [see Courant-Hilbert, vol. I, p. 167] is 
that: 


| Vy, siti » Yin) dt = E(y, ae » Yin-1))s (9) 


vt 


i.e., that w is independent of the path and hence does not affect the optimum. 


We have proved the 
Theorem I: A necessary and sufficient condition that the extremal corresponding to 


(2) satisfies assumption (A) is that: 


F = | E + Vbyw + DD anove | dt + &(y, *-* 5 Yn-))- (10) 


t=O 3=0 


Frgm (7) we can deduce: 
Theorem II: In >>?" a,Yc), all a, must vanish for odd g. 
Proof: Since all summations in (7) run from zero to n, for each term (—1)* (aj, + ai) 


Yccae) We have a corresponding term (—1)' (ai; + @ix) Yass - If (¢ + &) is an odd number 
these terms will have opposite sign and, being otherwise identical, will vanish. 
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Theorem II has an important consequence. Form the characteristic equation for (7): 


> anp™ = 0. (11) 
h=0 


Because this is an even function, if p, is a root, —p, is also a root. Hence, either (1) 
all the roots are pure imaginaries, or (2) there is at least one root with positive real part. 
If we interpret the system described by the Euler equation as a linear “filter” or decision 
rule, the behavior produced by application of the rule will be dynamically unstable. 

The two theorems we have established may be interpreted as follows: 

1. From Theorem I it appears that linear filters will have desirable qualities when 
the norm or “error’’ we wish to minimize is quadratic. If the norm is decidedly non- 
quadratic this suggests that we can improve filter performance by introducing appropri- 
ate non-linearities. 

2. From Theorem II and the consequence derived from it, it appears that straight- 
forward application of the calculus of variations to the filter design problem leads to 
the prescription of an unstable filter, and hence is not practicable. 

It can easily be shown that point (2) is related to the fact that the Euler equations 
give only a necessary and not a sufficient condition for a minimum. Hence it does not 
follow that a path, y(¢), that satisfies (3) will thereby minimize (1). The specification of 
appropriate initial and terminal conditions to guarantee a bona fide minimum requires 
information about the future of the signal, and therefore, in the face of an incompletely 
predictable signal the method breaks down. 

The difficulty cannot be avoided by minimizing in the domain of the Laplace trans- 
form of ¢ instead of in the time domain. For, in general, by Parseval’s theorem:’ 


a oO 


[ (o(t)]? dt = | [¢*(p)]’ dp, (12) 
where ¢*(p) is the transform of ¢(t). 
Hence, any function that is a stationary value for the right-hand side of (12) is 
simply the transform of the function that is a stationary value for the left-hand side. 
It is beyond the scope of the present note to discuss newer methods that avoid this 
difficulty without reverting to cut-and-try procedures or the use of arbitrary figures of 


merit. 
*Professor A. Charnes pointed out to me the relevance of Parseval’s theorem for this problem. 
Minimization in the domain of the transform is essentially the method of Ritz for the “direct” solution 


of variational problems (Courant-Hilbert, I: p. 150). 


BOOK REVIEWS 


History of strength of materials with a brief account of the history of theory of elasticity and 
theory of structures. By S. P. Timoshenko. McGraw-Hill Book Co., Inc., New York, 


Toronto, London. 1953. x + 452 pp. $10.00. 


Professor Timoshenko’s original contributions to the subject of Strength of Materials are outstanding 
in contemporary literature; the same is true of his various books on Theoretical Mechanics. It is therefore 
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with considerable interest that a new book by him is read by persons sharing his interests. This is perhaps 
even more true of the present book. The author has delayed writing it until he had the opportunity to 
devote an uninterrupted period of time in uncovering all the necessary historical material. 

The book is founded upon lectures on the History of Strength of Materials which the author has 
given over the last twenty-five years to engineering students with some knowledge of the strength of 
materials and theory of structures. In writing this book the author has had in mind not only such students 
but also other people whose immediate interests are not specifically in this field. The histories of the three 
subjects, strength of materials, theory of elasticity, and theory of structures, are so closely interwoven as 
to be almost inseparable from each other. Much has been written in this historical field from various 
points of view. The author has followed the example of Saint-Venant in his “Historique Abrégé. . .” 
rather than that of Todhunter and Pearson in their ‘‘A History of the Elasticity and Strength of Mate- 
rials”. The result is an informal account of the history of strength of materials. Included are closely- 
related parts of the histories of the theory of elasticity and the theory of structures, purely mathematical 
developments being omitted from the former and purely technical developments from the latter. The 
presentation, with few deviations, is chronological. Seven chapters are devoted to the strength of materi- 
als, two to the theory of structures, and four to the theory of elasticity. These last chapters omit the 
mention of some well-known names, such as Volterra, and should not be read separately. 

A brief introduction describes the empirical approach to the strength of materials current from 
Ancient Times right down to the Renaissance, and includes comments on the failure of early Renaissance 
engineers to appreciate da Vinci’s work. The story proper begins with Galileo and unfolds through many 
equally famous names right down to the present time. The wealth of material may be judged from these 
statistics. About 675 names are mentioned in the text, and biographical material is given for nearly 80 of 
these. There are 245 figures, and about 40 of these are reproductions of likenesses of famous men. The 
reader can hardly fail to be stimulated by the biographical accounts often enlivened by anecdotes. The 
presentation inevitably presents a more difficult task as knowledge accumulates through the efforts of an 
ever-increasing body of research workers. This is especially true of the period 1900-50 dealt with in the 
final three chapters. This period is of course the one with which the author has been intimately concerned, 
and the developments are ably described. Here the coverage is extremely wide; thus fatigue of metals, 
approximate methods of solving elasticity problems, and the theory of ship structures are some examples 
of topics discussed. Many original references are cited throughout the book, and discussion is often given 
of the influence of great teachers. Comments are also made on the circumstances that either stimulated 
or retarded progress in the strength of materials in different countries over similar periods of time. 

A stranger to the subject would not suspect, although outside the covers of this book it is well- 
recognized, that the genesis of the modern development of the strength of materials in the United States 
is largely attributable to Professor Timoshenko himself. In the text the author simply refers briefly to 
himself as a pupil of Prandtl; his name does not occur in the index, and is relegated exclusively to footnote 
references to original papers. 

The reviewer agrees with the author that there is a place for formal lectures on such historical devel- 
opments in any well-balanced engineering curriculum. For pedagogic reasons it is often not possible to 
follow the chronological development of a subject in a course of lectures to students. Engineering students 
must therefore often fail to appreciate the painstaking way in which our present knowledge of the strength 
of materials has been gradually built-up since Galileo’s time. The current plea is a reflection of the need 
for the increased teaching of the history and philosophy of the natural sciences in undergraduate curricula. 
The reviewer regrets the author’s strict preference for the presentation of the research worker’s results 
in a manner directly of use to the engineer. 

The reviewer has detected but few misprints and errors in statements, although some must be almost 
inevitable in a work of this kind. Kelvin’s college is stated to be St. Peter’s; the present name is Peter- 
house. 

H. G. Horxrns 
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Anniversary volume on applied mechanics dedicated to C. B. Biezeno by some of his 
friends and former students on the occasion of his sixty-fifth birthday, March 2, 1953. 
N. V. De Technische Uitgeverij H. Stam, Haarlem, 1953. 328 pp. and 2 fold-in 
plates. $5.60. 


The beautifully made-up volume contains a brief biography and a list of publications of C. B. Biezeno 
and the following articles: W. Boomstra, Triangles équilatéres inscrits dans une conique donnée; H. 
Bremekamp, Sur Ja théorie de Sturm-Liouville; Th. vy. Karm4n and G. Millan, The thermal theory of 
constant pressure deflagration; J. M. Burgers, Some remarks on detonation and deflagration problems 
in gases; C. Koning, Some interference problems; R. V. Southwell and Gillian Vaisey,.A problem sug- 
gested by Saint-Venant’s Mémoire sur la torsion des prismes; R. Grammel, Nichtlineare Schwingungen 
mit unendlich vielen Freiheitsgraden; R. J. Legger, The d’Alembert principle; A. D. de Pater, La stabilité 
d’un dicéne se déplacant sur une voie en alignement droit; A. van der Neut, The local instability of 
compression members built up from flat plates; J. A. Haringx, Stresses in corrugated diaphragms; 
A. van Wijngaarden, Ut tensio sic vis; D. Dresden, Shrink-fit used to transmit a torque; W. T. Koiter, 
On partially plastic thick-walled tubes; G. G. J. Vreedenburgh and O. Stokman, Some new elements in 
the calculation of flat slab floors; J. P. Mazure, Statical problems in the code of practice for stee! windows; 
J. J. Koch, The Laboratory for Applied Mechanics at the Technological University of Delft; R. G. Boiten, 


The design of diaphragms for pressure measuring devices based on the use of wire-electrical straingauges. 
W. PRAGER 


The dynamics and thermodynamics of compressible fluid flow. Volume I. By Ascher H. 
Shapiro. The Ronald Press Co., New York, 1953. xiii + 647 pp. $16.00. 


Professor Ascher H. Shapiro of the Massachusetts Institute of Technology has undertaken a very 
important and useful task in writing a treatise on compressible fluids. As the title suggests, it is impossible 
in such a study to separate the dynamical from the thermodynamical point of view and this accounts for 
both the difficulty and the interest of this branch of mechanics of continua. The subject dealt with in this 
book is of great practical importance. At the present time, aeronautical engineers, mechanical engineers, 
chemical engineers, applied physicists and applied mathematicians frequently neéd to investigate and to 
apply the theory of compressible fluids. They will find in this book material of much interest. The author, 
a professor of mechanical engineering, has given a number of practical applications to various branches 
of engineering science in each chapter. He has planned this treatise in two volumes, the whole divided into 
8 parts and 28 chapters. So far oniy the first volume containing the first four parts and part of the fifth 
has been published. The size of the volume (647 pp.) gives some indication of the amount of information 
contained in it. 

Part I sets forth the basic concepts and principles from which the remainder of the book proceeds. 
The first chapter—foundations of fluid dynamics—contains definitions of a fluid and of the continuum 
properties (density, pressure, viscosity) and discussions of the conservation of mass and momentum 
theorems applied to a control volume. In the second chapter—foundations of thermodynamics—the 
author reviews the most important definitions and concepts of thermodynamics and states the two laws. 
The second law is introduced by means of the definition of the thermal efficiency of heat engines and is 
related to the impossibility of obtaining positive work from a system passing through a complete cycle 
while exchanging heat with only a single source. The thermodynamic properties of the continuum are 
then defined (internal energy, enthalpy, etc.) and the theorem of conservation of energy is obtained by 
application of the first law to a control volume. The chapter ends with the basic formulae for a perfect 
gas. In the third chapter which concludes Part I are introduced the principal ideas concerning compressi- 
ble fluids. The speed of sound is defined by its physical significance (velocity of propagation of a plane 
pressure pulse) and the physical differences between subsonic and supersonic flows are then emphasized 
by statement of the three Karman rules for supersonic flow and are discussed with the help of various 
schlieren and interferometer photographs. This chapter also includes a study of the similarity parameters 
and a survey of the optical methods which were used in the experiments (interferometer, shadowgraph 


and schlieren methods). 
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Part II is devoted to “‘one-dimensional flow’’; although it is restricted to the steady case (the un- 
steady case will be treated in Volume II) it is almost 200 pages long. In the first chapter of this part the 
fluid is assumed to be continuous and isentropic, and a complete discussion for converging or converging- 
diverging nozzles and supersonic diffusers is given, including a special study of the choking effect. In the 
next chapter, the theory of normal shock waves and the basic formulae relating to this theory are found. 
Special attention is given to the physical explanation of the formation of shock waves, and some discussion 
is given of the thickness of these shocks. The ideas developed in this chapter are applied to nozzles, super- 
sonic diffusers and the supersonic Pitot tube. Then a chapter is devoted to the detailed study of flow in 
constant-area ducts with friction. It is shown how friction can produce choking even though the duct area 
remains constant. Information on isothermal flow in long ducts is also given. The following chapter is 
devoted to flow in ducts with heating or cooling, i.e., to the case in which, though the area is kept constant 
and the friction is neglected, the stagnation enthalpy may be changed for various reasons. There exist 
many problems in engineering in which such a situation arises. Finally, in the last chapter of Part IT, 
more general problems are discussed in which two or more of the previously mentioned phenomena may 
appear. 

Part III contains an introduction to flow in two and three dimensions. It gives the basic definitions 
and theorems: circulation, rotation and their physical significance; Euler’s equation; Kelvin’s and Ber- 
noulli’s theorems; continuity equation; equation of velocity potential; etc. 

Part IV deals with subsonic flow. The author begins with the so-called linearized theory in which 
small perturbations are superposed on a uniform flow. After the classical example of flow past a wave- 
shaped wall, the fundamental similarity rules are given. These allow us to relate the subsonic compressible 
flow past a certain profile to the incompressible flow past a second profile derived from the first by means 
of an affine transformation. The famous Prandtl-Glauert rule is discussed and many applications are 
made, in particular to wind tunnel corrections. The next chapter is devoted to the hodograph method for 
two-dimensional subsonic flow, which in fact allows us to linearize the problem without introducing new 
simplifications. The tangent-gas approximation of Chaplygin is discussed and applications lead to the 
KArman pressure correction formula. In the following chapter miscellaneous methods and results are 
given—method of expansion in series in terms of the Mach number (Rayleigh-Janzen), method of expan- 
sion in series in terms of a shape parameter, relaxation method, etc. Comparisons of theoretical results 
with experiments are made. The last chapter of Part IV is devoted to three-dimensional subsonic flow— 
flow past spheres, ellipsoids, bodies of revolution, wings of finite span and sweptback wings. The extension 
of the similarity rules is made with special care. 

Part V deals with two-dimensional supersonic flow (only the first three chapters of this part are in- 
cluded in Volume I). The linearized theory is treated in great detail in order to facilitate its application 
to practical problems and to provide a good introduction to the following chapter which is devoted to the 
method of characteristics. This last method is presented as a natural generalization of the basic results 
of linearized theory. Application of the method is made to the jet problem and to the design of a super- 
sonic wind tunnel; the case in which rotation is present is also considered. In the last chapter, the theory 
of oblique shock waves is developed (shock equations, reflection and interaction of shocks). Many exam- 
ples of flow containing shocks are found, together with experimental results. The book ends with an ap- 
pendix devoted to a mathematical formulation of the theory of characteristics and with some useful 
tables (standard atmosphere, isentropic flow, normal shock, one-dimensional flow with friction or with 
change in stagnation enthalpy, hodograph characteristic functions for supersonic flow). 

Each chapter in the book begins with a clear introduction advising the reader of the aim of the study; 
it is followed by a nomenclature list explaining the significance of the symbols used in the chapter. 
After each article, “working formulae and charts’ are given which can be consulted when a specific 
application is planned. At the end of each chapter there are several well selected problems most of which 
are drawn from practical situations arising in engineering science. A selected bibliography is appended 
to each chapter. 

On the whole, the material contained in this book is well selected and clearly presented. The second 
part which is devoted to one-dimensional flow is particularly well developed. Because the book is written 
primarily for engineers, the viewpoint of physics and engineering predominates, emphasis being placed 
on results rather than on proofs. For the most part, the author attempts to make each chapter a self- 
contained unit with the result that some formulae or theorems are proved more than once. The order is 
not always the most logical. For example, the equation of continuity for two- and three-dimensional flows 
is given on p. 283—after Euler’s equation, Bernoulli’s theorem, Kelvin’s theorem and Crocco’s theorem; 
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the reviewer suspects that this equation may have been used implicitly before its formal presentation. 
Moreover, in discussing the linearized theory, the linearization of the boundary conditions does not appear 
to be carried out systematically. Also, it is not immediately clear why the author claims that the first of 
the three similarity laws for two-dimensional subsonic flow given on p. 323 is the most exact. However, 
only a few such criticisms can be made. 

The general level of the book is perhaps too high for it to be of value as a classical textbook but those 
who teach the subject of compressible fluid flow will undoubtedly find it immensely useful. Some years 
ago, Courant and Friedrich’s classical book provided a mathematical approach to the theory of com- 
pressible fluids; this new book provides a more physical and technical approach to the same subject. 
Together, these two extensive treatises constitute a firm basis for the study of this subject. 

P. GERMAIN 


Methods of theoretical physics. By Philip M. Morse and Herman Feshbach. Volumes I 
and II. McGraw-Hill Book Company, Inc., New York, Toronto, London, 1953. 
xxii + 997 pp. (Vol. I), xviii + 980 pp. (Vol. II). $30.00 set. 


In more than 1900 pages and two volumes the authors have packed a formidable amount of material 
on the subject of mathematical methods in physics. This treatise is largely concerned with fields and 
boundary value problems of all types. The discussion of techniques in connection with field problems 
naturally leads to the newer aspects of perturbation theory and variational methods; the reader’s atten- 
tion is called particularly to Chapters 9 and 11 in this connection (see topics listed below). The section 
on Green’s functions is really excellent as are many other sections of this work. Anyone with strong interest 
in the use of advanced mathematical methods will find this work extremely valuable. One of its advan- 
tages is that it is rather complete and self-contained. 

The advantages mentioned become disadvantages when one considers using this treatise as a text. 
On the other hand the work was the outgrowth of a course given for quite a few years by the authors at 
|, a 

The chapter headings are listed below together with section headings within each chapter and the 
space devoted to each chapter. Chapter 1 (117 pages) Types of Fields: Scalar Fields, Vector Fields, 
Curvilinear Coordinates, The Differential Operator Del, Vector and Tensor Formalism, Dyadics and 
Other Vector Operators, The Lorentz Transformation, Four-Vectors, Spinors. Chapter 2 (155 pages) 
Equations Governing Fields: The Flexible String, Waves in an Elastic Medium, Motion of Fluids, Dif- 
fusion and Other Percolative Fluid Motion, The Electromagnetic Field, Quantum Mechanics. Chapter 3 
(73 pages) Fields and the Variational Principles: The Variational Integral and Euler Equations, Hamil- 
ton’s Principle and Classical Dynamics, Scalar Fields, Vector Fields. Chapter 4 (133 pages) Functions 
of a Complex Variable: Complex Numbers and Variables, Analytic Functions, Derivatives of Analytic 
Functions, Multivalued Functions, Calculus of Residues, Asymptotic Series: Method of Steepest Descent, 
Conformal Mapping, Fourier Integrals. Chapter 5 (164 pages) Ordinary Differential Equations: Separ- 
able Coordinates, General Properties and Series Solutions, Integral Representations, Table of Separable 
Coordinates in Three Dimensions, Second-Order Differential Equations and Their Solutions. Chapter 6 
(115 pages) Boundary Conditions and Figenfunctions: Types of Equations and of Boundary Conditions, 
Difference Equations and Boundary ‘Conditions, Eigenfunctions and Their Use, Table of Useful Eigen- 
functions and Their Properties, Eigenfunctions by the Factorization Method. Chapter 7 (100 pages) 
Green’s Functions: Source Points and Boundary Points, Green’s Functions for Steady Waves, Green’s 
Function for the Scalar Wave Equation, Green’s Function for Diffusion, Green’s Function in Abstract 
Vector Form, Table of Green’s Functions. Chapter 8 (100 pages) Integral Equations: Integral Equations 
of Physics, General Properties of Integral Equations, Solution of Fredholm Equations of the First Kind, 
Solution of Integral Equations of the Second Kind, Fourier Transforms and Integral Equations, Tables 
of Integra] Equations and Their Solutions. 

Votvume II. Chapter 9 (172 pages) Approximate Methods: Perturbation Methods, Boundary Per- 
turbations, Perturbation Methods for Scattering and Diffraction, Variational Methods, Tabulation of 
Approximate Methods. Chapter 10 (158 pages) Solutions of Laplace’s and Poisson’s Equations: Solutions 
in Two Dimensions, Complex Variables and the Two-Dimensional Laplace Equation, Solutions for 
Three Dimensions, Trigonometric and Hyperbolic Functions, Bessel Functions, Legendre Functions. 
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Chapter 11 (252 pages) The Wave Equation: Wave Motion on One Space Dimension, Waves in Two 
Dimensions, Waves in Three Space Dimensions, Integral and Variational Techniques, Cylindrical 
Bessel Functions, Weber Functions, Mathieu Functions, Spherical Bessel Functions, Spherical Bessel 
Functions, Spheroidal Functions, Short Table of Laplace Transforms. Chapter 12 (173 pages) Diffusion, 
Wave Mechanics: Solutions of the Diffusion Equation, Distribution Functions for Diffusion Problems 
Solutions of Schroedinger’s Equation, Jacobi Polynomials, Semi-cylindrical Functions. Chapter 13 (143 
pages) Vector Fields: Vector Boundary Conditions, Eigenfunctions and Green’s Functions, Static and 
Steady-state Solutions, Vector Wave Solutions, Table of Spherical Vector Harmonics. 

With each chapter a set of problems is offered. 

The seventeen tables at the end of Volume II have the following headings: Trigonometric and Hyper- 
bolic Functions; Hyperbolic Tangent of Complex Quantity; Inverse Hyperbolic Tangent of Complex 
Quantity; Logarithmic and Inverse Hyperbolic Functions; Spherical Harmonic Functions; Legendre 
Functions for Large Arguments; Legendre Functions for Imaginary Arguments; Legendre Functions of 
Half-integral Degree; Bessel Functions for Cylindrical Coordinates; Hyperbolic Bessel Functions; 
Bessel Functions for Spherical Coordinates; Legendre Functions for Spherical Coordinates; Amplitudes 
and Phases for Cylindrical Bessel Functions; Amplitudes and Phases for Spherical Bessel Functions; 
Periodic Mathieu Functions; Normalizing Constants for Periodic Mathieu Functions and Limiting 


Values of Radial Mathieu Functions. 
Roun TRUELL 


The aim and structure of physical theory. By Pierre Duhem. Translated from the French 
by Philip P. Wiener. Princeton University Press, Princeton, N. J., 1954. xxii + 344 
pp. $6.00. 


This book is a translation of the original work of the author entitled La Théorie Physique: Son Object, 
Sa Structure (Marcel Riviére et Cie, Paris, 1914) and is the first of Duhem’s writings on the philosophy 
of science to appear in the English language. Actually the French version from which the translation 
has been made is the second edition of a still earlier work, published in 1906. 

Pierre Duhem (1861-1916) was a well known French thermodynamical physicist of the late nine- 
teenth century and made nofable contributions to the application of thermodynamics to physical chem- 
istry in developing the ideas of Gibbs, which indeed he was among the first to popularize in France. 
At the same time he took a profound interest in the history and philosophy of physical science. Though 
somewhat overshadowed in this field by his more brilliant and versatile contemporary, Henri Poincaré, 
he nevertheless established a well justified reputation as an incisive critic of physical methodology along 
with men like Mach, Helmholtz, Clifford and Pearson. 

The methodology of physics has been the subject of serious investigation during the past half century 
and to a certain extent this inevitably serves to date the treatise of Duhem. Nevertheless, it appears 
he still has something helpful to say to his successors and his essentially pragmatic, positivistic viewpoint 
of the deductive character of physical theory will prove very stimulating to the present day reader. 

The author’s treatment is divided into the two main sections reflected in the title. In the first he 
discusses the relation between physica] theory and metaphysics, physics and the ‘‘natural’’ classificatio n 
of experience, the history of the development of theories, and the relation between abstract theories 
and mechanical models. The second part is devoted to an analysis of the structure of a physical theory 
from the standpoints of the notions of quantity and quality, the use of mathematics, the relation of 
theory to experiment, and the choice of hypotheses. There are two appendices, the first of which on 
“The Physics of a Believer’, though not at all a detailed discussion of the relation between physical 
science and religion presents some interesting side lights on this perennial question. 

The clarity and vividness of Duhem’s style have been very ably caught and preserved by the trans- 
lator, and the thoughtful reader will never mistake the author’s meaning. Though in general rather 
judicious in his estimates, he was sufficiently strong-minded to develop a few obsessions which give a 
personal flavor to his work. For example, Duhem firmly believed that English physicists approach the 
building of physical theories in quite different fashion from their Continental colleagues. His arguments 
are persuasive but will scarcely stand up in the light of more recent history. Moreover his fervent opti- 
mism that physical theory progresses by steady evolution to an ideal end, the “natural classification’”’ 
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of experience, apparently failed to foresee the unparalleled creation of new experience in modern physical 
laboratories which is giving the present day theorist nightmares. Of course we can all agree with Duhem 
that physica] theories will continue to be constructed and that some of them will be successful, but it 
seems undesirable to hamper ourselves by any preconceived metaphysical notion of an “ideal end’’. 
This book can be recommended in unqualified terms to all who are interested in the logical structure 


of physical science. 
R. B. Linpsay 


Applied elasticity. By Chi-Teh Wang. McGraw-Hill Book Co., New York, Toronto, 
London, 1953. ix + 357 pp. $8.00. 


This elasticity text is clearly intended for the use of the advanced engineering student, as is stated 
in the preface. The classical theory of elasticity is well presented in a manner that presupposes no ex- 
tensive knowledge of advanced calculus. Where necessary, an attempt is made to introduce whatever 
higher mathematics is needed. For the most part, this attempt is successful so far as both rigor and tech- 
nique of solution is concerned, keeping in mind that the main purpose is the latter. However, the intro- 
duction to the theory of complex variables is, I believe, much too hurried and sketchy; and at some 
points certainly does not possess the rigor the author hopes is not sacrificed by concentration on methods 
of solution. 

Chapters I and II develop the concepts of stress and strain, respectively, in straightforward fashion. 
Chapter III is concerned with stress-strain relations and the concept of strain energy. Chapter IV sets 
up the plane-stress and plain-strain problem in terms of Airy’s stress function and solves problems by 
choosing solutions to the biharmonic equation. Rotating disks and thermal stresses are also discussed. 
Chapter V treats the torsion problem in terms of the warping function. The stress function is introduced 
to make use of the membrane analogy in the solution of the torsion problem of thin open sections and 
thin tubes. 

Since one of the laudable purposes of the book is to provide techniques of solution where exact 
analytic solutions are intractable, much space is devoted to numerical methods. Chapter VI discusses 
finite-difference approximations and relaxation methods, with application to the torsion problem as an 
example. In Chapter VII, energy principles and variational methods are discussed analytically. The 
results of Chapters VI and VII are used to obtain numerical solutions of buckling problems in Chapter 
X, the analytic formulation and solution of which are discussed in Chapter IX. 

Chapter VIII uses complex variables to obtain solutions to torsion problems and plane problems. 

Chapters XI and XII discuss bending and buckling of thin plates, and the theory of thin shells, 
respectively. The latter is preceded by a comparatively extensive discussion of the differential geometery 
of a surface. 

The book is clearly written, easy to read, and on the whole does, I believe, fulfill the author’s purpose 
“... to provide the student with the necessary fundamental knowledge of the theory ... to acquaint 


him with the most useful analytical and numerical methods (of solution).” 
Harry J. Wetss 


Stability theory of differential equations. By Richard Bellman. McGraw-Hill Book Co., 
New York, Toronto, London, 1953. xiii + 166 pp. $5.50. 


This attractive little book gives an original and useful survey of various aspects of behavior, as 
t — o, of solutions of systems of ordinary differential equations. After an introductory chapter intro- 
ducing the convenient vector-matrix-norm notation for treating such systems, the author defines the 


sense in which he uses the word stability, as follows: 
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“Definition. The solutions of 
dy 
(7) —~ = A(t 
2 (Oy 


are stable with respect to a property P and perturbations B(t) of type T if the solution of 


dz 
dt 
also possess property P. If this is not true, the solutions of (7) are said to be unstable with respect to 
property P under perturbations of type 7'.” 

This definition can be applied to many problems, concerning non-linear (when rephrased) as well 
as linear differential equations. In many cases (e.g., if A(t) is constant), the pattern followed by the 
author consists in showing that the qualitative asymptotic behavior of known special solutions is un- 
affected by ‘‘small”’ perturbations of the coefficients—Liapounoff’s famous stability theorem is a special 
case. In other cases, striking counterexamples to plausible guesses are given. In still others, ingenious 
isolated results are “rescued from oblivion” by displaying them in easily accessible form (e.g., the 
Fowler-Emden equation, to which Chap. VII is devoted). 

The reader who is looking for results of the type just described should first consult Bellman’s book. 
If he does not find them there, he will probably locate them by consulting the well-organized bibliography. 

On the other hand, there are a number of important topics, whose omission is not suggested by the 
somewhat misleading title of this book. For example, the basic Routh-Hurwitz stability criteria are 
nowhere mentioned, being replaced throughout by the phrase “Let all solutions of dy/dt = Ay tend to 
zero as t — .”’ Again, no discussion is given of topological arguments or of the stability of limit cycles 
in the Rayleigh-van der Pol equation. In treating the simple second-order equation (ku’)’ + ru = 0, 
deeper methods (saddle-point method, asymptotic expansions in , etc.) are not explained. 

Perhaps because of his exclusive concern with a single problem and with formal methods, the author 
achieves an admirable clarity and uniformity of style. The result is a very useful survey of the qualitative 
asymptotic theory of ordinary differential equations. 


(8) = (A(t) + Bid)z 


GARRETT BIRKHOFF 


Handbook of elliptic integrals for engineers and physicists. By P. F. Byrd and M. T. 
Friedman. Springer Verlag, Berlin-Gottingen-Heidelberg, 1954. xii + 355 pp. $8.58 
(paperbound), $9.44 (clothbound). 


This useful handbook contains over 3000 integrals and formulas and 28 pages of numerical tables 
designed to aid engineers and physicists in the evaluation of elliptic integrals that occur in practical 
problems. The explanatory material is written on an elementary level and does not require previous 
acquaintance with elliptic integrals or functions. The notations of Legendre and Jacobi are used in pref- 
erence to that of Weierstrass. The typographical presentation, of particular importance in a work of this 
character, is excellent. A considerable saving of space has been achieved by the use of 9-point type in 
numerators and denominators of fractions. The resulting compactness of formulas greatly facilitaties 


the use of the handbook. 
W. PRAGER 


Proceedings of the Eastern Joint Computer Conference. The Institute of Radio Engineers, 
Inc., New York, 1954. 125 pp. $3.00. 
This report contains the Papers and Discussions presented at the Joint I.R.E.-A.I.E.E.-A.C.M. 


Computer Conference held in Washington, D. C., in December 1953. The Conference was the third of 
its kind and the theme chosen for it was, “Information Processing Systems—Reliability and Require- 
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ments.”’ The papers can be classified in three groups: (1) Statements of data processing and calculation 
requirements. Descriptions of requirements for the Life Insurance Business, Numerical Weather Pre- 
diction, Air Traffic Control and Linear Real-Time Systems are given. (2) Papers describing the operational 
and performance characteristics of existing machines. Operating experiences with the OARAC, the Los 
Alamos 701, a large REAC installation, the UNIVAC, the SEAC and the ORDVAC are described. 
Two papers deal with performance tests, the ‘“National Bureau of Standards Performance tests” and 
the “Acceptance Tests for the Raythorn Hurricane Computer.’”’ Two further papers can be placed in 
this group. They are descriptions of the MIT Magnetic-Core Memory and the Electrostatic Memory 
of ILLIAC. (3) The third group of papers is descriptive of work carried out in connection with the 
reliability of components; Magnetic Tape, Electrolytic Capacitors, Resistors, Electron Tubes, and 
other electronic units are discussed. 

Two very interesting papers, one by H. H. Goldstine and the other by J. W. Mauchly, do not con- 
veniently fall into one of the above catagories. Goldstine gives a brief resumé of von Neumann’s work 
on logical design and Mauchly discusses the advantages of built-in checking. 

The report as a whole is well printed and readable. 

J. FouLKEs 


Linear operators spectral theory and some other applications. By Richard G. Cooke. 
MacMillan and Co., Ltd., London, 1953. xii + 454 pp. $10.00. 


The topics presented in this book range over the spectral theory of operators in Hilbert space, an 
introduction to quantum mechanics, a valuable discussion of Banach algebras and certain ideas on 


sequence spaces. 

The Hilbert space theory is introduced by a discussion of ¢ and of abstract Hilbert space. The 
usual elementary operator theory is covered and the spectral theory due to von Neumann using the 
Cayley transform. The Stone theory of spectra and resolvent, a development of the Hallinger theory 


due to Lengyel, the Cooper theory based on integrating the differential equation for the Stone group, 
and the Riesz-Lorch proof are given. 

The presentation tends to be non-geometrical, with a preference for analytical and matrix methods. 
This permits a simple and direct connection with quantum mechanics. The harmonic oscillator is treated 
both in the Heisenberg and Schrédinger formulation and perturbation theory is given also in matrix 
terminology. 

The analytic and non-geometric viewpoint is further developed in Chapter 6, which is entitled 
“Projective convergence and limit in matrix spaces and rings.” The basis of this chapter is certain defini- 
tions of limits for sequences, which are given in another book of the author, i.e., R. G. Cooke, “Infinite 
matrices and sequence spaces,’”’ MacMillan and Co., (1950), and not repeated in the present work. 
A number of topologies for sequence spaces are introduced and these in turn lead to topologies for the 
matrices associated with linear transformations between sequence spaces. The work of H. 8. Allen is 
compared with that of Kéthe and Toeplitz and Banach’s results, specialized to sequence spaces. Greater 
generality is claimed for Allen’s notions, but they apply, of course, only to the sequential situation. 
However, the discussion of Chapter 6 seems to the reviewer to be of a far more specialized character 
than that of the other chapters and of less general interest. 

Chapter 7 is on Banach algebras and is quite fascinating. The initial sections contain Zron’s Lemma, 
the Hahn Banach Theorem, a discussion of weak convergence of linear functionals and an algebraic 
discussion of maximal ideals. The objective is the development of the Gelfand theory in the case of 
commutative Banach algebras. This theory is applied to yield certain results of Wiener’s on the conver- 
gence of the Fourier series and the Fourier integral for a reciprocal of a function. It is also used as a tool 
in the investigations of topological spaces by means of the set of continuous functions defined on them. 
This chapter also contains a discussion of Radon measure and a proof of Wiener’s theorem on the density 
of the translations of a function in L. 

It should be clear from the above that this represents on the whole an interesting and useful collection 
of topics in the mathematical theory of operators. Most graduate students should be able to read the 


book without difficulty. A useful set of references for further reading is also given. 
F, J. Murray 





























